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Abstract 



The ground state energy of a quantum field in the background of clas- 
sical field configurations is considered. The subject of the ground state 
energy in framework of the quantum field theory is explained. 

The short review of calculation methods (generalized zeta function and 
heat kernel expansion ) and its mathematical foundations is given. 

We use the zeta-functional regularization and express the ground state 
energy as an integral involving the Jost function of a two dimensional scat- 
tering problem. We perform the renormalization by subtracting the contri- 
butions from first several heat kernel coefficients. The ground state energy 
is presented as a convergent expression STiitcd for numerical evaluation. 
The investigation for three models has been carried out: scalar quantum 
field on the background of scalar string with rectangular shape, spinor vac- 
uum polarized by magnetic string of the similar shape and spinor vacuum 
interacting with the Nielsen Olesen vortex. Using the uniform asymptotic 
expansion of the special functions entering the Jost function we are also 
able to calculate higher order heat kernel coefficients. 

Several features of vacuum energy have been investigated numerically. 
We discuss corresponding numerical results. 
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1 Introduction 



1.1 General remarks 

The vacuum energy is a spectacular manifestation of the quantum nature of 
the real objects, in particular the fields of matter. In the framework of natural 
science it is also a powerful tool for investigation and of fundamental properties 
of our world. 

It is a really non-trivial fact, that quantum field can be observable even in the 
case that there are no real particles of it present. The reason for this phenomena 
is one of the main properties inherent basically to any quantum system, namely 
the existence of the ground state. Indeed, it is well known from such a simple 
example of a quantum mechanical system as the linear harmonic oscillator with 
the eigenfrequency tu, that the energy spectrum of this system possesses the lowest 
state and the energy of this state E'^ is not zero but is proportional to the Planck 
constant 



This fact is generally also an intrinsic feature of any quantum object. Here we 
mean first of all quantum fields, namely matter fields (fermionic fields) and fields 
of interaction (bosonic fields). It means, one can establish the existence of a 
quantum field in the nature without to observe any real particles (quantums) 
of this field, but only through non zero vacuum values of observables, that are 
measurable, for example the vacuum energy. 

The main point of this phenomena is that the energy of the ground state of a 
quantum system is not defined absolutely, that applies by the way in general to the 
energy itself as a physical feature. It is defined generally with a relation to some 
other state. That is why the discussion about the vacuum energy of quantum 
field makes sense in correspondence to some external condition, under which it 
becomes visible (or measurable in physical sense). Such an external condition 
can be some boundary or some classical field interacting to the quantum one 
(background) 

The most famous manifestation of the vacuum energy was the Casimir ef- 
fect (H.B.G.Casimir 1948), which consist in an interaction between two parallel 
conducting plates playing the role of boundaries for the electromagnetic field 
in vacuum state. This effect was predicted theoretically in jJEI and measured 
experimentally 10 years later. 

Sometimes, all the problems concerning the calculation of vacuum energy 
under external conditions are being called "Casimir energy", including also the 
problems of vacuum polarization by classical backgrounds. Actually, the vacuum 
energy of quantum spinor field polarized by a special homogeneous magnetic 
background was calculated by W.Heisenberg and H.Euler in 1936, 11 years be- 
fore the discovery of the Casimir force. 



E° = -huj. 
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1.2 Motivation 



The investigation of vacuum effects in general and vacuum energy in particular 
has a number of applications to wide variety of physical problems (see recent 
reviews fl^ I65| 1^)- Indeed, in a lot of problems concerning the experiments in 
particle physics, certain effects caused by interaction with vacuum fluctuations 
must be taken into account as radiation corrections. For example in the bag 
model of hadrons in Quantum Chromodynamics the vacuum energy of quark 
and gluon fields makes essential contributions in the total nucleon energy. As 
further examples of such measurements one can refer to the experiments for the 
exact mass, isotopic charge, anomalous magnetic moment and so on. 

Generally speaking, the vacuum corrections are more or less relevant in all 
problems of particle and high energy physics. 

Moreover, one can imagine a situation that only the sole possibility to establish 
the existence of a quantum field in the nature, (since it is present in the model) is 
the interaction of its vacuum with an observable classical background (and mea- 
surable effects reasoned thereby). Actually, the interaction of the classical field 
configuration with quantum fields being in vacuum state leads to modification 
of macroscopic physical features of this classical object. In particular the energy 
and its dependence on the external parameter of this configuration undergoes 
an alteration. This interplay between the quantum vacuum and classical field 
is also called sometimes "Casimir effect". The modification of classical energy 
by vacuum corrections for real objects is usually negligible small. However, one 
can imagine some special field configurations such that the quantum corrections 
become comparable with the classical energy of the background. 

Therefore it is a reasonable question to ask: could the quantum fluctuations 
change the main physical properties of classical object? First of all, the question 
concerns the stability (or instability) of field configurations. 

For example, one can look for a field configuration, that acquires stability 
through an interaction with the quantum vacuum, although it would be instable 
in pure classical approach. 

In particular, the issue of stability of strings is of interest. In electroweak 
theories, and in QED in particular, the coupling is small. Hence the vacuum 
energy (being a one loop correction to the classical energy) is suppressed by 
this coupling and indeed in most of the cases is not sufficient to change the 
stability properties. While in QED a magnetic string is intrinsically unstable, in 
electroweak theory there are unstable and stable configuration. Here quantum 
corrections may become important for the stability. A question of special interest 
is whether a strong or singular background may have a quantum vacuum energy 
comparable to the classical one. Such an example is considered in the present 
work for two kinds of field. 
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The vacuum correction could provide stability to classical instable objects, but 
vice versa it could also destabilize classical stable objects. The good illustration 
for this speculation can be provided by the model with the vortex of the Nielsen- 
Olesen type. This configuration is known to be stable for any integer vorticity n, 
if the parameter jS < 1. If the (3 > 1, then only the vortex with n = 1 remains 
stable, all configurations with the vorticity n > 1 decay into n simple vortices 
with n = 1. The question arises there: could the quantum vacuum fluctuation 
modify the parameter f3 in such a way that transition from stable configuration 
to instable one takes place? The corresponding model will be also considered in 
the following sections. 

The results and techniques gained here can be also applied to the problems 
of cosmology and to Condensed Matter Physics. 

In Mathematical Physics the investigation of vacuum polarization has stimu- 
lated the development of powerful regularization and renormalization techniques 
based on the use of zeta function and heat kernel expansion, which are also con- 
sidered below. 

The subject of the work is restricted on the investigation of ground state 
energy of vacuum fields in presence of classical static field configurations with 
cylindrical and translational symmetry. 



1.3 Historical review 

Some people claim, that the subject of vacuum energy has no sense since this 
value is originally infinite. This point of view is in principal true, but the fact 
of Casimir force between two plates suggests the simple thought: if two plates 
are positioned with the distance d from each other, the photon vacuum possesses 
really an infinite energy, and for the distance d' the vacuum energy is infinite as 
well, but not the same, and crucial for the problem is that the difference between 
these two values is finite and moreover measurable. 

The same concerns the vacuum energy of quantum field, polarized by some 
external background. The first known example was the calculation of the vacuum 
energy of a spinor field in the background of a magnetic field. This problem 
was first considered as early as in 1936 by W. Heisenberg and H. Euler |53j . 
They were interested in the effective action in the background of a homogeneous 
magnetic field. It turned out that the global energy of spinor vacuum in presence 
of constant homogeneous magnetic field (the whole space is filled uniform by 
the field of magnitude B) depends on this magnitude non trivialy. In particular 
the characteristic behaviour for the vacuum energy in the asymptotics of strong 
potential 5 — >• oo in the form 



lO 



BHnB 




was obtained. 
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Since the classical work j^Hl of H.B.G. Casimir, where the energy E'^ of vacuum 
polarized by two conducting planes was calculated analytically with the finite 
result 



E 



lO 



(1.3) 



720 a3 



a number of similar problems was investigated for various external conditions 
(boundary conditions, background potentials etc.). One may consult, for exam- 
ple, in the recent review j23j or the books • 

The discovery of the Aharonov-Bohm effect jH] gave rise to a series of works 
researching various approaches to the vacuum in presence of configurations similar 
to the Aharonov-Bohm string, so the vacuum in gauge theories claims its non 
trivial geometrical and topological structure, that is closely connected with the 
structure of gauge groups. 

In this way local properties of photon and electron vacuum have been actively 
investigated. There it has been shown that the structure of vacuum state pos- 
sesses a more reach variety of features as merely the global ground state energy 
jlHj- In particular the subjects of interest were various local effects such as dis- 
tribution of the vacuum energy density [73], vacuum currents |3S], |SH] and the 
backreaction (so called "feedback effects") |35] . 

The other important application of the external field problem is the problem 
of non stationary vacuum under external conditions. This activity includes all the 
effects of particle creation by time dependent boundaries 81J, or non conservative 
backgrounds, e.g., electric field. The most concerned problems were investigated 
in early 70's, like |66|, the results were restricted and confirmed for reduced space 
dimension and generalized for arbitrary high dimension This interest 
was also motivated to explain of Hawking radiation of black holes. 

The class of investigated backgrounds is usually restricted by the complication 
and solvability of a certain problem. In applications of bag models |37|j, ^27|j or in 
cosmology the spherical symmetric backgrounds are considered ^Tj |H3 1^ , for 
string like models the cylindrical symmetry is present [HZl 01201 • 

Pure analytical calculations are only possible for simplest configurations, e.g. 
for the considered classical sample of homogeneous field [Sni;|Sn] or a very specific 
ones such as the transparent refiectionless domain wall [21]. For more complicated 
cases a number of alternative numerical approaches has been developed |68| [33 




1^ . In the present work we restrict ourselves on the problems which can be 
treated analytically, (or analytically to some degree, demanding finally numerical 
evaluations [H [201 IZH 123 (Zj) 

Although a variety of problems for various backgrounds and in many QFT 
models have been successfully investigated, and a corresponding experience in 
handling these problems have been gained, the general understanding of the prob- 
lem is so far not complete, and many features are still unclear. 
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The main open question is that no general rule for the dependence of the 
vacuum energy on the background properties has been found so far. In particular, 
it is unknown how to forecast the sign of the energy. Some special cases only 
has been recently investigated. An interesting approach is that used e.g. in [10] 
where the issue of the sign and of bounds on fermionic determinants in a magnetic 
background had been considered. 

1.4 Examlpes 

In QED, for the background of a flux tube with constant magnetic field inside the 
ground state energy was calculated in [20]. It turned out to be negative remaining 
of course much smaller than the classical energy of the background field. 

The remarkably simple case of a magnetic background field concentrated on a 
cylindrical shell (i.e., a delta function shaped profile) was considered in Here 
different signs of the vacuum energy turned out to be possible in dependence on 
the parameters (radius and flux). However one has difficulties with a physical 
interpretation of this result since the model possesses an infinite classical energy 
density. So one needs a consequent interpretation in a suitable limit, as for exam- 
ple was done for an infinite thin Aharonov-Bohm string ^421 ^ . Here a somewhat 
intermediate model between the rectangular potential and 5-shaped potential 
[7^ is investigated in the present work. 

An interesting approach is that used in jJU] where the issue of the sign and of 
bounds on fermionic determinants in a magnetic background had been considered. 

Presently, the interest is shifted to string like configurations. In [55 the 
contribution of the fermionic ground state energy to the stability of electroweak 
strings was addressed. In [30! the gluonic ground state energy in the background 
of a centrum-of-group vortex in QCD has been considered. 

The ground state energy for vortex type solutions in gauge models with min- 
imal coupling has been calculated implementing special numerical and partially 
analytical approaches to U(l) [HZI and more complicated SU(2) and SU(3) theo- 
ries pni- In the case of supersymmetry the sufficient simplification is established 
through the high symmetry, that allowed to calculate recently a non-zero vacuum 
correction to the mass [7^] and to the central charge [nH] of the SUSY N=2 vortex. 
In the present work the result of |T^ are extended to the vortex of Nielsen-Olesen 
type in the simplest U(l) abelian Higgs model. 

1.5 The aim of the work 

The main goals of the present research are: 

Methodologically, an improvement of the existent technical tools to calculate 
the ground state energy in classical backgrounds is of interest as well as the 
applicability of these technics for more complicated configurations. 



10 



Physically, the influence of the quantum corrections on the background, in 
particular on the stability properties of the fleld system at all, and the asymptotic 
behaviour of the one-loop vacuum contribution to the total energy, are to be 
investigated by using the applied method. 

1.6 The method 

The class of considered problems to which we are restricted in the present work 
excludes basically all the problems with non stationary vacuum, since the method 
used below is not developed for such kind of models. 

All the considered models contain one quantum fleld in vacuum state polar- 
ized by a classical background with special symmetry properties that simplify 
sufficiently the calculation procedure. 

The expression for the ground state energy contains the ultraviolet diver- 
gences. To make the results of this calculations physically meaningful one has to 
extract the finite value from this infinite result, in other words it is necessary to 
subtract the ultraviolet divergences using some renormalization prescription. 

There are various methods developed to get rid of non-physical divergences. In 
the present paper the zeta function regularization procedure has been used. This 
approach is well known and we follow in general Roughly speaking one has to 
subtract the contribution of the first few heat kernel coefficients (oq through 02 in 
the given case of (3+1) dimensions). After that one can remove the intermediate 
regularization and one is left with finite expressions. However, in order to obtain 
these finite expressions in a form suitable for numerical calculations one has to go 
one step further. As described in [2011111, see also below in Sec l5.4^ one has to add 
and to subtract the certain part of the asymptotic expansion of the integrand. 

In the present work we generalize the analysis done in [201 of rectangu- 
lar shaped background fields (scalar and magnetic) and, furthermore, for mixed 
background given in the from of numerical distribution. 

The proceeding has been arranged consequently through the sections 3-6 
simultaneously for the following three considered models. 

We consider the vacuum energy of a spinor field in QED for a rectangular 
background, the scalar background with the same shape interacting to the scalar 
fiuctuations, and finally the spinor vacuum polarized by the background config- 
uration in form of the Nielsen-Olesen vortex in the abelian Higgs model. 

The problem of the rectangular magnetic background reproduces in principle 
the main features of the one considered in [20j5 it is however technically more 
involved and allows progress in two directions. First, it allows to refine the 
mathematical and numerical tools for such problems and, second, it allows to 
address the question how the vacuum energy behaves for an increasingly singular 
background (making the rectangle narrower). So this model interpolates to some 
extend between the fiux tube with homogeneous field inside in |20] and the delta 
shaped one in [721 El • 
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The basic principles of this analytical calculation procedure are described in 
sections 3-5. We start with the well known zeta functional regularization. The 
regularized ground state energy is represented as a zeta function of a hamiltonian 
spectrum and treated in terms of the heat kernel expansion. A representation of 
the regularized ground state energy as an integral of the logarithmic derivative 
of the Jost function for wave scattering problem on the external background is 
obtained. The divergent part is identified as that part of the corresponding heat 
kernel expansion which does not vanish for large m (mass of the quantum field). 
After the subtraction of this divergent part the remaining analytical expression 
must be transformed in order to lift the regularization . A part of the uniform 
asymptotic expansion of the Jost function is used for this procedure. The Jost 
function is obtained from the exact solutions of the equations of motion (these 
are the Klein-Gordon and Dirac equations respectively), derived in SecOl 

The explicit form of exact and asymptotic Jost functions is considered in 
Sec J5.2l and 15.31 In Sec 15. 11 15.41 the representation useful for further numerical 
evaluations is derived. Sec ]5.6l is devoted to the calculation of the heat kernel 
coefficient for the special magnetic background. 

Finally the analytical expressions for the ground state energy are evaluated 
numerically. SecEl contains some numerical evaluations of the finite ground state 
energy for the considered models. 

The ground state energy for the Nielsen-Olesen vortex background cannot 
be performed analytically since the background itself is given as the numerical 
solution of Nielsen-Olesen equations . In this case some analytical manipulations 
applied for the rectangular background have been replaced by the corresponding 
numerical procedures. 

Additionally, the paper contains an introduction to the subject of ground 
state energy of quantum field in Sec ]2.2[ Some technical details are contained in 
the Appendix. 
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2 The ground state energy of quantized field 



2.1 The effective action 

In the generalized approach the contribution of the one-loop ground state 
energy appears as a first correction to the classical action in the so called "effec- 
tive" action. In order to derive this representation we start with the basic object 
of the quantum field theory in the functional integration approach namely the 
generating functional. 

Let the action S of the considered system consist of the classical field $ and 
let the quantum field be defined in the canonical way as 

S[<l>, 0] = S^lassm + ljdx< <P{x)Km<P{x) >, (2.1) 

where $ is a classical field configuration called up from here and below "back- 
ground", and is a quantized field, so called "fluctuation". The symbol K[^] 
relates to the kernel of the corresponding equation of motion, (the so called kernel 
of free action). It means particulary K[^] = □ + + V[^] for scalar field with 
scalar potential V and = ij'^^df^ — ieA^) —m for spinor field with the gauge 

background A^. The edge brackets denote the full scalar product on the space 
of 0. Proceeding from the defined action S'[$,0] (j2.ip . we construct the related 
generating functional for Green functions of 

Z[$,J] = J P0exp{i5[$,0] +i j (l){x)J{x)dx}, (2.2) 

where the explicit dependence on the fluctuation does not appear, since it is 
integrated out. Here we have to notice, that we work in this paper with abelian 
background fields only. For non-abelian models the formula ()2.2|) is in general no 
more valid j43j. 

If we restrict the problem only to connected Green functions, the correspond- 
ing generating functional iy[<l', J] is related to the J] by 

Z[(^,J]=exp{tW[(^,J]}. (2.3) 

One can formally restore the function 0(x) from the W[^, J] by the functional 
derivation 

and express the source parameter J{x) with respect to the field variables 

J = J[$(x),0(x)]. (2.5) 
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This change of variables corresponds to the Legendre transformation {$(x), J{x)} 
{^{x),(j){x)}. The effective action r[$(x),0(x)] is defined in terms of these new 
variables as 

r[$, (J)] = ly [$, J[<l>, </)]] - J dx< ^(x) J[$(a;), (pix)] > . (2.6) 
Using the Gaussian integral for the functional integrals of kind 

'^\e-«(^) (2.7) 



/27r 

for any quadratic form Q(x) = | < x, Kx > + <b,x> + cso that 

a^^-ii<.,K.>+<b,x>+c) ^ g^p|i ^ ^-1^ > -c}(det ir)-^ (2.8) 
v2vr 2 

for the functional integrals ()2.2j) . one can calculate the generating functional 
Z[$(a;), J(x)] explicit: 

Z[^,J] = e*^='-[*](detir[<l>])^5e-5<-^^~^-^>. (2.9) 

Here the det K[^] is the functional determinant of the kernel of the free action 
(e.g. □ — m^, more examples for certain models are listed below, (SecOl). 

The choice of the sign of the power (det i^'[$])^2 corresponds to the bosonic/ 
fermionic fields (p respectively. Using the known representation of the logarithm 
of determinant through the trace of logarithms 

Indeti^' = Tr InK, (2.10) 

one can arrive for the connected generating functional W^[$(x), J{x)] at 

Wmx),J{x)] = S,iassm±'-Tr \nK-^<JK-'j> (2.11) 
and for the effective action ()2.6|) 

Tmx),<j){x)] = S^iassm ±'-Tr \nK + ^< <f)K~'<j) > . (2.12) 

If the quantum field (j) is considered only in the ground state, the corresponding 
effective action reads 



r[$(x),0(x)]U=o = r[<l'(a;)], (2.13) 

whereat it is supposed that there is no spontaneously broken symmetry of the 
ground state. 

The zeta-functional regularization formalism can be applied to the effective 
action in the form ()2.13|) . The corresponding zeta function - regularized effective 
action is considered in the (Sec. Ej) 
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2.2 The ground state energy in the framework 
of canonical quantization 

The explicit expression for calculation of the ground state energy can be derived 
also from the canonical quantization formalism, that we even use. In the following 
work two kinds of quantum field are discussed-the Klein-Gordon scalar field and 
the Dirac spinor field. 

2.2.1 Quantized scalar field 

Consider a simplest case of the real massive scalar field with the lagrangian 

L = \{d^^d^^-m\''} (2.14) 

leading to the Klein-Gordon equation 

{□ +m2}(^(f,t) = (2.15) 

All the further discussions and calculations are concerned with the Minkowski 
space-time g^y = diag{l, —1,-1,-1}. Then the ()2.15|1 will be rewritten as 



- A + m^}ip{x,t) =0. (2.16) 



Since the equation is of second order with respect to time t = we can define 
so called "positive-" and " negative- frequency" solutions ip~^{x,x'^) and (p~{x,x^). 

The further explanation follows the textbook ^1]. Let the equation of motion 
have a continuous spectrum of solutions ip^{x,x^). The (p{x) are representable 
as the Fourier integral: 

oo 

^(x) = — i— / 4^ [e^-^'^^^'^e^^^^^^Vl^cfc) + e(^'=°^")e-^(*^^^V+(fco,^) 



(27r)3/2 J L 

— oo 

(2.17) 

The normalization of the solutions is kept to be the one used in JT], but 
for our problems it is unimportant. The sign switch of the frequency ^(/cq) is 
supposed to be contained in the corresponding Fourier images (p^{ko, k). 

We consider the problems with the translational invariance along the z-axis. 
It means, the z-variable is not present in the equation of motion explicitly. The 
solution of the equation ()2.16|) may contain only the s z-dependent factor of kind 
g±(ipz2) Q^-^^ ^]^jg factor can be separated for the next step. 

The energy of quantum field will be understood by means of the 00-component 
of the canonical energy momentum tensor T^^ , that is defined for a given la- 
grangian C as 

dC 

^'^' = ^7^^>-^'''^ (2.18) 
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For the zero- zero component of T^^" (the hamiltonian density) of the neutral 
scalar field we have 

nix) = T°°(x) = ^ {[do^ix)f + [Vipix)f + mV(x)'} (2.19) 

By (I2.19|) we have derived the explicit from of the hamiltonian density (energy 
density per unit volume) of the field (p. Then the resulted function of coordinates 
T^^{x) describes a local spatial distribution of the energy of scalar field ip 

The energy contained in the finite volume V can be obtained through the inte- 
gration over d^x. Suppose we are interesting for the energy of the field contained 
in the layer of finite thickness L between two plates transversal to the axis z, say 
z = and z = L. The global field energy contained in this layer is the integral 



L 



Er = T'^ = l dz d'x 



2 i J J 2y/koPo 



ip+{p)ip-{k)e-'^''^''''-P''^e''-''-P^^ - ip-{p)p'^{k)e'''°^^^^ 

The first and the second summand turn out to be equal zero if the integra- 
tion over d^x is performed, in the remaining two terms the part, containing the 
dependence on z- variable can be separated and the integration over dz builds the 
Fourier representation of delta-function 6{pz — kz), that reads 

oo 

S{x) = [ ^e-'P"" (2.20) 

J 277 

— oo 

Finally one arrives at the form: 



El = jT'' = \L j'^ j d'p Po[^^p)p--{p) + ^-{p)^-'{p)] (2.21) 

where the remaining components of the momentum space are denoted as d'^p. 
Divided by L this value means the energy density per unit length. 

The canonical secondary quantization procedure (for bosonic scalar field) 
means the decomposition of the arbitrary solution (p{x, t) of ()2.15|) onto the eigen- 
functions 

/rpk 
J^^^Hk)e-''^ + aHk)e^'^^] (2.22) 
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with the coefficients a{k),a^{k) called the annihilation and creation operators 
with the momentum k. The integration is carried out with the invariant measure 

on-shell,[ni[Zn!, cok = VP + m?. 

The commutation relations between a{k), a^{k) implied by the normalization 
conditions in the Fock space |llj are 



[a{k), a+(A;')] = 5{k - k'); [a+(A;), a+{k')] = [a{k), a{k')] = (2.23) 
Within these relation the hamiltonian 7i takes the form 



N{k) + I 



(2.24) 

where A^(A;) = a^(k)a^k) is the operator of number of particles with the momen- 
tum k. 

The vacuum state is understood as the lowest state of the quantum field in 
the Fock space is the vacuum state |0 > defined by 

a(fc)|0>=0, (2.25) 
and the corresponding co-state < 0| respectively 

<0|a+(A;)=0 (2.26) 

Now we define the vacuum energy as a vacuum expectation value of T^^ in 
the vacuum state 

E^ac = I < 0| r°° |0 > . (2.27) 



All problems considered below in this paper concerns with the backgrounds 
possessing translational symmtery along the z-axis, it is integrated over the z- 
component of the momentum, p2/(27r) and the transversal component of the 
momentum, p (say, for cylindrical coordinates p = Pnp), can be quantized, so 
that it acquires discrete values p^n). Then, by means of the integration measure 
of Fourier transformations J8Q , the integration over ^ will be substituted by 
summation over p(^n) . Since we mention the vacuum state as a state without 
real particles, it leads to the relation 

1 f dp, f (Pp f dp, ^ 

= 2 y 27 y (2^^- ^'^t^ 

The ()2.25I2.26|) . substituted into the ()2.21|) . provide the final expression for vac- 
uum energy density S per unit length of the translational invariant field config- 
uration 

(n) 
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where ujn = a/ Pn^ + m? are the eigenvalues of the hamiltonian ()2.19|1 . 

Here it is suitable to notice that the solutions ({>{x) of the Klein-Gordon equa- 
tion should be orthonormalized with respect to the inner product for scalar field, 
but the normalization coefficients of the eigenfunctions (x) does not appear 
in the final expression of the ground state energy because we handle here with 
the eigenvalues only; thus the discussion concerning the normalization coefficients 
can be omitted for our problem. 

Thus we obtained that the vacuum energy defined by 12.271 is equal to the 
half of sum over the spectrum of eigenvalues of hamiltonian, that has been first 
defined as the ground state energy energy. Therefore the both definitions coincide 
with. However, it is necessary to notice that this equivalence is true for a system 
of not interacting ("free") fields. It will be mentioned further in SecEJin context 
of the effective action. 

2.2.2 Quantized Dirac field 

Let the quantum field il){x) now be the 4-component spinor field with the free 
lagrangian 

£ = 1 {^{x)-^^'d^il){x) - d^i){x)-i^'^{x))] - m^{x)il){x), (2.30) 
and the spinor function iIj{x) obeys the Dirac equation 

{i-1^d^ - m)i){x) = (2.31) 

One can use however the simplified lagrangian that also leads to the equation 
dTSU), namely 

C = ij{x) {i^d^, - m} ip{x) (2.32) 
The canonical energy-momentum tensor corresponding to the fl2.30|) reads 

T^"" = ^ {4j{x)Yd''^{x) - d''^{x)Y^{x) } - g'"'C (2.33) 

For simplicity we work here with the short form ()2.32j) . The corresponding 
energy-momentum tensor T'^'^ is 

T'^'^ix) = tlj{x)Ydf,^{x) (2.34) 

and the 00-component of them, that we are interested in, is the corresponding 
hamiltonian density 7i(a;), that is defined as 

nix) = T^\x) = ^(x)7°9oV'(a;) (2.35) 



18 



Let the solutions of the Dirac equation be ip'^ and ip , which are normahzed 

by 

J d^xip^{x)ip{x) = 1, (2.36) 

what will be no further discussed since the normalization of wave functions does 
not play any role in the resulting ground state energy (the remark in the end of 
the previous subsection). 

Now use again the decomposition on the states with positive and negative 
frequences for each of solutions 'ip{x),il!{x) as in the scalar case above, 

^±(x) = J d'ke^^'^ Yl ^iikX^ik) (2.37) 

i^^ix) = J d'ke^^'^ J2 «^(^)^'^^(^) (2-38) 

the index z/ enumerates the two states with different spin projections, v^^lk) and 
v^^{k) arc the 4-componcnt spinor amplitudes, spinor indices arc omitted. 
The resulting hamiltonian density expressed through the amplitudes reads 

yoo = / d^k k^e^'^'' [at{k)a-;{k) - aZ{k)a+{k)] (2.39) 

The secondary quantized spinor field in the canonical formalism will be rep- 
resented as 

'il){x) ^'il)^{x)+'il)-{x), ^{x) ^ilj+{x)+ijj-{x) (2.40) 

where the creation and annihilation operators a^,a~ and a~,a^ obey now the 
anticommutation relations in order to have a positive definite form for the energy. 

{a;(k),a-(k')}^5^J(k-k') (2.41) 

{a-{k),at{k')} = 5^J{k-k') (2.42) 

{a+(fc),a+(fc')} = {d-{k),a-{k')} = (2.43) 

{d+{k),d-{k')} = {a;{k),a-{k')} = (2.44) 

By applying the creation/annihilation operators and complex (but not 
Hermitian) conjugated of them to the vacuum state 

< 0|a+ = a-|0 >= (2.45) 

< 0|a+ = a-|0 >=0, (2.46) 



19 



one arrives at the vacuum expectation value for this hamiltonian 



(n) 



(2.47) 



For translational invariant problem the corresponding energy density per unit 
length is 



Thus the hamilonian density, that involves the summation over 2 states with 
different sign of energy ( fermoin and anti-fermion energies ) distinguished by the 
sign of mass, and over 2 spin projection states, compared with this one for the 
scalar field ()2.29l of the previous subsection) acquires the resulting multiplicative 
factor -4. 

In the following sections the formulae ()2.48l and l2.29|l will be used as a starting 
point for calculations of the ground state energy . 




(2.48) 
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3 The equations of motion and solutions 



We consider basically the field system which consists of quantum fields (fluc- 
tuations) interacting with classical fields $i of fixed geometrical configurations 
(backgrounds). 

The corresponding lagrangian is constructed as 

£ = + + Cint-, (3-1) 

where Cmt is an interaction between cf) and $j. The background must not neces- 
sary satisfy any field equations. This simplified approach is somewhat physically 
incomplete but it allows to reduce the problem to solutions for quantum fluctua- 
tions only, that are supposed to obey the corresponding equations of motion with 
the background field as parameter. Generally speaking, the introduction of an 
external background field in a quantum field model is not trivial and not unique 
(see for example [7B]). 

A more consistent approach treats any field as consisting of a background 
part and fluctuations jSO], [22] • 

The total energy density of the system is the sum of the classical and quantum 
parts 

S = S^^ + S^, (3.2) 

where E'^^ is the local energy density as a function of coordinates of backgrounds 
$i defined as the 00-component of the canonical energy-momentum tensor 

i i 

and the quantum part is the vacuum expectation value of the canonical hamil- 
tonian of quantum field cf). It is assumed that a considered problem is static and 
the vacuum state in the Fock space is defined in such a way that the constructed 
object is local: 

S^{x) = £t =< 0\n{x)\0 > . (3.4) 

The object defined in such a way is generally speaking non-local, because 
it depends on the background potential in the whole space. The integral of it 
over the space provides however a well defined global object, which we are even 
interested in, see below. The S"^ has been shown in Sec. ()2.2|) to be proportional 
to the sum over the spectrum of Ti 

Here the spectrum is defined as usual 

n^nix) = a;2$„(x), (3.6) 
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and is supposed here to be discrete for technical reasons. 

The energy £ is treated as the density per unit volume. Since we are interested 
in the global energy E of the field system, it should be integrated over the volume 
V, 

E = J d^x S{x). (3.7) 

V 

3.1 The scalar field with the scalar background 

Consider first a simple model described by the lagrangian density 



V{x) = a^x)\ (3.8) 

where the field is a quantum fluctuation interacting with a scalar background 
field $(a;). The lagrangian can also contain the self-interaction term A$^. It will 
be shown later that such a term is needed for the interpretation of a counter term 
arising from the quantum vacuum fluctuations and leading to a renormalization 
of the self-interaction coupling constant A. 

The corresponding classical energy density of the background and the vacuum 
energy can be obtained in a standard way and they are 

= l{9o$9°$ + (V$)2 + M2$2} (3.9) 

= 1 < 0\dod° + + + V{x)\0 > . (3.10) 

The equation of motion for the fluctuations corresponding to the (j3.8p is the 
Klein-Gordon equation : 

(^0^°- A + m2 + 1/(x))0(x) = (3.11) 

From now on we restrict ourselves to cylindrically symmetric potentials V{x) = 
V{r) 

For the scalar field 0(x) in cylindrical coordinates {cp, r, z) it is usual to make 
the ansatz 

0(x) = 0(r)e-'P°"''e-''^e-*f^^ (3.12) 
which leads to an ordinary second order equation for the radial part: 

+ + (^^ - ^(0 - -i) I 0(r) = 0, (3.13) 
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with k"^ = — m^. This is an ordinary equation of the Bessel type, known to 
have an analytical solution representable in terms of the Bessel functions. 



The rectangular shaped scalar potential 



Suppose the construction of the background potential looks as follows: 

V{r) = ( ^° =^'''' < ' < (3.14) 
U otherwise 

Then fj3.13|) has the following solutions in the domains restricted by the cylinders 
r = Ri and r = R2: 



I. r<R, ^li{r) = Ji{kr) 

II. Ri<r<R2 (pi'^iir) = Ci' H^j^\qr) + Ci' H'f\qr) 

III. r>R, ct>i'J{r)=C[''H['\qT) + Ci''Hf\qT) (3.15) 

Here q = y/k'^ — Vq, Ci are some normalization coefficients. If the solution 
4>ii{r) is chosen to be regular in r — 0, then the coefficients C|^^ and Cf^^ are 
the corresponding Jost function fi{k) and its complex conjugated, fi{k), (which 
will be discussed later, see Sec J5.2|) . 



3.2 The Dirac spinor field with the magnetic background 

The next model to consider is the spinor electrodynamics: 

C = -^F^.F^-' + tP[tY{d, - leA,) - m]7A, (3.16) 

where the is the classical electromagnetic background, 

F^, = d^A, - d,A^, (3.17) 

as usual. We consider only a static problem. The background cannot depend on 
time and must be pure magnetic, the electric field is excluded, E = 0. 

Like the scalar case considered above we derive the corresponding classical 
energy of the background and the vacuum energy of the spinor quantum fluctu- 
ations. 

£'^1 = ^F^^F"^ + F^^F^ = ^B^, (3.18) 
where the B is the magnetic fleld, 

=< o|7%{7(-2V - eA) + m}|0 > (3.19) 



23 



Let us chose the potential to be of a rectangular shape (like the scalar one 
discussed above): 

r $ a(r) ^ 

A = -^e-;, (3.20) 

so that it contains only the A^p component and possesses a cylindrical symmetry 
again. The radial part profile function a(r) is taken to be 



0, 

1 { 



r < Ri,R2 
1, r>Ri,R2 



(3.21) 



Then the profile function for the magnetic field B(r) = ^h{r)ez is 



r or 



0, r<R,,R2 
^Ri<r< R2 
0, r>Ri,R2 



(3.22) 



with K = '^r,-2 ^ ■ The shape of the background can be interpreted geometrically 
as an infinitely long flux tube empty inside. 

We proceed now with the Dirac equation for spinor field tplx) of mass m and 
charge e 

«7^^ -m- e^AA ^/j = 



(3.23) 



where is given by (j3.2(Jj) . The gamma matrices in our representation are 
chosen as in [221 IS] '■ 



7 



7 



^3 
-^3 

—iai 
i(7i 



7 



7 



10-2 












/ 








(3.24) 



(Tj are the 2x2 Pauli matrices, / is the 2x2 unit matrix. 

Now we follow the standard procedure and separate the variables. Using the 
ansatz 

(/.(r, if, z) = e-^P«^°e-^?'^"^(r, ^) (3.25) 

0(r, if) 



we obtain the equation for the 4- component spinor \l/(r, cp) 

Po + L- mas Pzl 

pj po + L + mas 



^(r,v9) = 0, 



(3.26) 
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where L = i ai{di + ieAi), 0, x are the two- component spinors. 

i=l 

Respecting the translational invariance of the system it is convenient to solve 
the equations only for pz = 0. Then the system ()H.2(j|l decays in two decoupled 
equations, for (p and x respectively, which are distinguished only by sign of mass 
m. 

We consider the first of these two equations. By meaning of the standard 
ansatz 



(f){r,ip) 

I is the orbital quantum number), it reads: 



(3.27) 



dr 



Po — m 

l+l-Qajr) 
r 

Here the radial part $(r) denotes 



_a l-na{r) 

dr r 



Po + m 



(3.28) 




and 



2^' 



(3.29) 



(3.30) 



In order to use later the symmetry properties of the Jost function we redefine 
the parameter / in ()3.27|) (orbital number) as u according to 

/ + i /or / = 0,1,2,... 

(3.31) 



-/ - i for I 



-2,... 



For the given construction of potential ()3.2H) we get three equations and three 
types of solutions for 3 of space respectively: 



Domain I: r < Ri. The free wave equation (the magnetic flux is zero), the 
solutions are the Bessel functions J^^i{kr) 



(r) 



$+(r) 



1 



-VPo + mJ^_i{kr) 
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\/n^^Ju-^{kr) 



This solution of equation ()3.28p in the domain I is chosen to be the so called 
regular solution which is defined as to coincide for r — > with the free solution. 



Domain II: Ri < r < R2 The equation with a homogeneous magnetic field 
^h(r) = has the solution 



$7,(r) = c; 



c- 



(po+m) a+1 F Q _ hi a + 2- 

2(q+1) ' ^ l-f^U-L 2fl' " ^ ^' 



2/3' 



2 ^ 

r'^e-^iFi(-|,a + l;^) 



1 2a ( /3r2 - 



PQ—m r ^ 2 



V — — — 

e 4 



2/3 



- /3r^ - 
> 2 ' 



e 4 ^Fi(-f^-a,l-a;i|-^ 



(3.32) 



^}j{r) is the same, but a replaced by a 
here: a = z/-| + ^for/>0 

a = ^ - z/ - I for / < 

/5 = % 

iFi is a confiuent hypergeometric function P, [iH] . 

The coefficients Cj, are some constants which will be unimportant for ex- 
pressing the Jost function. The indices u and / denote "upper" and "lower" 
components of spinor respectively. The lower index "i" or "r" corresponds to 
"regular" or "irregular" part of solutions dependent on the behaviour of the 
function if continued to r = 0. If the external background vanishes {Q — > 0), 
contributions of irregular parts to the solution disappear. 



Domain III: r > R2 The free wave equation outside of the magnetic fiux has 
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the solutions 



iii\ 



i.i- 



2 VPo - 



riJiAr) 



^(1) 



(A;r) 



^(2) 



7] 



2 v^po - m 



riiiAr) 



+ \f^{k)—^ 
2 v^po - 



l^f^k) 



+ ^f;^ik) 



and fAk) is the corresponding Jost function in accordance with the definition 
SecEa (EIHl)- 

3.3 The Dirac spinor field with the abelian Higgs back- 
ground 

Let us consider the abehan Higgs model with the Yukawa type coupling. This 
model is known to be the simplest model with spontaneously broken symmetry, 
where the initially massless fermion field obtains a mass through the Higgs mech- 
anism during the U(l) symmetry breaking. It is interesting since it is also the 
simplest model possessing the solution of Nielsen-Olesen (NO) vortex type. This 
vortex will be treated in the considered context as a certain classical background. 
The complete lagrangian of the model reads 



C — C^j^ + + C 



(3.33) 
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where 



^A(|$r-^f, (3.34) 



(3.35) 



-f(/e|$|+e7M^)^. 



(3.36) 



The covariant derivative is defined by = d^ — ieA^, = A^{x) -the gauge 
field (magnetic field), $ is a complex scalar field (Higgs bosonic field ). 

The interaction lagrangian in (j3.33p is chosen in this form to simplify the 
equations of motion for the spinor quantum fiuctuations It will be discussed 
below in Sec. ()3.3.2|) . 

The classical energy of this scalar-electromagnetic background is in general 



d X 



2 \ 2 -I 

|Do$l' + lA$l' + A(<l>t$-^j +-[E2 + H'] 



(3.37) 



where the E, H are components of tensor F^j, corresponding to classical electric 
and magnetic field respectively. 

Since we consider this field system on a particular background, namely the 
NO vortex, it will be useful for further calculations to rewrite E'^'' right now in 
the terms of the Nielsen- Olesen solution 



<l> = ^/(r)e*"^ qA^ = nv{r), A^ = A = A, = 0, 
V2 



(3.38) 



that provides: 



1]^ 
~2 



n 

lyt '2i 



and finally we arrive at 



E 



d 



TT 



r dr 



(3.39) 



(3.40) 

The ground state energy of the spinor quantum fiuctuations is derived similar 
to the previous case, as the vacuum expectation of the hamiltonian 



S'^ =< 0|7°{7(-^V - eA) + ^/(O^lO >, 



(3.41) 
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and the equation of motion for the spinor field can be obtained also in a 

similar way for the radial component \l/(r). 



— Po + ^f{r) 



^(r) = 0. (3.42) 



dr r ' V2 

Comparing with the case ()3.28|) we see that 

^/(r) ^ /z(r) (3.43) 
plays the role of space depending mass for the fermionic field. 
3.3.1 Nielsen-Olesen vortex 

The Nielsen-Olesen vortex is a soliton stabilized by a topological charge. Here 
we describe briefly the derivation of the NO equations , following |2|. After 
that SU(2) invariance of the electroweak model is broken, we have a Lagrangian 
possessing the relic U{1) -symmetry, i.e. the invariance under U{1) gauge trans- 
formations 

^ e*?^(^)<l>(x), A^{x)^A^{x)+d^x{x). (3.44) 
We start with the U{1) -invariant Lagrangian 



2 \ 2 



(3.45) 



where is the covarinat derivative = — iqA^, and the tensor Y^^ is like 
an electromagnetic field strength Y^y = d^Au{x) — duA^{x). By applying the 
Euler-Lagrange formalism as usually we obtain the equations of motion for both 
fields $(x) and A^{x) 

D^D^^ix) + 2A (^$"f(x)$(x) - ^x) = 0, (3.46) 

and 

d^'Y^^ix) = -iq^^{x) A, ^x) (3.47) 

The Nielsen-Olesen ansatz [HZj is introduced in cylindrical coordinates (r, r, if, z) 
according to ()3.42|) . where f (r), /(r) have to satisfy the boundary conditions 

/(O) = t;(0) = 0, f{r),v{r)^l as r ^ oo (3.48) 
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Figure 1: The profile functions of the Nielsen- Olesen vortex as functions of the 
radius r for qr] = 1. 

Substitution of the ansatz ()3.42|) into the ( I3.46l3.47jl leads to a system of two 
second order equation with respect to functions v{r), /(r): 

/"(r) + ^ - n'^[l - v{r)Y + \v\l - f{rf)f{r) = 0, 

v"{r) - ^ + gV/(^)'[l - v{r)] = 0, (3.49) 
r 

which are in general not to be solved analytically, but it is possible to obtain 
several numerical evaluations for v{r) and f{r). 

If we are only interested in the shape of the solutions v{r), f{r), it is convenient 
to introduce a re-scaled variable r, r = p/ [qr]), so that dependence on parameters 
in equations (j3.49j) can be reduced entirely to one parameter /3 = Now, 
since the parameter /? is specified, the corresponding shapes can be reconstructed 
using the boundary conditions ()3.48|) and the asymptotic behaviour of solutions 
f (r), /(r) at r — i> 0. We put 

/(r) = ar + . . . v{r) = br'^ + . . . , (3.50) 

where a and b are some constants, which are fixed by the asymptotic condition 
(j!riH|l at r ^ oo. 

Two examples are shown in Fig l3.3.l1 for qi] = 1. The corresponding values 
of the constants a,b are (a, 6) = (0.26817,0.17481) for (3 = 0.3 and (a, 6) = 
(0.79958,0.42848) for (3 = 6. As the asymptotic values for r oo are reached 
with exponential speed it is in fact sufficient to consider r up to ^ 15. 

3.3.2 Interaction with the spinor field 

The interaction term C^j^^ in the form ()3.36|) leads to the equation of motion 
for spinor field like the ordinary Dirac equation with magnetic background, that 
allows to save a lot of analytical calculations by using the results of the solved 
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problem of magnetic string, discussed above. The interaction ()3.36|) possesses 
the U(l) gauge invariance and the couphng constant /e is dimensionless, what 
should serve to keep the renormalizability of the model. However, the latter is 
violated through the interaction of the spinor with the absolute value of the scalar 
background $, that is originally complex, thus the term 

/e|$| = ^/(^WTWy (3.51) 

has a non polynomial structure. 

As a consequence, the renormalization of the one-loop vacuum energy requires 
the introduction of an additional counter-term, as considered below, Sec J5.5l 

The correct renormalizable form of the interaction lagrangian was supposed, 
e.g. in [7^1 inn] for N=2 supersymmetric extension of this model. 

Variation of this lagrangian (|3.36p with respect to the spinor \[' leads to the 
equation of motion: 

{i^D^- femx)\}^ix)=0. (3.52) 

One can see that the scalar background ^{x) plays the role of the variable 
space-dependent fermionic mass. It suggests that this problem can be treated in 
much the same way as that one for the massive spinor field in magnetic back- 
ground. 
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4 The zeta functional regularization of the ground 
state energy 

The ground state energy of the quantum field derived in Sec. ()2.2| ^ as the spa- 
tial integral over the infinite volume of the sum over the spectrum of hamiltonian 
is not well defined mathematically, since the spectral sum is not necessary con- 
vergent and the spatial integral as well. The one loop vacuum energy is known 
to contain ultraviolet divergences. Even in all further applications the energy 
defined as (HTMfUP^ is divergent. 

Quantum field theory operates with a variety of regularization procedures to 
treat infinite objects reasonably. We use the so called zeta-functional regulariza- 
tion. 

This procedure is based on the representation of divergent sums and series 
through the zeta function and is described in the present section below. 



4.1 The Riemann zeta function 

The Riemann zeta function (or simlpe "zeta function ") is the one of the central 
objects of the number theory, defined as 



oo ^ 

C(^) = E;;i- (4-1) 



n 

n=l 



It was introduced by Riemann in 1859. It is known that for real s the series ()4.ip 
converges for s > 1. 

The next important step to proceed with the zeta function is to gain an 
integral representation of the infinite sum (j4.ip . We perform the standard repre- 
sentation of the Euler Gamma-function in the way 



r(s) = j dt f-'e-' = n' J dt f-'e"""'. (4.2) 





Then in the product 



„ oo 
oo „ 

as)T{s) = J2 dtf'\--\ (4.3) 

n=l{ 

the summation over n and integration over t can be interchanged under assump- 
tion that s > 1, that admits the resummation of e~"* as a simple geometrical 
series, and finally we obtain the desired integral representation 



oo 

s-1 
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We can proceed from this expression with the analytical continuation of ({s) 
in the complex plane. It allows i.e. to establish a number of important properties 
of ({s) as a function of complex variable s. Here are referred to some useful ones. 

• The zeta function is an analytical one-valued everywhere in the complex 
plane excepted the point s = 1 where it has a simple pole with residue 1. 

The corresponding Laurent expansion around the s = 1 reads : 

(is) = + 70 + 7i(5 - 1) + l2{s - If + ... , (4.5) 

s — 1 

where the irrational coefficients jk are to calculate as 



lim 

7fe =n -> oo 



(4.6) 



and the 70 = 7 is the Euler-Mascheroni constant; 

• For the s = we have 

C(0) = -\ (4.7) 

• The real zeros of the zeta function are located at negative even integers 

C(-2n) = 0; (4.8) 

• Furthermore the zeta function has infinitely many complex zeros in the 
strip < 3fJs < 1 but no one outside of it. 

The one of the most important generalizations of the Riemann zeta function 
is the spectral zeta function of some operator V . This kind of spectral function 
(the function defined on the spectrum of operator) is known in mathematics since 
1949 

It is defined as 

Cv{s) = Tr (4.9) 

In particular if the spectrum A„ of V is discrete, the ()4.19p means a sum over 
the spectrum A^^. 

Cr{s) = TrV-' = J^i^n^ (4.10) 

n 

For example, the problem of a linear oscillator or a charged particle in homo- 
geneous magnetic field in quantum mechanics [56- leads to the spectral problem 
of the operator 

1 d'^ d 
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on the space of bounded functions H{z). The solutions of this problem are 
the Hermite polynomials Hn{z), and the corresponding spectrum A„ is n = 
0,1,2,3.... If we redefine Vz as follows: 

V+ = V, + 1 (4.12) 

to exclude n = from the spectrum, the corresponding zeta function Cy+{s) will 
be even the Riemann zeta function 1)4.11) . Obviously there are more different 
spectral problems possessing the same spectrum and the same spectral functions, 
as e.g. zeta function respectively (so called isospectral problems). 

Analogous to the Riemann zeta function , the generalized zeta function (con- 
sidered below, I4.19|) can be analytically extended on the complex plane s as a 
meromorphic function possessing simple poles [JAi- In particular (-pis) is regular 
at s = if P is of Laplace type. 

One can define also a determinant of the operator formally as 

detP = J](A„). (4.13) 

n 

It can be represented through the zeta function of the operator V in the 
following way 

exp l-^^CAs)] = n(^")^'"^"^- (4-14) 

For s = it gives 

(4.15) 

s=0 

4.2 The generalized local zeta function 

One can also extend the definition of zeta function corresponding to the spec- 
trum of operator P and define a bilocal function depending on coordinates of two 
points (x,y). 

If we have a typical spectral problem for eigenvalues A(„) and eigenfunctions 
</)(„) (x) of the positive defined elliptic differential operator P: 

P0(„)(x) = A(„)0(„)(x), (4.16) 

where the eigenfunctions (x) are the functions of coordinates (x) on the man- 
ifold M with boundary and obey the boundary conditions 

i30(„)(x)|a^=O (4.17) 



detP = J](AO = exp |-^C7'(s)| 
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( B is some boundary operator), then we can define 



^ , I . ^ 0(n)(x)0^„)(y) 

Cp(x, y|s) = 2^ . (4.18) 

(n) ^("O 

Thereby the original zeta function is naturally constructed as a derived global 
object being the trace of the local two-point function ()4.18|) 

Cp(s)=TrCp(x,i/|s), (4.19) 

that coincides with the definition ()4.10|) (the trace means summation over all 
discrete indices and integration over continual variables {x,x)). It is assumed, 
that for the operator P such a trace exists. For simple examples, where the 
spectrum of P is known explicitely, we suppose the sum in (j4.1(jp to be convergent 
for 3fJs being large enough. More precisely, this assumption is based on the Weyl's 
theorem, which presumes, that for a second-order elliptic differential operator 
on the D-dimensional flat euclidean manifold Ai the eingenvalues A„ behave 
asymptotically for n — > cx) as 



^ r(D/2) 



^- (4-20) 

volM 

(We work basically with this class of operators only). In particular, this formula 
provides the convergence of the sum ()4.in|l for 3?s > D/2. [55] . 

Here it is also suitable to mention, that the local zeta function possesses the 
following useful feature 

PCp(x,y|s)=Cp(x,y|s-l), (4.21) 

that can be checked immediately and is useful when applied to problems con- 
cerned with the renormalization of Green function . For s = one obtains from 
(j4.2ip the initial conditions 

Cp(x,y|s = 0)=(5(x-y), (4.22) 

where it can be seen explicit, that the local zeta function for s = coincides with 
the corresponding Green function . 
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4.3 The heat kernel 

The heat kernel provides a number of tools to use within the zeta function reg- 
ularization method. Mathematically this object is closely related to the general- 
ized zeta function which is defined on the spectrum of some elliptic self-adjoined 
operator. 

One defines the simplest (global) heat kernel as a sum over the spectrum A(„) 
of some operator P: 

irp(t) = J]e-*^("). (4.23) 

in) 

We can also define the local heat kernel K{x,y\t): 

i^p(x,y|t) = 5^0(„)(x)0^„)(y)e~*V), (4.24) 

(n) 

where are the normalized eigenfunctions of the operator P. 

The trace of Kp{'x, y\t) (here trace means summation over all discrete indices 
and integration over continuous variables x), 

Kp{t) = Tr Kp{^,y\t) (4.25) 

obviously gives the corresponding global object ()4.23|) . The kernel ()4.24|) satisfies 
the equation 

|^ + p|A:p(x,y|t)=0 (4.26) 
with the initial condition at t = 0: 



ir(x,y|t = 0) = 5(x-y). (4.27) 

A natural example of the operator P is the minus Laplacian —A. For this 
case the fl4.26|) becomes the usual heat equation 

^-A|ir_A(x,y|t) = 0, (4.28) 

where the name "heat kernel " comes from. In the applications of mathematical 
physics P means usually the hamiltonian of field configuration. Some suitable 
examples are considered in further subsections. 

Transition from the heat kernel Kp{'x,'y\t) to the local two-point dependent 
zeta function 

Cp(x.y|.)^^ ^-''y'"''^' (4.29) 

(n) (") 
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can be made by using the standard representation of F-function 

oo 

r(s) = 1- j dt t'-^e-'^. (4.30) 



This equation together with (j4.29j) yields: 

Cp(x,y|.) = I dt^i^^(x,y|t). (4.31) 



Taking the trace of both sides of this identity we obtain the connection be- 
tween global objects. 

oo 

Cp(^) = / dt^^^Kpit). (4.32) 



Another important spectral function also representable through the heat ker- 
nel is the Green function G{x,y) of the operator P, that is defined to obey the 
equation 

PG'(x,y) = 5(x-y) (4.33) 

(it will be considered in Sec J5.2l in context of the Jost function ). 
The function G{x,y) is usually constructed as 

G(..,).^WkM, ,4.34) 

(n) 

(that obviously satisfies the ()4.33|) ). Then representing 1/A(„) as the integral 

1 



oo 

tx 



A{n) 



rfte~*^w, (4.35) 







we get 



oo 

G(x,y) = j dtK{^,y\t). (4.36) 



In the historical overview it was the first physical application of the local heat 
kernel , where the Green function was proposed to be represented as integrals 
like ()4.3fjj) of a kernel satisfying the heat equation over an auxiliary coordinate 
("proper time"), introduced by Fock, 1937. 
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In appllications to the ground state energy it is important that the heat kernel 
has the power-law asymptotics for small t. For all strongly ellilptic boundary 
value there is an asymptotic expansion 



K{x,y\t) 



1 



(4.37) 



n>0 



and the asymptotics of the global kernel K{t) (or heat trace) reads respectively 



as t — s> 0. The d denotes the space dimension and the coefficients called 
the "heat kernel coefficients ". It will be shown below that the zeta-regularized 
ground state energy for massive fields in the limit of large mass m — >■ oo ()4.115| 
SecE} can be expressed in terms of these coefficients. 

4.4 Heat kernel coefficients 

There are various methods to calculate the heat kernel coefficients a„ One of 
this mostly applicable in physics is the consequent iteration of the equation for 
the local heat kernel K{x,y\t) (j4.26j) with an initial condition (j4.27p using the 
ansatz with the free solution KQ{x,y\t). This method goes back to DeWitt [21] 
and is called in some sources as the "DeWitt iterative procedure" . 

This procedure was originally developed for manifolds without boundaries [21] 
with non zero curvature and can be generalized also for manifolds with boundaries 



The short description of it, following [TH] is concluded to the recursive solution 
of the chain of equations 




n>0 



(4.38) 



[nmnn]. 



(j + {D^'a)D^)a2j + Pa^j.i 
that are derived by substituting the DeWitt ansatz 







(4.39) 




(4.40) 



n>0 



into the equation ()4.26|) with the initial condition 



ao = 1 



(4.41) 



and K^^\x,y\t) obeys the equation ()4.26|) as t = 0. 
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4.4.1 Scalar quantum field in scalar background 

Let (f){x) be a quantum scalar field of mass m, interacting with the classical 
background potential V{x). The operator P is a hamiltonian Ti. of this system 

n = -A + m^ + V{x). (4.42) 

The variable x includes now all spatial components, A = d^di. The equation 
flOH|l reads: 

1^ - A + m2 + \/(x)|ir(x,y|t) =0. (4.43) 
The ansatz for the local heat kernel expansion must be chosen in the form 

m 

K{x,y\t) = K^^\x,y\t)Y,an{x,y)e; (4.44) 

n>0 

where an{x,y) are the local heat kernel coefficients which depend on two points 
x and y. Here K^^\x, y\t) is the zeroth order term satisying the initial condition 
at t = and 

for the operator ()4.42|) the K^^\x,y\t) is chosen to be 

K^'\x,y\t) = J^2^^P{-^^^ - ^-'} • (4-45) 

Considering the limit of t — > in 14.451 we assure that this ansatz 
The kernel fl4.45j) obeys also the initial condition ()4.27|1 . Now substituting the 
ansatz E. 441 into ()4.43|) one obtains 

^ - d^d^ + V{x) + m") 

(4.46) 

We equate the coefficient in front of each power of t to zero. This gives the 
following chain of equations 

(x - yyViao{x,y) = 

Aao{x, y) + Vao{x, y) + ai(x, y) + (x - ?/)*Viai(x, y) = 
Aai(a;, y) + Vai{x, y) + 2a2{x, y) + {x - y)*Via2(x, y) = 

(4.47) 

with the initial condition ao{x,x) = 1; Vj = 



(47^t)«!/2 



exp 



yy 



At 



0. 
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Now let us apply the differentiation on the first equation of ()4.47p we 
obtain, that 

29'=9fcao(x, y) + {x- yfdAa^ix, y) = 0. (4.48) 

If the coordinates x and y coincide, this equation gives Aao(a;, x) = 
A similar procedure applied to the second equation of (j4.47|) results in 

A2ao(x, y) + 2Vao(x, y)VV{x) + ao(x, y)AV{x) + Aao(x, y)V{x) 

+3Aai(x, y) + ((x - y)V)Aai(x, y) = 0, (4.49) 

where Aao(a;, y) is a solution of ()4.48|1 and A^ao(a;, y) satisfies a similar equation 

A2ao(x, y) = -i((x - y) V) A2ao(x, y). (4.50) 
In the limit x = y { 14.491 ) simplifies 

Aa,{x,y)\,=y = -^AVix) (4.51) 

and 

ai{x,y)\,=y = -V{x) (4.52) 

(the laplacian A acts on the x-variable) 

Using these results one obtains from the third equation ()4.47|1 for a2{x,x) 

a2{x) = a2{x,x) = {Aai{x,y)\^=y + V{x)ai{x,y)\^=y) . (4.53) 
Finally we get 

a2{x) = a2{x, x) = ^AV{x) + ^V{xf. (4.54) 
Further coefficients 03, 04, . . . can be calculated in a similar way jlHl El E]- 
4.4.2 Scalar field in constant vector background 

Following the procedure described in the previous section, we construct a chain 
of recursive equations for ai{x,y). 

It starts as previously with the same ansatz ()4.44|1 into the equation ()4.26|1 . 
The operator P reads now 

P = -A + m^ + V{x), (4.55) 

where the Laplace operator A is modified now through the gauge field y4^(x), since 
it is constructed of gauge covariant derivatives A = g'^'^Df^D^, = dfj, — ieAfj_{x) 
One obtains a similar chain of equations for ai{x,y) : 
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(x - yfD^ao{x,y) = 
{-D^ + V{x))ao{x, y) + {l + {x- yYD^)a^{x, y) = Q 
i-D^ + \/(x))ai(x, y) + {2 + {x- yYD^)a2{x, y) = 



and arrives finally at the following relations: 



aQ{x,y) = exp 



y 

-ie J A^{z)dz^ 



(4.56) 
(4.57) 
(4.58) 
(4.59) 



(4.60) 



what yields for the coincidence limit x = y again ao{x,x) = 1, Fp^, is as usually 
the electromagnetic field strength tensor dpA^, — d^Ap] 



ai{x,x) = {D^ - V{x))ao{x,y)\ 



x=y 



-V{x) 



(4.61) 



D'ai{x,y) = -{D'D' - D'V{x))ao{x,y)l,=y 



-^e'Fix)' - ^AVixy, (4.62) 



a^ix^x) = ^{D' - Vix))a,{x,y)\,=y = -^Fixf + ^Vix)' - ^A\/(x). (4.63) 

4.4.3 Spinor field interacting with gauge vector field 

If we are interesting in the heat kernel expansion for the Dirac field interacting 
with external gauge vector field. 
We start from the Lagrangian 



£ = --Fp^F^" + ij[ir{dp - leAp) - m]^. 

The equation of motion for the spinor field is the Dirac equation 

[i-l''{dp - ieAp) - m]ij{x) = 0. 
The corresponding Hamiltonian is: 



(4.64) 



(4.65) 



(4.66) 



(Here the hamiltonian Ti is not to be confused with the definition ()2.35[ Sec l2.2|l . 
and is to understand by means of primary quantization, namely it provides the 
eigenvalues po of energy, acting on the static states ipi^^, = 6"*^°^ ip{^)- K can 
be obtained immediately from the Dirac equation ). 
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The energy spectrum A„ for static problems will be defined from the equation: 

HMx) = KMx), (4.67) 

the eigenf unctions ipni^) of the operator H depend on spatial coordinates x only. 

The global heat kernel K{t) cannot be constructed properly for this operator. 
Namely since the operator 7i is a fisrt order operator, it possesses a symmetric 
spectrum, that broads from — oo to oo and can get some discrete values between 
— m and m corresponding to bounded states. Thus K{t) in the form ()4.23|) is no 
more "well defined". 

To avoid this difficultness, one can use the operator Ti^ instead of Ti. In other 
words we can deal here with the sum over squared energies. Indeed, one can 
apply the operator H to the relation 14.671 with the resulting equation 

n^M^) = xlM^) (4.68) 

with 

n^ = -A + ^cr'^Fij + m^ (4.69) 
and by using the recursive procedure to the one described above one obtains 



ai{x) = ai{x,y)\x=y = (4.70) 
a2{x) = a2{x,y)U=y = -^F^, + ^P' + ^AF. (4.71) 

Here the notation F stays for a^'^Fij, where cr*-' = |[7*7-'] as usual, and the 
spinor indices are omitted. 

In applications of global ground state energy the trace of heat kernel coeffi- 
cients (14.7111 is used: 



02 = Tr / d^x a2{x) = ^F^.. (4.72) 
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4.4.4 Spinor field interacting with an electroweak background 

The electroweak model with spontaneous broken symmetry is represented by 
the Lagrangian SecO Eq. (j3.33|) 

The heat kernel coefficients ai and 02 can be calculated by using the formulae 
for the squared Dirac operator, ()4.7H) established above. 
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{ilp+ /i(r)) {ip- /x(r)) = -D^ + - ^/eT'^^ I $ I -/e I ^ ? 

= -D'^ + V (4.73) 

with cr^*^ = f[7'^,7'^], /i(r) corresponds to the definition of radial-dependent 
fermionic mass in Nielsen-Olesen notations Sec. I^ ()3.43|) . The general expressions 
for relevant heat kernel coefficients are 



ai = Tr y d^x {-V), 



«2 = Tr y d'x {-^F'^u + \V' - ^AFj , (4.74) 

where the trace is taken over the spinor indices. Substitution of the expression 
for V defined by ()4.73|) yields 



7C dr r (/i(r)^ — mg) 



a, = 8n Idrr ( + ^ (/i'(r)^ + (M^)^ - ^D') 



(4.75) 



4.4.5 Other approaches to the calculation of the heat kernel 
coefficients: the expansion of the exponent 

This approach makes the structure of the global heat kernel coefficients quite 
transparent. Its advantage compared to the DeWitt procedure is that it works 
with the global objects a„ immediately without need in intermediate relations for 
two-point dependent coefficients an{x,y). 

The approach is based on the generalized representation of Kp{t) ()4.25p in 
the symbolic form: 

Kp{t) = Tr L2e~'^. (4.76) 
It is convenient to replace this object by a more general one: 

Kp{f;t)=Tr L2{f e-'^}, (4.77) 

where a smooth function / plays the role of smearing function and keeps Kp{f; t) 
finite even in an infinite space. 
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Now consider the calculation of the generalized heat trace ()4.77|) . To this end 
we have to define the basis in space where the trace is defined. Following the 
[TH] we use the simplest basis of the fiat waves in cartesian coordinates. 

I >= e'*^" (4.78) 
The operator P is some second order operator of the Laplace type 

P = -{g^"'D,D, + E), (4.79) 

containing some (possible matrix valued) function E. The metric g^'^ is fiat. D 
is the covariant derivative 

D = d + uj, (4.80) 

that includes the connection u (general notation for geometrical or gauge con- 
nection). Carrying out the calculation for heat trace ()4.77|) of this operator in 
the fiat space we obtain: 



Tr f e-*^ > 



[ (Tx [ e"^'^^ tr f(x) e"*^ e^^^ 

J J (27r)« ' 



f J^tv fix) ^ti(D,+ik,)(D^+^kn+E)_ (4 8^) 

J (27r)" ''^ ' ^ ' 

The notation tr denotes now the trace over all remaining indices (spinor, 
colour etc.) 

The exponent entering here can be represented as a product e"*'^ gt{«fc^i5M+«fc^-D''+-DM'D^+-B) 
and then the second exponent can be expanded in the limit t — > with the result 



^ti^k^D,+^k,D^+D,D^+E) ^ ^ ^ ^^2 ^ ^) _ 2t\kD) + 

4+3 

t\D^D^ + D^E + ED^ + E^) - —UkDfE + E{kDf + {kD)E{kD)] - 

6 

^[{kDfD' + D\kDf + {kD)D\kD)] + ^-^{kDf 
b 24 

Keeping in mind that the integrals of terms with odd powers of k'^ disappear, 
these terms are already omitted. 
Using the integrals of type: 

! -^e~"''k^'k^'\.k^^^ = , \ ^^rMi,M2,.../.2,i 
} (27r)" (47rt)"/2 2j(t)J 

J = 0,1,2,..., (4.82) 
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where [. . . ] denotes the full symmetric combination of indices /li, one can carry 
out the integration over k and get 

-{D'^D^ + D'^E + ED^ + E^) {D^E + ED^ + D^^ED^ + 2D'^D^ + D^D^D^) + 

2 3 

-{D^'D'^D^.D^ + D^D^ + D^'D^D^) + 0{t^)]. (4.83) 

Finally, introducing the corresponding notations for commutators of connec- 
tions appearing thourgh commutations of covariant derivatives 

^f,u = d^ujy - dyUJ^ + uj^ujy - ujyuj^,, (4.84) 
we rewrite the result in the form 

^P^^) = J^^2 I ^"^tr f{x)[l + tE++t\^E' + ^E,,, + ^n,Mn + 0{t% 

This result is very general comparable with the explicit formulae for certain 
cases obtained above by the DeWitt recursive procedure. In particular if the 
Maxwell tensor F^u corresponding to the electromagnetic potential is treated 
as an endomorphism E while the potential itself plays at the same time the role 
of the U(l) gauge connection, we identify the contribution of order t"^ in ()4.85p 
with the contribution of 02 (I4.7ip . The F = F^^a^^ corresponds to the E and 
F^u = d^Ay — dyA^ to the Vt^y respectively, the covariant derivative ;^ in flat 
space is the ordinary one d^. 

4.4.6 Some general properties of heat kernel coefficients 

We conclude this review of the heat kernel coefficients with some general rules 
inherent to construction of this objects. 

We restrict ourself again on certain class of problems related to the physical 
models being handled in this paper. We consider in general the field configura- 
tions in the fiat Minkowski space-time without boundaries. However a boundary 
can be introduced in intermediate steps with suitable boundary conditions, that 
will be removed later by the transition from the finite volume inside the boundary 
to the infinite one. Then in each case it must be respected carefully that the heat 
kernel expansion is transformed properly by this transition. 

The general construction of the heat kernel coefficients a^, = 0, 1, 2, . . . for 
some operator P of the Laplace type (j4.79|) with the smooth connection VL and 
smooth endomorphism E acting on the manifold M without boundaries reads 
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1. ak{f,P) = for k = 2j + l 

2. afc(/) P) for k = 2j are locally computable in terms of geometric invari- 
ants 

afc(/,P)=tr j (rx^{f{x)ak{x,x\P)} = Y,^v j d^xf{x)u'Ai{P), (4-85) 

M ' V 

where Al are all possible local invariants of the dimension k constructed from 
E,Q, and their derivatives. The constants appearing here contain some ratio- 
nal constants of kind p/q, p,q E N and powers of n. These constants define the 
algebraic structure of heat kernel coefficients . In fact since the general structure 
of the coefficient is known (|4.85p . the problem is reduced to the calculation of 

only. It can be reached by using relations between . To establish the cer- 
tain relations one introduces some additional parameter and applies the known 
variational properties, described more detailed in [7HJ ES] . 

The important remark about the smoothness of E must be mentioned. If 
the function E{x), (for example the scalar potential V{r), SecOJ is not infinite 
smooth, E{x) ^ C°°, then the statement 1. is no more valid. So for the magnetic 
background with a step function the coefficient 05/2 appears to be non zero. It 
will be calculated in Sec. 15.61 
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4.5 Zeta functional regularization 

The application of methods using zeta function and heat kernel to the regular- 
ization problems of quantum field theory was introduced independent by number 
of authors 1211 (HH] ■ The generalized zeta function was used for example in 
explicit form in problem of renormalization of cosmological constant The 
usage of heat kernel expansion for regularization of effective action comes back 

to ISB. 

The regularization procedure using the zeta function can be considered to 
some degree as a certain extension of the wide used method of the dimensional 
regularization in the perturbation theory. 

The main advantage of this approach is that it allows to identify and separate 
divergences as a pole terms in power expansion for small regularization argument 
s 0. 

The zeta functional regularization is allied closely with so called dimensional 
regularizations while it uses an additional space dimensionality that will sent to 
zero and thus detect the divergent terms by arising poles at s ^ 0. 

The basic idea of the zeta functional regularization method is the represen- 
tation of the regularized object as the generalized zeta function of an (elliptic) 
operator P. 

The essence of the dimensional regularization consists in the introduction 
of some additional (in fact subtractional) dimension d in order to reduce the 
dimension of integral measure to make the integrals "conditional" convergent. 
After that the performance back to the limit d —>■ allows to get the separated 
divergent contributions. They appears as a terms containing pole, that will be 
cut out (renormalized out). 

The zeta functional regularization procedure is also based on the introduction 
of the additional dummy dimension s that will be abrogated afterwards and the 
pole-containing terms will be regularized out, that brings the final results, coin- 
ciding with the one given by dimensional regularization if the extra dimensions 
are flat [S2]- 



4.5.1 The space dimension and structure of poles 

In the review of this topic we follow ^01 EH] and start with the definition of 
zeta function in the form SecEJ ()4.19p 

CAs) = Tr '{if.-'Vr} = (^) . (4.86) 

(n) ^ ^ ^ 

This definition is generalized compared to the ones ()4.19I4.10|) by the introduc- 
tion of the dimension-adjusting constant fi to keep the zeta function dimensionless 
for all s. The spectrum of V can possess some values A(„) equal to zero, which 
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make the Cv{s) ill-defined for positive s. To correct the situation we suppose all 
zero-modes to be excluded in ()4.86p . that is noticed by the prime. 

Let the operator V again be a second order elliptic differential operator Dd 
of the Laplace type of dimension d (Sec J4.Hj) . (Further we discuss here generally 
this class of operators only). Second, we will work in the Minkowski space-time 
d = 4 and assume that only the second order time- derivative is present. Thus it 
can be separated explicit 

D, = + D,. (4.87) 

The integration over x° in the global zeta function C4{s) can be carried out as 
it was made in Sec J2.2l 

C4(.) = ^^^^C3(.-b, (4.88) 
V47r r s 2 



(here time is intergated oud and T means therefore the age of the system). 

Now we use the relation to the global heat kernel ()4.32|) together with the 
expansion of the heat trace obtained from ()4.H7j) 



(4.89) 

Since the diagonal part of the heat kernel contains exponentially suppressed 
terms of kind e~'^^'^^/*,jTU] it provides a non-zero contribution on manifolds with 
boundary as a boundary term present in 14.891 The (j4.32|) together with the 
recently mentioned ()4.89|) yields: 

where the sum runs over half-integers, 

Cn = ( / aMd^x + / hMd''-'y\ (4.91) 

VJn Jan ) 

and f{s) is an entire analytic function of s. 

Thus it can be seen that the location of the poles of zeta function depends 

on the dimension d of space-time. It is accentuated by the index d at the zeta 

function . In particular for d = 4 the zeta function (^{s) has no pole at s = 0, 

that allows to define the functional determinant in the form ()4.15p . 



4.5.2 The effective action 

This application of zeta function to the regularization of the one loop effective 
action goes back to the paper of Hawking jH21, where it was developed for the 
generalized case of curved space-time. 
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The one loop Feynman graph of perturbation theory is represented as the 
determinant of an operator as follows: It starts from the main object of the field 
theory in functional formulation, namely the generating functional for Feynman 
diagrams 



VgV<Pe'''^^^^\ (4.92) 

that is constructed as a functional integral over the product of spaces of all 
possible (admissible) metrics g{x^) and matter fields (j){x^) with measures Vg 
and Vcj) respectively. 

The 5'[(y', 0] is the classical action, that is also a functional of function sets 
g{x^) and 0(x'^). 

We restrict the discussion on the case of fiat space and omit the metric g in 
further calculations. Now assume that the quantum fields are the perturbations 
of classical ones (backgrounds) and represent the field variables in the form 

= 00 + 0. (4.93) 

In other words we decompose the arguments of the functional S" as a sum of 
classical, 'main' fields 0o and the quantum additions (fluctuations) 0. The latter 
are considered to be definitive small compared to the main classical ones. 

The backgrounds 0o are still supposed to satisfy the classical equations of 
motion and boundary or periodic conditions and provide the leading contribution 
to the action S. 

Respecting the smallness of fluctuations one can represent the complete action 
5'[0] as an expansion at 0o with respect to the fluctuations as 

S[g,cP] = S[g^,cP,] + sSl (4.94) 

Where the first corrections to the main action S'[0] are quadratic in the fluc- 
tuations; it follows further from ()2.2j) .that 

In Z = i5[0o] + \n j V4)e'^^^^\ (4.95) 

The general form of the corrections to the action generated by the fluctuations 
is supposed to be 

SM = -I J d^x 4>*V4>, (4.96) 

where the operator V, is the kernel of free action for the corresponding field; it is 
of second order for bosonic fields (e.g. □ + m^- Klein-Gordon kernel for spin 0), 
and of the first order for fermionic fields as well (e.g. Dirac operator i'y'^dfj, — m). 
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Suppose, the operator V to be elliptic and self-adjoint with the spectrum of 
eigenfunctions 0„ and eigenvalues A„ 

V<Pn = A„0„. (4.97) 

The eigenfunctions 0„ as a functions of coordinates will be orthonormalized 
as usual 

j d^X(l)l(l)m = 5nm. (4.98) 

Then any function can be decomposed in the basis 0„ 

4>=J2(^n(pn (4.99) 

n 

and the measure of the field space can be expressed through these coefficients a„ 

V4) = Y[fidan, (4.100) 

n 

fit out with the normalization parameter fi to keep the correct dimension. Then 
one arrives at 



that confirms the similar fact that the contribution quadratic in fields reduces to 
the evaluation of determinant. The logarithm of Z, treated e.g. in thermody- 
namical approach as the free energy reads now in terms of zeta function and its 
derivative 

lnZ[0] = iCp(O) + ^In Cp(0). (4.102) 

4.5.3 The ground state energy 

The simplest classical example of the zeta regularization is the zeta-regularized 
ground state energy considered in Sec. (j2.2p 

To consider this application of zeta function regularization procedure in details 
we follow again the [TUj . 

Here it would be suitable to mention three objects that are in principle dis- 
tinguished 



(4.101) 
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1. The vacuum energy, already defined in Sec J2.2l as the vacuum expectation 
value of the T^^ (and therefore called sometimes also vacuum energy): 

E,ac= [ <0|T0°|0> (4.103) 



2. The ground state energy (or Casimir energy), defined in ( SecOl Sec J2.2j) 
as the half sum over spectral values A(„) of hamiltonian 

^o = ^$^A(„) (4.104) 

in) 

3. The one-loop effective energy, being defined in view of the Sec. fl2.ip 

Eeff = ^Seff = ^ Indet D, = -^C(O) (4.105) 

with T = J dx^- age of the system, the index 4 indicates the 4-dimensional zeta 
function and the prime means that the possible zero-modes are not included. 
These objects under circumstances of problem are not equal ^U]. Following the 
notations of we use for further discussion in this section the dimensionless 
zeta function ^ of operator D which is related to the zeta function ()4.101 14.19P 
defined in Sec.4.1, by 

as) = tr '{(fi-'D)-^} = J^if^-'K)-^ = /i^^5^'(A„)- = fx'Xnis), (4.106) 

whereat D is supposed to be an usual second order operator. It has therefore the 
dimension [length]'"^ and the adjusting constant must be of dimension [length]'^ 
or [mass]. 

The 3-dimensional zeta regularized ground state energy is defined as 

£ris) = f 5^(A„/.-^)^/^- = l^sis - I) (4.107) 

n 

with the dimension-adjusting constant to keep the zeta function dimensionless 
for all s. Then it follows from the expansion ()4.9()j) with d = 3, that the £q'^^{s) 
is a meromorphic function of the complex argument s which has a pole in the 
regularization limit at s = 0. The residuum at this point is 

Res ..„{£„"»} = -\c^(<M)IH7r)-' = 



Because of the pole at s = the regulator s cannot be removed in a unique 
way. It means that the vacuum energy contains an ambiguity proportional to 
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the geometric coefficient C2. The minimal subtraction scheme adopted in [TU] 
supposes that 



So = ^PPUs-^), (4.109) 
where PP denotes the "principal part" and means 

PP Us -l)= lim ^ |C3(S - ^) + Csi-s - ^)} . (4.110) 

Another subtraction scheme will be considered below in Sec 15. 41 
Now since we establish the relation between heat kernels in 4 and 3 dimensions 
respectively 

K4{x,x,t) = -^Ks{x,x,t), (4.111) 
Ant 

then under assumption of ()4.88p we can apply the expansion ()4.90|) to both zeta 
functions with d = 3 and d = 4 and finally we obtain the relation between the 
ground state energy and the one-loop effective one. 

Eeff = Eo + ^mi) - ^(1/2)]^. (4.112) 

Here ipi-s) = '^^^fj'^^ is the polygamma-function arising from the expansion of 
F-f unction . 

Thus we can conclude, that the difference between Egff and Eq as well as the 
ambiguity of the Sq is caused by the non-zero geometric factor C2. If the latter 
disappears, the zeta function regularized ground state energy becomes unique and 
finite, and Eo = Eeff. However this is not the case for all problems considered 
below, since the heat kernel coefficient 02 is always non-equal zero. 

The vacuum energy defined by ()4.103p in the case that all interactions are 
switched off so that — > T°°gg then 

Eo = j d^x <0\T^^J0>, (4.113) 

as it has been assumed in Sec J2.2l 

Thus the ground state energy has been defined to be the sum over the spec- 
trum of eigenvalues Un of the hamilton operator TC (speaking strictly, eigenvalues 
of H are the squares of eigenfrequences u'^ ) 

n<^n{x) = ^^<f„(x) (4.114) 
and the regularized ground state energy can be just represented as 



Ens) = ^Cnis - i) = ^ = ^ E(^W + ^^y^-'- (4-115) 

(n) in) 



This expression will be used as a starting point in further calculations of 
vacuum energy for different models discussed below. 
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4.5.4 The large mass asymptotics for the zeta-regularized 
ground state energy of massive fields 

The regularized vacuum energy has been represented above Sec. (l2.2|l as the 
generalized zeta function of the corresponding hamiltonian Ti. Here it will be 
shown how the ultraviolet divergent contributions in the ground state energy can 
be identified. 

We start here with the expression for the zeta-regularized ground state energy 
of quantum field, derived in Sec. ()2.2|l . 

oo 

Sris) = -^Us -\) = -^j ^^^^W (4-116) 



and use the heat kernel expansion in the form 

^^('^^Ti^rtW^I]^"^"' (4-117) 

^ ' n>0 

that corresponds to the global heat kernel for the operator 7i 

H= -/\ + m^ + V{x), (4.118) 

d is the space dimension, n = 0, |, 1, |, 2, ... 

Then substituting this into £l^^{s) and re-expanding at s — one obtains 
the key formula for the renormalization of ground state energy . After collecting 
of equal powers of m together it acquires the form 



^ = -64^U^^"2j"°-24^"^/^ + 
/ 1 4/i^ \ m 



327r2 J IGvr^/^ 

--^(^- + \n^-2)a, + .... (4.119) 

The parameter t has the dimension thus the expansion in powers t at 

t — corresponds to the 1/m expansion. The first several terms of this expansion 
quoted here contain all ultraviolet divergences as poles at s ^ 0. 
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5 Analytical treatment of the 
ground state energy 

5.1 The regularized sum over the spectrum of hamilto- 
nian 

The problem was already formulated in Secs l2.2l and El We recall briefly 
the notations introduced in SeclHl The total energy of the system of fields is 
considered to consist of two parts 

where 

E^' = j dx£'\x) = J2j dxS'^^ix) = J2j dxTi%x) (5.2) 

i i 

denotes the complete global energy of the classical fields (backgrounds), and 

El"" = j dx£^{x) = j dx< 0\n{x)\0 >= Et (5.3) 

is the energy of quantum fluctuations in the vacuum state which is defined 
as the state without real particles. We work here with the stable vacuum, the 
problem is formulated in the infinite volume without boundaries, thus the global 
quantum vacuum energy E^ has been shown in the Sec. fl2.2|) to be equal to the 
half sum over the spectrum of eigenvalues ci;(„) of the hamiltonian Ti for the 
corresponding quantum field 

(n) 

with the symmetry (and statistics) factor C (e.g. for 4-component spinor fermionic 
field C = —4). U(^n) obey the equation 

n^nix)=ul-^^nix), (5.5) 

and by using the zeta function regularization procedure described in SecjUE'Q is 
represented the operator ?i as a generalized zeta function the operator H, ()4.115|) 

(n) 

The parameter fi appears here as a usual regularization parameter in all pro- 
cedures and serves to adjust the correct dimension of energy. 

It will be shown further, that one does not need to solve the equation (j5.5j) 
to obtain the eigenvalues explicitly and summarize them in sense of ()5.6|) . 
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By usage of the Jost function formalism developed in the scattering theory it is 
possible to avoid this troublesome procedure and express the £q^^{s) in terms of 
the Jost function for the corresponding scattering problem. 



5.1.1 Scalar field 

Consider the Klein-Gordon equation of motion for the scalar field 

{d^- A + m^ + V{x)}^x) =0 (5.7) 

with some cylindrically symmetric background V{r), depending in coordinates 
((/?, r, z) on the radial coordinate only, and use the ansatz 

^x) = (l){r)e-'^^e-'P^'. (5.8) 

In accordance with the formula for the zeta-regularized ground state energy, 
derived in Sec. IH we have 

Eris) = lf^'^j2^ki)+'^"y~"'^ (5-9) 

(n) 

— * 

k(^n) is the spatial momentum of the (n)-th wave being an eigenfunction of the 
corresponding hamiltonian. 

Because of the translational symmetry of the problem the variables z and kz 
appear in ()5.8|) only in the exponent, then it follows: 

oo 

Ens) = -^/i- I ^ Y.^'^l + kU + m't^-'. (5.10) 

-oo (") 

The symbol ^ still includes the summation over all eigenvalues of momenta 

(n) 

/c(„) transversal to the z-axis and the integration over space coordinates, excluding 
the z-coordinate: this integration has been executed and the resulting global 
energy will be understood in sense of energy density per unit length fSec ]2.2|) . 

Thus we state that the momentum k^ in (|5.10p can be separated and integrated 
out explicitly in the analytical form: 



The sum ^ means now the summation over all remaining quantum numbers. 

(n) 

In the case of cylindric symmetry these are: ki the radial momentum and / the 

oo 

orbital momentum number, then the ^ denotes ^ 

(n) l=—oo k„ 
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For scalar field the ansatz ()5.8p leads to the radial wave equation 

"I Sr^ ~^ r dr ~^ ^ '^^ ^ ' ^ ^ 

where the is in general r-dependent and denotes = k"^ ~ V{r). To proceed 
with the formalism we have to start with a discrete spectrum of momenta kn- 
To this end we suppose that fluctuations 0(r) are quantized in the finite cylinder 
of radius r = R with the suitable boundary condition. The spectrum kn is then 
discrete. 

For example we may use the Dirichlet boundary condition on the boundary 
r = R, such that ipkiR) =0. So we can consider 4>k{R) = as a function 
■= of ^) and the set of solutions fc„ satisfying this condition is the 
desired spectrum kn- Now since we know the function (^{k) we can obtain the 
sum over the spectrum kn- Note that we don't need to obtain all the values 
of kn to summarize it over the whole spectrum. Instead of that we can use 
the approach well known from the scattering theory, namely we use the integral 
properties of special construction ^ln0(A;), so called "logarithmic derivative" of 
the eigenf unction (j){k)- 

Since the eigenvalues kn correspond to zeros of (f){k) they correspond also to 
poles of the logarithmic derivative ^ In (f){k) as well, because 

and the first derivative of (f){k) is finite for all real k as supposed. The residuum 
of ()5.13p in each point kn is equal to the order of zero of the in kn- Then it 
means that the sum over kn is equal to the contour integral 

J2{kl + mY-' = ^ I dkik' + m'Y-^^lnm, (5-14) 

in) i 

where the integration contour in the complex k-plane encloses the real positive 
axis. The contour integral chosen in such a way contains all the contributions of 
residua of the function -^\n(f){k), which are placed on the real half-axis. 

The function (j){k) can also contain zeros in the k-plane aside from the bound- 
ary conditions on r = i? that correspond to the bound states and build some 
subset of the discrete spectrum. They can be taken into account by addition of 
the sum over ki 

^EE(-'-^'^)"^ (5-15) 

/ i 

to the integral over k ()5.14j) . For the further transformation of the spectral sum 
we note that the solution (pki'f') of the wave equation in cylindrical coordinates 
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(r, (fi) may be chosen by the behaviour of this solution at r — > 0. For our pur- 
poses it is suitable to require this solution to be the so called "regular" one [77j . 
that becomes proportional to the free solution (of equation ()5.12|1 with potential 
V{r) = 0). For cylindric coordinates it means 



fkTi^) ~ Mkr), r ^ 0, 
where Ji{kr) is the (cylindric) Bessel function, having the expansion 

1 /^^ " 
r(z/ + 1) V2. 



(5.16) 



JJz) 



(5.17) 



for small arguments z — > 0. 

For the solution chosen in such a way the asymptotic behaviour at r 
can be described in terms of the Jost function so that: 



oo 



1 . ^ 



r) ~ i-{mHl'\kr) + mH\'\kT)} 



(2), 



(5.1^ 



where the Hi ' [kr) are the Hankel functions, that satisfy also the wave equation 
where the potential becomes constant, (as we even suppose always for r — >■ oo ). 
These solution correspond to so called "in-" and "outcoming waves" respectively. 
The arising coefficients fiik) and fiik) are the Jost function and its complex 
conjugated. For the real k the following useful symmetry property 



h{-k) = m 



(5.19) 



applies, that is the expression of the Schwarz reflection principle known from the 
function theory and applied to the Jost function in the scattering theory 

Using the representation ()5.18|) with the property ()5.19p we can proceed with 
the transformation of the spectral sum of H expressed as a contour integral (j5.14j) 
for R —>■ oo. Rewrite the logarithm of the regular solution 0''^^ using the (|5.18|) : 



ln(f)'-^^{k) = ln{Mk)Hi'\kr) + fi{k)H\'\kr)} 



In fi{k) + In H}^\kR) + In 



In fi{k)+ In Hl^\kR) + In 



Mk)Hl^\kR) 
!i{k)H\^\kR) 



(the dropped constant — In 2 is unimportant for further calculations) and consider 
the behaviour of ()5.20|) at i? ^ oo. 
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(1 2) 

The Hankel functions ' {z) are known to behave for a large argument 

2; oo as 

^ /Ae^(-^-f) (5.20) 

V TTZ 

^ /Ae-(-^-f) (5.21) 
V vrz 

Keeping in mind these features we use the first expression of (I5.2()j) for the 
upper branch of the integration contour 7 {Qk > 0) and the second expression of 
fl5.20|) for the lower one {Qk < 0). This approach provides that the contribution 
of the third summand in ()5.20j) (the term of kind ln(l + ...)) vanishes at R 00. 
The construction of the background potential is supposed to be of kind ()H.14j] 
as considered in SecO] namely the non zero potential is concentrated in some 
environment of the coordinate origin r = and vanishes at infinity. Then the 
second term of ()5.20|) . namely the logarithm of Hankel function does not depend 
on background properties at i? — >■ cxo and can be treated therefore as a pure 
(infinite) contribution of the empty Minkowski space. It can be dropped on this 
step under this consideration. Thus we are left with the logarithm of the Jost 
function only, that contains all the necessary informations in this approach. The 
following terms remain 



The sign of the integral is changed because of direction of moving around the 
contour 7 namely 00 — > above and ^ 00 below the real axis. 

This representation can be further simplified. In this order we turn the inte- 
gration contour towards the imaginary axis. This action involves two transfor- 
mations: 

1. The upper branch of the contour turns to the positive imaginary axis: 
k ^ ik in the positive direction, so that k e*'^/^/c; 

2. The lower branch of the contour turns also to the positive imaginary axis in 
the opposite direction;since some bound states can be present, that corresponds 
to simple zeros of the Jost function on the imaginary axis. Therefore we must turn 
the upper branch of the contour to the imaginary axis in the negative direction 
to avoid the integration contour crossing the possible poles of ln//(A;) in order 
to keep the value of integral unchanged during this transformation. Thus the 
argument k must be transformed as — > e~'^'^^'^k. It corresponds to the change 
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of the integral over the lower branch according to 







-oo 
oo oo 

- j dk{e + my-^-^\nM-k) = - I dk{k' + mY-^^\nMk), (5.23) 



where the Schwarz reflection principle ()5.19j) has been used. Thus we have for 
the contour integral: 

oo oo 

j dk{k^ + mY~'\nfi{k) + J dk{k^ + m^y-'\nfi{k) \ = (5.24) 
j 

f oo oo 

-il [ dki^k"^ -m^f-'e'^^^-'Hnfiiik) - I dk{k'^ - m'^Y-'e-'^^^-'Hn fi(ik) 



Using the symmetry property ()5.19|) of the Jost function again we get for the 
integral 



COSTTS 



TT 



dk{k' -my-'lnfiiik) (5.25) 



If some bound states are present, that are placed on the imaginary axis be- 
tween and im, then by moving of the contour up to the point k = im we exclude 
the contribution of these states, that coincides exactly with the sum over ku once 
being taken into account in the 15.221 Thus the both contributions cancel each 
other and the final formula looks: 



oo 



^^^cosvrs I ^^^^2 _^2Y-s^^f^^^^y (526) 



This is a very useful representation of the regularized ground state energy. 
It means, we need only the Jost function fi{ik) (speaking exactly, the logarithm 
of it) of the problem to obtain the regularized spectral sum. Finally we have to 
recall the starting expression of the regularized ground state energy to get 



2s 



40Fr(s 



r(s-i) 

2^ l=-oo 



COS TCS 



271 



dk{k^ 



m 



2\1- 



•Infiiik). (5.27) 
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5.1.2 Dirac spinor field 



Consider the same problem for the Dirac spinor field with the radial dependent 
potential ~ a(r). To this end we recall some formulae from the Sec®. There 
it has been shown that the equation for 4-component spinor function is reduced 
to two equations for two 2-spinor functions and x 



X 



(5.2^ 



It is enough to consider here merely one of these equations , say for 0: 

^^r"" ^- l0(x) = O (5.29) 
Using the suitable ansatz 



0(x) 



V''(r)e-^'^ 



(5.30) 



we arrive at the radial equations for two-component spinor ^jj'"'' (r). 



d l-a(r) 

Po-m ^ ^ 

7^ Po + m 



V'(r) = (5.31) 



The analytical solutions of this system for the special potential a(r) ()3.21|) are 
obtained in the Section.®. Here we use the results of the previous subsection to 
construct the spectral sum via Jost function as it has been done above with the 
scalar field problem, namely we proceed as follows: 

The component Pz of the momentum is separated and integrated out again in 
the similar way: 

oo 

Ens) = -^/i^^ / ^ Y.(P' + ^(n) + ^)^-'- (5-32) 

Here E^^^^ has again the meaning of density per unit length. Taking the states 
for +pq and —po energy (upper and lower half of the spectrum respectively) and 
two different projections of spin into account we get the factor 4, since all these 
states possess the same energy. The general sign of the energy is minus because 
of the anticommutation relations for the Dirac field operators, as pointed out in 
Sec.(j2l2}. 

The explicit integration over pz/ (27r) leads to 
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Further we proceed like it has been done in |20]. A calculation for a similar 
problem has been done in [71 |T7j , however without usage of the Jost function for- 
malism. The advance of the latter approach is a possibility to gain the necessary 
uniform asymptotic expansion of the regularized energy in the elegant way by 
using the Lippmann-Schwinger equation (see Sec J5. 31 below). . 

For technical reasons we assume again that the system is contained in a large 
but finite cylinder of radius R in order to get the discrete spectrum of the transver- 
sal momenta The suitable boundary on the surface r = R conditions must 
be imposed. Though we cannot use the Dirichlet conditions for all 4 components 
of the spinor \l/(x) as in the case of scalar field because the spinor satisfied the 
Dirac equation with this condition is identical zero. One of the possible solu- 
tion is the so called "MIT bag conditions", originally proposed for bag models 
of hadrons j2H] with the constant homogeneous background field B inside of the 
"bag" 



n,. 



_d_ 

dxi 



iljijj 



2B 



(5.34) 



To proceed with these boundary conditions we redefine the integer orbital 
momentum number / of the solutions by the half-integer u according to 1)3.311) in 
Sec.©. 

Choosing the normal unit 4- vector in Minkowski coordinates (t, x, y, z) as 
= (0, — cos ip, — sin 0) and the representation of 7 -matrices ()3.24j) we obtain 
the spinor matrix n^7^ as 



-up 



-vp 







(5.35) 



Then the application of the first condition ()5.34|) to the spinor 



i)y{R,(p) 



with 



/2 



WPO + "^^?+l/2 

VPO - "^^2 1/2 

(2) 



^(i'2) = H^^2'XkR)e-'P^ 
for positive u on the cylindric boundary r = R gives 



WPO - m%+l/2 



(5.36) 



(5.37) 









ai-J 



$+(i?) 
S+(i?) 



0, 



(5.38) 
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where 



1RH^^\ ,(kR) 



and can be obtained from through the substitution m —>■ —m. The 

(Ti denotes the Pauh matrix 

" 1 
1 



0"! 



(5.39) 



and / is the identity matrix, therefore the determinants of di — / as well as (Xi + / 
are zero. It means that the first couple of 4 equations (|5.38p are in fact identical 
to each other. The same concerns the second couple of (j5.38p . 

Finally the boundary condition ()5.34l ) provides only two independent equa- 
tions : 



- 

f^ik) {y^^^T^Hf^,_^{kR) - V^^^Hl^l,_^ikR)) = 0. 

(5.40) 

Here we introduced the notations for boundary functions S^^-* for positive and 
negative values of the energy po respectively. 

To summarize the eigenvalues (m^ + k'^n)) ^^^^ boundary problem that are 
implied by the boundary conditions ()5.34|) we proceed in the similar way as in 
the previous subsection for scalar field. The eigenvalues k(n) are being viewed as 
the zeroes of the boundary functions (|5.4U|) . We use again the properties of the 
logarithmic derivative and write: 



J2{kl) - m')'-^ = / dk{e-my-^^[\nBi^^\k) + \nBi-'^\k)], (5.41) 

in) ; 

where the contributions of zeros from B^^^ and B^~^ are summarized independent. 
Finally we get for the contribution of positive orbital momenta z/, that 

E^oHs) = E / ^{k' + mr-Wn[Bi^^\k)Bi-^\k)]. 

(5.42) 



62 



As it follows from the ()3.3H) of SecOl it can be easily changed from the positive 
orbital momenta to the negative ones through the substitution of — z/ instead of 
V. Note that the Jost function in this case is different from that considered above, 
and will be denoted as f^{k). 

The change to the representation via Jost function is similar to the one de- 
scribed above for the scalar case. This procedure is described in details in the 

App.dHia) 

Now after the subtraction of contribution from empty Minkowski space (iden- 
tified as the term divergent at i? — >■ oo, which is independent of the background 
and corresponds to the contribution from the heat kernel coefficient (which 
is simple the integral over empty space), and rotation of the path towards the 
imaginary positive axis as described above one obtains pUj 



oo °° 

oo 

A merit of this representation through the Jost function of imaginary argu- 
ment is the absence of oscillations of the integrand for large arguments (in con- 
trary to the real Jost function which is known to have an oscillating behaviour)^. 



^Note that bound states are also taken into account in Eq H5.43|I . In the considered problem 
bound states appear if the flux is larger than one flux unit. Strictly speaking, these are zero 
modes located on the lower end of the continuous spectrum, i.e. at k=0 (this is known from 

El)- 
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5.2 The exact Jost functions 



The Jost function for problem with cyhndrical symmetry was defined in the 
Sec J5.1l as follows: The regular solution 

^'''{r) ~ J.(gr) (5.43) 

of radial equation is defined to behave in the coordinate origin r — > as a free 
one. The solution chosen in such a way is known [77^ to have at infinity a 
representation 

^''''{r) ~ \[f{k)Hl}\kr) + f{k)H^:^\kr)]. (5.44) 

Ju{qr), Hy{kr) are the cylindrical Bessel/ Hankel functions respectively. Then 
the coefficients f{k) and f{k) are the Jost function and its complex conjugated. It 
is assumed that the scattering potential V{r) entering in the radial wave equation 
becomes constant V{r) — > V for large r fast enough, so that the relation 15.441 is 
valid. Practical this means usually that the scattering area is localized around 
the symmetry axis of the cylindric background. 



5.2.1 Scalar field with rectangular scalar background 

The scalar background V{r) is considered to have a rectangular shape fSec ]3l3.T^ 

y(^\ ^ / ^ = co'^st Ri <r < R2 
[ otherwise 

and the solutions 4>i^i{r) , (j)jJi{r) , drCpl'^^l (r) of the equation 13. 131 are chosen in the 
form of ()3.15|) for three space domains under consideration of the potential shape. 

Now we require the wave function to be continuous and to have a continuous 
first derivative. It imposes the matching conditions on the domain boundaries, 
i.e.: 



dAi\Ar) = dAi'J {r%=R,. (5.46) 

Thus no more undefined coefficients in (j3.15|) remain and the resolution of this 
system of algebraic equations with respect to f{k) provides us the Jost function 
of the real momentum k expressed in terms of the Bessel functions J and H. In 
order to proceed with the spectral sum we need to get the Jost function f{ik) of 
the imaginary argument. The transformation k ^ ik leads to the replacement of 
the Bessel functions H and J by the corresponding modified Bessel functions K 
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and /. The resulting expression allows to be simplified by Wronski determinants 
of K and / and finally one arrives at 

fi{ik) = R1R2 X (5.47) 
{ [qIi{kRi)K[{qRi) - kl[{kRi)Ki{qRi)\ [kIi{qR2)K[{kR2) - ql[{qR2)Ki{kR2)\ 
- [qIi{kRi)l[{qRi) - kl[{kRi)Ii{qRi)\ [kKi{qR2)K[{kR2) - qK[{qR2)Ki{kR2)\} 

Here primes denote the derivation with respect to r. It is suitable to check 
the useful symmetry property of the Jost function . Taking in mind that the 
/ is integer and using the explicit symmetry relations for the modified Bessel 
functions K and / [lU] : 

Ki{z)=K.i{z) (5.48) 

we can see that 

fi{ik) = f.i{ik). (5.49) 
This relation allows to consider further the Jost function only for / > 0. 



5.2.2 Spinor field with rectangular magnetic background 



The radial Dirac equation ()3.28|) corresponding to this problem with the chosen 
special form of background ()3.21l3.22j) . has been already solved in the SeclHl The 
solutions ^iii in three space domains (jSI21 - EiSSl) are obtained in terms 

of Bessel and confluent hypergeometric functions . 

To get a Jost function we need to impose certain matching conditions like in 
the scalar problem considered above. It is known that for a continuous potential 
a(r) as we consider here, the spinor wave function must be continuous too. Thus 
for the spinor wave function ()3.2I - I3.33|) we demand only to be continuous at the 
boundaries between I-II and II-III domains. It means: 



r=Ri 

r=R2 
± 



(5.50) 



Resolving these equations with respect to we obtain, that 



A 
B 



o A + B 
WniWn 



Ri 



II.r> 



l± 



R2 
R2 



(5.51) 
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The denominator of this expression can be written using the Wronskians of 
hypergeometric and Bessel functions as follows 



Win = WnAii.r - ^ni.rin.r) 



,Po + m 1 4i 

j2 



R2 po — m 7ikR2 



W, 



II 



Wl.r^n. - 4l.r^IU 



a/3 ,I3R 



1 



Po — m 



(5.52) 
(5.53) 



To calculate Eq we need the Jost function of imaginary momentum ()5.43p . so 
we need to replace k by ik. As a result we obtain the new expression for f{ik) 

that contains now modified Bessel functions instead of J„_i_i and modified 

'^^2 2 

Bessel functions Kj^_^i-^g (McDonald functions ) instead of 



TikRiR2e~^/''^ 



X 



2v+l 



29 + 1 



iFif z/ H , 

R2 4 2/?' 2 



2/3' 2' 

fc2 1 



X 



/3Rl 



K.sikR2)- 



{kR2 



/c iFi z/ H -, u\ 

^2 2/5' 2 ' 2 ■ 

^ 1^1(2 " + 2;^' 2 ~ 2 ^ h+^S^Ri) + 



2^^^^^^+2:^' ^^+2' 2 



3 (3Rl 
k^ 



K.,{kR2) + 



1 PRh^^ 



{kR2 



where v = v + -sr-fe, and 
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TikRiR2e~^/^' 



kRi 
2v-\ 



+ V 



iFi(z/ + - + 



(5.54) 



1_ ^R\. 



irRi) 



1 e 3 (3Rl, ^ 



{kRi 



X 



2z/- 1 

i?2 



lFl(l 



^ 3 
2^' 2 



_ (3Rl 



k iFif— -, z/: 

^2/5' 2 



PRl 



{kR2 



2z/ 



(3Rh 
2 ' 

1 A;^ 



(A;i?2)+ 



3 pRl. 



{kR, 



X 



where v 



V 



SRl 



5.2.3 Spinor field with background of Nielsen-Olesen Vortex 

The technical distinction of this case from those recently considered above is the 
not explicit background function , given by two functions /(r) and v{r), being the 
solutions of the system of two second-order equations ()3.49|1 - Nielsen-Olesen (NO) 
equations . This system has no known solutions representable in analytical form. 
Therefore one cannot apply the experience of two previous models to obtain some 
analytical expression for the Jost function . Nevertheless the system of Nielsen- 
Olesen equations is known to be solvable at least by numerical methods. Let we 
have the radial part of 2-component spinor function $i,A:(?") to be a (numerical) 
solution of the system 



d_ 

dr 



l+l—nv(r) 
r 



(5.55) 
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with 

$(r) ={ V , (5.56) 




indices I and u denote "upper" and "lower" components of spinor respectively. 
These equations can be solved numerically for any specified functions v{r) and 
fir). 

In order to define the Jost function properly we start as usual with the defini- 
tion of the regular solution for the system ()5.55p . The behaviour of the potentials 
/(r) and v{r) is known for r — and r oo, since they define a vortex solution. 

The function ^f{r) entering here plays a role of space dependent fermion 
mass and for further purposes it will be denoted as /u(r). 

The usual boundary conditions for ()H.28j) are generally speaking 

f{0)=v{0) = 0; /(oo) = t;(oo) = 1; (5.57) 

for the Nielsen-Olesen vortex. It corresponds to the situation that the fermionic 
field in the center of the vortex become massless ( ^fir) = fJ-ir) = 0), and 
outside of the vortex (far enough from the vortex core) has the observable mass 
rrte. In order to check the symmetry properties of solutions better, we suppose 
for a moment that 

/i(0)=mo, (5.58) 

mo is some constant, that will be later set to zero. 

The behaviour of the f{r),v{r) for vorticity n = 1 closed to zero follows from 
the corresponding asymptotics of the NO-equations ()3.49p and is known to be 



0. (5.59) 



/(r) ~ Ji{r]VXr) 
v{r) ~ 

Further we must take respect to this fact, that the fermionic mass distinguishes 
for r — i> and r ^ oo and are equal to ttiq and rrif. respectively. Note also that 
< rriQ < nif,. Thus we have to define two radial momenta 



q = .ipl-ml 



k = y pI — ml. (5.60) 

Now the definition of the Jost function demands some generalization. Here 
we follow the approach derived in j^. The regular solution of ()5.55|1 will be 
chosen to be proportional to 



f VPo + TTT'O Ji-n+i{qr) 
^-^(r)~(-) { } (5.61) 

y/po - mo Ji-n{qr). 
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Further we represent the regular solution ip^'^^[r) at arbitrary r 7^ in the 
form: 

mr)oo = I {/KA:,r)V^(A;,r)« + /.(A;, r)V^(fc, r)(^)} , (5.62) 

where 



(1 2) 



|(1,2) 

'00 



(5.63) 



Here the coefficients fi depend now in general on the radius r of the cylin- 
drical shell. The background potential enclosed inside of this is considered as the 
scattering potential. 

Then it follows from the asymptotic behaviour r ^ 00 of the proffie functions 
/(r),t>(r), that for the radius r large enough the functions fi{q,r), fi{q,r) do 
not depend more on r and coincide therefore with the Jost function fi{k) and 
its complex conjugate fi{k) respectively for this scattering problem. Thus we 
define the function fi{q,r), which for large r approximates the Jost function . 
The detailed motivation of this assumption is given in |21j . 

The / is replaced up from here by z/ as usual by means of ()3.311 SecOJ. Thus 
for i> = 1/2, 3/2, 5/2, ... we have two solutions for positive and negative orbital 
momenta, denoted below as i>+ and i>— respectively. 



1 



/^(fc,r) 



+ 



VPo + me Hl^li/2-nikr) 



y/PO - nie HI 



(2) 



(kr) 



ipu-{k,r) 



fuik,r) 



VPo + me H^ll+i/2-Skr) 



(kr) 



+ 



L-{k,r) 



VPo + H^^l+i/2-Skr) 
-VPo - me H^^l^y^^^{kr) 



(5.64) 



(5.65) 



The spinor functions ip{k,r)'^^ are the (numerical) solutions of ()5.55p . 
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Thereby the generahzed (r-dependent) Jost function by means of ()5.62|) can 
be calculated as follows 



-nrr 



(5.66) 
(5.67) 



r(l) 



with the initial condition 



kV^-- 



VPo + mo J^±i(gr) 



(5.68) 



r=0 



The spinor upper and lower components functions il)^^''^^{k, r) obey the equa- 
tions : 



Po - ^^{r) 

Q v+\-nv(r) 



Q i/—^—nv(r) 
dr r 

Po + Kr) 

Q —v—^—nv{r) 



Po - f^ir) ^ — 

\ Po + /W(^) 



dr 




(5.69) 



However we are interesting again for the Jost function of the imaginary mo- 
mentum i/c, so we make the transformation k —>■ ik, like the one considered above 
for magnetic background. By means of ()5.60|) this implies 



Po P = ^/ m'^ — k'^ = i a/A;^ — = ip 



q = \J k"^ — + ml <ij ml — k"^ — uIq = i^J k'^ — [ml — m^) = iq, (5.70) 
resulting finally in 
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/+(zA;,r) = -{-T)-\-'-r | ^ip - m,K^_._Skr) C+(2A;,r)+ (5.71) 

i^ip + m^K^^^_^{kr) ^^+(iA;,r)| 

/-(zfc,r) = {-ir+'^^-r[^/W^eK^iUkr) r.-{^Kr)+ (5.72) 

i^ip + meK^_i^^{kr) i)[-{ik,r)^ 

The functions il)lf''^^{ik,r) entering here obey the modified Dirac equation 
which is to obtain from the original one (j5.69p through the replacement of Pq by 
ip according to (j5.7U|) . 

From the ()5.73|) together with the equations ()5.69|) modified by means of 
(I5.70|) follows in particular, that for the Jost function of imaginary momentum a 
change of the sign of scalar potential /i corresponds to complex conjugation and a 
change in the sign of the magnetic background v{r) to opposite one v{r) —v{r) 
corresponds to an exchange of positive and negative orbital momenta f^^f~- 
It leads to the simplification of the resulting formula for the regularized ground 
state energy , which includes the sum of logarithms of and /~. The odd 
powers of v{r) disappear and pure imaginary contributions of as well. These 
features are described and argued in details in [14 . 

Using these formulae we could be able to calculate the needed Jost function 
numerically. To this end we must have numerical solutions ip'^'^{ik, r) of the Dirac 
equation with the profile functions f{r),v{r) which are calculated from the NO- 
equations ()3.49|) numerically as well. 



Numerical calculation of f^{ik,r) 



The practical experience shows the need of higher raised precision to get 
correct results for ip'jf'^^lik, r) by usage of "NDSolve" package of " Mathematica" 
for the radial Dirac equation directly in the form 



+ Kr) 



dr 



and starting with the initial condition following from (jf 




y/p + mo 

tWP - ^0 



(5.73) 



(5.74) 



r=0 
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To avoid this trouble we apply the representation of regular solution at the 
point r = ()5.74|) in the form 



C^(fc,r)=/o(z/)C±(A;,r), 
where Jq are defined from the asymptotic behaviour of 



I.{kr) 



for small arguments, namely: 



{kr 



kr 



2u+ 1 

Thus we can rewrite the initial condition as 



ip {r)\r=o 



-iIo{u + 1) 



^+(r)|.=o 



iloiu + 1) 



(5.75) 



(5.76) 
(5.77) 

(5.78) 



Then for r ^ the initial condition for numerical iterations of the system 
()5.73|) in terms of (r) looks as 



^^(r)|.=o 



(5.79) 



Now the final formulae to calculation of the Jost function in terms of the 
redefined wave functions looks 



f+{tk,r)= (5.80) 

-/o(z/)(-z)'^~^'V I y/ip-m^K^_i^{kr) ^\ + i^/ip + meK^^^{kr)^^+ 

f~{tk,r)= (5.81) 

which are to evaluate numerically. 
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5.3 The uniform asymptotics of the Jost function from 
the 

Lippmann-Schwinger equation 

The present section contains the exphcit derivation of the uniform asymptotics 
of the Jost function that is necessary to construct the regularized ground state 
energy E^*^^ has been used in the previous Section ()5.4|) 

The uniform asymptotics of the Jost function fu{ik) is the power expansion in 
the combined hmit k,i'—>-ooa.tk/i' = z = const. It can be obtained in various 
ways. One can make for example the explicit substitution of asymptotic expan- 
sions for Ki{kr), Ii{kr) and other functions that the Jost function fuiik) consists 
of. Another approach is the iteration of the Lippmann-Schwinger equation for 
the corresponding scattering problem. 

This procedure is the most laborious analytical topic of the method. 

5.3.1 The Lippmann-Schwinger equation 

The Lippmann-Schwinger equation is one of important tools of the scatter- 
ing theory. This is simple an expression of that fact, that the scattered wave 
function (j){x) = 0°"*(a;) observed far enough from the scattering area can be rep- 
resented through the initial (unscattered) wave function (p^ = and integral 
of interaction with the scattering potential. 

It can be shown to follow from the theory of ordinary linear differential equa- 
tions : 

We suppose to be the general solution for the homogeneous equation 

(or "free solution" by notation in scattering theory) 

P0O(x) = (5.82) 

and 0^(x) to be any particular solution for the inhomogeneous equation of kind 

P(j)^{x) = V{x)(j)^{x) (5.83) 

with the V(x)- generally speaking some function of coordinates x, in our case it 
plays the role of the potential. Then we have that the general solution (j){x) of 
the inhomogeneous equation (j5.83|) reads 

(f){x) =(f)\x) +(j/{x), (5.84) 

in accordance to the statement about general solutions of inhomogeneous equa- 
tionEia (e.g. IHOl). 

On the other hand the particular solution 0^(x) of the inhomogeneous equa- 
tion EUSl can be expressed by using the Green function g{x,x') of the operator 
P, satisfying the definition 

Pg(x,x') = 5(x-x'). (5.85) 
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Suppose such a function is known. Then the solution of any inhomogeneous 
equation ()5.83|) can be restored from the inhomogeneous part (right-hand part) 
by the integration 



<P''{z) = J g{x,x')V{x')(j){x')dx'. (5.86) 



Thus we obtain the integral equation for the general solution (j){x) 

X 

(j){x) = (p^ix) + / g{x,x')V{x')(j){x')dx', (5.87) 





that calls in applications to scattering problems "Lippmann-Schwinger equation" 
We consider below the several cases of usage of the equation l5.87l for potentials 
V{x) with cylindric symmetry. 

5.3.2 Scalar field 

The equation of motion for the massive scalar field ip{x) of the mass m in 
the scalar background V{x) depending in cylindrical coordinates on the radial 
variable r only, has been reduced to the one for radial part (SecOj) 

|^ + -^+ffc^-^)U(r) = 0. (5.88) 
ar^ r or \ r'^ J ) 

Then the Lippmann-Schwinger equation ()5.87|1 reads: 

oo 

Uz) = <P%z) + j dz'Gi{z,z')V{z')<f)i{z'), (5.89) 



where the Green function Gi{k,r,r') of this equation is constructed according to 
the known prescription [ZZl EOj as 

G,{z, z') = '-^[Mz)hI'\z') - Mz')hI'\z)], (5.90) 

and z denotes z = kr. 

It can be easily proved by explicit substitution that ()5.90|) satisfies the defi- 
nition of the Green function ()5.85|) apphed to the equation ()5.88|) . 

Further, substituting the Green function (j5.9(J|) into the Lippmann-Schwinger 
equation (j5.89|) and choosing the solution 0°(-2) of the homogeneous equation (i.e. 
at V{z) = 0) as 

= M^), (5.91) 



74 



we have for the Lippmann-Schwinger equation 



(f)i{kr) = Ji{kr) + — J r'dr'V{r')[Ji{kr)H\'\kr') - Ji{kr')Hl'\kr)](l)i{kr'). 



(5.92) 

Now we need to obtain the expression for our required Jost function fi{k) 
defined in fj5.93p . Let the solution of the inhomogeneous equation (f^ {kr) 15.831 
far from the scattering area be denoted as (f^^^ikr). Since the solution at r — > 
has been chosen to be a "regular one" fSec ]5.2l Eq ]5.43n . it can be represented 
at r ^ oo (SecEH Ea lOil) as 

= \[mH['\kr) + mH'i''\kr)]. (5.93) 

Substituting this relation into ()5.92j) and using the certain relations between 
Bessel and Hankel functions ( Ji and H^'^), ^\ 

Mz) = l{Hi}\z) + Hi'\z)} (5.94) 

we arrive at the relation 



1.^ 



^[Mk)Hi'\kr) + Mk)Hl-^>{kr)] 



r(2)/ 



Hl'\kr) 



Iff (2) 
2 

lH;'\kr) 



1 + y / r'dr'V{r')H[^\kr')(t)i{kr') 



+ 



1 — — j r'dr'V{r')Hi'^\kr')(f)i{kr') 



whereby the Jost function fi{k) is left to be identified as 



fi{k) = 1 + ^ y r'dr'V{r')Hl^\kr')(j)i{kr') 



and for fi{k) one obtains 

oo 

f^(k) = 1 - ^ y r'dr'V{r')H[''\kr')<Pi{kr') 



(5.95) 



(5.96) 



(5.97) 



respectively, that is really not surprising, and confirms the supposition once more, 
that fi{k) and fi{k) are indeed complex conjugated to each other, what was not 
initially adopted. 
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Going on with the regularization procedure we have to transform the ob- 
tained real Jost function fi{k) to the complex one, changing the argument k to 
ik respectively. It provides [IHl [l] 

oo 

fi{ik) = 1 + y r'dr'V{r')Ki{kr')(t)i{ikr'), (5.98) 



where K[^\kr) is the modified Bessel function (McDonald function ). The wave 
function (j)i{ikr') entering here obeys now the equation derived from ()5.92|) by the 
substitution k ^ ik 



(f)i{ikr) = i^Ii{kr) + j r'dr'V{r')[Ii{kr)Ki{kr') - Ii{kr')Ki{kr)](j)i{ikr'), (5.99) 



{Ii{z) is the modified Bessel function , [13 CI)- On comparing ()5.98p with ()5.99p it 
can be seen that the factor i^' provides only the redefinition of the wave function 
and is irrelevant for further calculations. We omit these factors in both equations. 

Solving the equation ()5.99|1 by iterations of (f)i{ikr') and beginning from the 
(f)f'\ikr) = Ii{kr) as a zero order approximation of 0/ one obtains after two 
iterations: 



(j)[^\ikr) = Ii{kr) + j r'dr'V{r')[Ii{kr)Ki{kr') - Ii{kr')Ki{kr)]Ii{kr') 



r 

<j)f\ikr) = Ii{kr) + / r'dr'V{r')[Ii{kr)Ki{kr') - Ii{kr')Ki{kr)]Ii{kr') + 



r'dr'V{r')[Ii{kr)Ki{kr') - Ii{kr')Ki{kr)] x (5.100) 



r"dr"V{r")Ii{kr')Ki{kr")Ii{kr") - I r"dr"V{r")Ki{kr')Ii{kr"f 



The considered regularization method doesn't operate with the Jost function 
immediately, but with it's logarithm In fi{ik) (See Sec 15.11) . Then by using the 
standard series representation of logarithm: 

ln(l + x) = X - -x^ + -x^ - ... (5.101) 
2 3 

we can derive from ()5.98|) by means of ()5.101|) that 

In fi{ik) = In f[^\ik) + In ff\ik) + In ff\ik) + In ff\ik) + ..., (5.102) 
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where the terms In //^■* (ik) denotes the resuhing expressions, that after substi- 
tution of (|5.1(J1|) into (|5.98|) and re-expansion and re-ordering of terms contain 
integrations n times respectively (or the same, the ra-th power of V). They are 
represented in Math.App. ()8.1.3p . 

The given approach is though apphcable for all models considered in the 
present work. However it will be shown later in Sec 15. 41 that for regularization 
procedure the only minimal subtraction scheme is used. So for the scalar field we 
need only the contributions up to In fi'^\ik) to take into account in In/"*. 

The procedure of simplification of In //^'' (ik) fro this case is described in details 
PU] and partially in Math.App. (|8.1.3p . Finally we obtain the result for / > in 
the form 

fimr^^ I rdrJ^Xn/^ (5.103) 

ri 

in terms of polynomials of the variable t = , , Here X„ , are: 

Xi,i = l, X3,3 = i(l-2rVo), 

^3,5 = -f, Xsj = l- (5.104) 
For / = we have additionally at oo 

/o(^A;)~-^(ri-r2), (5.105) 
the derivation is presented in Math.App. ()8. 10118. l2j) . 

5.3.3 Spinor field in cylindric magnetic background 

To obtain the uniform asymptotic expansion of the Jost function , we proceed 
in much the same way as in the previous section for the scalar field. 

It starts again with the integral Lippmann-Schwinger equation for the general 
solution 0(r) 

r 

0(r) = 0°(r) + J g(r,r')V{r')(f)(r')r'dr' , (5.106) 



where the 0(r) is a two-component spinor function , satisfying the equation 

P0(r) = V{r)(l){r) (5.107) 

with 
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The Green function g{r,r') is chosen in the form 



g{r,r') = -^[0°,(r)0?,^(r') - /.(r'X^W]- (5-109) 

Here the 0j(r) and 0j(r) are solutions which behave at r — ^ as ones of the 
free equation , it means the equation ()5.1U7p with ^ = 0: 



«('■) - ( ) for I < 0. (5.111) 

and the corresponding 0^^(r') is obtained from 0j(r) by the substitution Jp{kr) 
through H^p\kr) and transposition. It leads to the equation for the real Jost 
function fi{k) 

oo 

fi{k) = l-^J r drct>Y{r)V{r)(t>{r), (5.112) 



where 0(r) obeys the full equation ()5.107|) . 

The orbital quantum number / is substituted by u by means of ()3.3H) . Sec 
El and the radial momentum k is turned to ik respectively. In contrary to the 
recently considered scalar field problem, in the case of spinor field we have to 
retain the terms up to In/i^^ in the expansion (|5.1U2l / z/), in order to take 
into account all powers of v up to needed for the minimal subtraction. This 
procedure can be look up in details in [201 especially for the problem of the 
spinor field in magnetic background. (The similar calculation is given in the 
present paper in details for the spinor field in background of NO-vortex, see the 
next subsection and Math.App l8.1.5j) . 

As a result, the expansion in powers of v for logarithm of asymptotic Jost 
function of imaginary momentum is obtained in the form (see|20j) 

^ r A J. fj 
n=l j=l -l 

where Xnj are represented in terms of r, f2, a(r) and their derivatives: 
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^l,l = -Xi,3=X2,6 = ^(afi)2 

X2,2 = ^^^(ct^ ~ raa ) 
X2,4 = ^Q'^{-3a'^ + raa) 
3 3 = -ii (a — raa + -r aa — -\l a 



■'"02/2 1.^2 II -'-02 4\ 
-\l [a —raa H — r aa \l a 

X3 5 = — —a^ + 7raa' — r'^aa" + 6i7^a^ 



4 ' 2 2 

lo2/ ^„2 , v^„„/ _2„„// I cri2„4N 
2 



—35 

^3,9 = ^^'a', (5.114) 



fcrl2\-^ 



This expansion had been obtained in ^H] by iterations of Lippmann-Schwinger 
equation up to the order z/~^. In general it is possible to obtain higher orders 
using this formalism. However as the calculations [12] showed, the complication of 
the involved expressions increases very fast. It is remarkable that this expression 
does not contain a term with power z/~^. In the finite part of the energy the 
corresponding term is cancelled in the sum of terms corresponding to positive 
and negative orbital momenta u as well. The absence of the power is a 
succession of the zero heat kernel coefficient a3/2- Also it is a non-trivial fact that 
both, the fourth power of the magnetic flux Q and the second one is present. 

It can be checked numerically that Zn/"* (j5.113|) is indeed the uniform asymp- 
totic expansion of logarithm of (j5.54l I5.55j) for k,^- large, k/u = 2;-fixed. The 
numerical examples will be shown in Sec®. 

5.3.4 Spinor field in arbitrary abelian Higgs background with 
cylindric symmetry 

As it was pointed out in SecOl the abelian Higgs background distinguishes 
from the case of pure magnetic background considered above, since the spinor 
field has a space-dependent mass, coming from interaction with the scalar Higgs 
condensate. It gives reason to discuss this problem more detailed, following in 
principle [H]. 

We start with the Lippmann-Schwinger equation in the form ()5.87|) . and in 
order to proceed with the expansion for Jost function in the form ()5.1021 we have 
to define an expression for the Green function in ()5.87p in a suitable form. 

Consider the equations of motion for the spinor field In SecOthey have 
been reduced for the special case of cylindrical symmetry in cylindrical coordi- 
nates (r, <f, z) to the form 



79 



Po + L- mas Ps^^s 

PsCTs po + L + mas 




(5.115) 



with 

2 g 

L = iJ2(^k{-g^ + ieAk) (5.116) 
fc=i 

for a wave function ip{r) of the radial coordinate only. Now let the potential V{r) 
be the solution of the Nielsen-Olesen vortex equation, (that we namely need) and 
is expressed in terms of functions v{r) and /(r) respectively. 

For the radial equations for two-component spinor $ = [tp^ {r) , tp^r)) we 
obtained: 



_d_ _ l+l~nv{r) 
dr r 



with 




$(r) = (5.117) 



$(r) = <^ } . (5.118) 



Indices I and u denote "upper" and "lower" components of spinor respectively. 
These equations can be solved numerically for any specified functions v{r) and 
/(r). The integer n is the vorticity of the vortex configuration (SecOJ. AH further 
calculations are given for n = 1. 

For boundary values of potentials /(r) introduce the notations: 

-I-f[r)=mQ,r^O -^/(r) = m^, r ^ oo, (5.119) 

and for the v{r) potential is still valid 

v{r) = 0, r ^ v{r) = 1, r ^ oo. (5.120) 
Further we introduce the notation for scalar potential 

^/(r) = Mr). (5.121) 
Now consider the construction 

g{ry) = -^^mr)^Vir') " ^'j{r')^Vir)}, (5.122) 



80 



where 



^Po + mo H[]^l\pr) 



$0 



(5.123) 



with p = a/po ~ ""^0' assure, that it satisfies the definition of Green function 
for the equation ()5.117j) at the hmit of r — > 0, which we assume to be a homoge- 
neous equation by means of (j5.82p . Sec ]5.3l Here we must take into consideration 
that the integral over r' in the general solution for the inhomogeneous equation 
f|5.117p in view of ()5.83p . Sec ]5.3l is to be taken 



r 



g{r, r')\>(r')$^'''(r'), 



(5.124) 



in accordance with ()5.87|) . Sec 15. 31 where the integration is to carry on in limits 
from r' = to r' = r. This is equivalent to the choice of the Green function in 
the form 



g{r, r') 



T/ 
H K 



(5.125) 



Basically the method presupposes, that the Lippmann-Schwinger equation 
does not "know" about the behaviour of the solution of the ()5.117|) outside of 
the close neighbourhood to the point r = 0, where the functions ()5.132j) satisfy 
the equation , respectively thereby the constructed Green function (j5.122j) as 
well. These informations are all to be contained in the integral kernel V{x) of 
()5.124p that is even the scattering potential of Lippmann-Schwinger equation . 

Then one can proceed with iterations of the Lippmann-Schwinger equation 
in the same way as it has been done for the spinor field in magnetic background 
(previous subsection) and represent the logarithmic Jost function immediately as 
the series 



InMtk) = \nfj^'\tk) + \nff^{tk)+\nf^'\tk)+\nfl'^{tk) 



<-{2) 



(3), 



P(4)/ 



(5.126) 



Here, as well as in the previous case, we need at least four initial terms of the 
uniform asymptotic series for In fuiik) in order to keep all powers of u up to z/~^. 

In order to use the Lippmann-Schwinger equation properly we solve the equa- 
tion of motion ()5.117|) with a little bit generalization. Suppose, the scalar poten- 
tial fi{r), representing the local space depending mass of the spinor field $(r) to 
acquire the value mo (not necessary equal zero) at the coordinate origin r = 
and me outside of the vortex core at r ^ oo respectively. 



81 



Here we must handle the solutions of ()5.117p carefully. First of all rewrite 
this equation for spinor function $ as 



{P-V}<i> = 0, 



(5.127) 



with the separated "source part" V which contain the dependence of the poten- 
tials /i(r) and v{r). Introducing the redefinition of fi{r) 



fj,{r) = fi{r) - mo + mo = jl{r) + mo 
we rewrite the fl5.117p in form ()5.127|) with 

Po-mo dr- ^ 

-dr - ^ Po + mo 
-7 /^W 



(5.128) 



V 



(5.129) 



(5.130) 



Then we can assure that for vanishing matrix potential V the following spinor 
functions 



H 



^ 2 



(5.131) 



(5.132) 



with go = a/po ~ ""^0; really satisfy the homogeneous "free" equation at r ^ 0. 
Indeed the 2-component spinor functions satisfy exactly the equation 

fl5.127|) without the right-hand side, i.e. at = 



P^u = P^j = 0. (5.133) 

Thus the form ()5.127|1 with P and V defined by ()5.129^ I5.13()j) is obviously 
appropriable to carry on iterations of Lippmann-Schwinger equation . Further 
steps are represented in Math.App 18.1.51 The final expression for In f^{ik) can 
be represented again in the form 



\nf^{ik) 



dr 



3n 



n=l j=l 



P 



nj 



V" 



(5.134) 
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Xn = ^v{ry + ^r\fi{rY -ml), 

X33 = ^v{r f -^r^v\rf -^v{r)^ -^r^v{rf{fi{rf - ml) 

39, .9 "^2//\2 ^/\4 ^2/\2//\2 2\ 

^35 = -v^v{r) + -r V {r) + -v{r) + -r v{r) {^{r) - m J 
lo 8 4 4 

3 2 / / \2 2\ 

-— r (/i(r) -mj, 

Xs7 = fvir)^-lv{r)^ + f^r\^,^-ml), 
3^ 

X39 = -777^(0'- (5.135) 
Id 

The dependence of Xnj on u(r) is the same as in the previous problem of mag- 
netic flux tube (obtained in |2I3, where it was denoted by a(r)), the dependence 
on n{r) is a new feature. We use integration by parts in Eq.f I5.l5^ in order to 
get the shortest representation for the coefficients Xnj. 
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5.4 The representation of the renormahzed ground state 
energy 

The results of Sec ]5.1l allows to represent the regularized ground state energy 
in terms of Jost function for corresponding scattering problem. These repre- 
sentations ()5.4H[ I5.27|l for each of considered models are constructed as a sum 
over orbital quantum numbers of integrals over radial momenta from logarithmic 
derivative of the Jost function . The infinite contributions of the empty Minkowski 
space are dropped through this performance. However the regularized energy still 
remains infinite through the ultraviolet divergences. 

Further in the Sec J4.5l the asymptotic expansion for the regularized ground 
state energy in terms of global heat kernel coefficients was obtained ()4.119|1 . 
It corresponds to the asymptotic of large mass m — > cxd for massive quantum 
fields. It was also pointed out that this expansion allows to identify all ultraviolet 
divergences of energy. 

Proceeding from the ()4.119|1 we define the terms with non negative powers 
of m (all the contributions up to 02) as divergent part i?'^™ and subtract it out. 
Thus the renormalized energy is defined 

^ren ^ ^reg _ ^div ^ (5.136) 

that now is finite, since all the present ultraviolet divergences are contained in 
the E"^™. 



^atv ^ _ + In ^ - - Uo ^ai/2 + 

/ 1 4/i^ \ m 



32n'^ \s rn? J IGvr^/^ 

It should be noticed that such a subtraction removes not only divergences 
being contained in the heat kernel expansion but also several finite contributions 
not related immediately to the ultraviolet divergences at (s 0). This action 
is dictated by the physical requirement, that the energy of quantum fiuctuations 
must tend to zero for large mass limit. This approach differs from the minimal 
subtraction (j4.1U9|) and is valid for massive fields only. 

5.4.1 The divergent part of ground state energy 

For certain models considered here the coefficients a„ entering in ()5.137p and 
therefore the divergent part of energy are calculated. 
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1. Scalar field with the scalar step-function background V{x) 

The problem is formulated in SecOl and the corresponding non zero heat 
kernel coefficients ai,a2 are derived in Sec J4.H^ Eg. ()4.52^ I4.52|l . Thus all the 
coefficients a„ in fl5.1H7j) are known: 

flo = 1 

0-1/2 = Os/2 = 

ai = 2n J —V{r) r dr 

a2 = 2n j (^^/\V{r) + ]^V'^{r)^ r dr (5.138) 

For the considered potential V{r) SecISl ()3.14|) the relevant coefficients con- 
tributing in E'^'^'" are 



ai = W^{Rl-Rl), 



a2 



(5.139) 



what provides 



327r 



s / V s 



(5.140) 



2.Spinor field with vector (magnetic) background A^{x) 

The heat kernel coefficients for the Dirac spinor field in background of vector 
filed are calculated in Sec 01 Eq. ()4.72|) . In fact only the coefficient 02 contributes 
to the E"^^": 



ao = 1 

^1/2 = ^3/2 = ai = 



(5.141) 



02 = Tr 



d^x 



-—F.^F'"' + -(r'"'F, ^2 
12 8 



\ I d^xF.^F'^'' 



(5.142) 

For the magnetic rectangular shaped tube ()3.2()f3.2H) a constant cylindric 
magnetic fiux Q SecEl ()3.3()j) it means: 



02 



d'xF^ 



Sir 



n"^ / drrh{ry 



3 i?2 — -^1 



2 • 



(5.143) 
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Moreover in the section ()5.6|) the next non-zero heat kernel coefficient 05/2 
has been calculated (that however does not contribute to the divergent part of 
energy ) 



05/2 = ^^e^^ 



2n2 (-^1 + -^2) 



^2) 

This coefficient (and further half-integer numbered coefficients) is present in 
the heat kernel expansion because of non-analyticity of the background Air) (the 
remark at end of Sec l4.3|) . 

The corresponding i?*^™ for this problem is therefore 

fi2 1 4„2 

^ = , ,02 p2, (- - 2 + In ^). (5.144) 
67r(i?2 — R{) s 



3. Spinor field with background of the Nielsen-Olesen vortex 

The only non zero coefficients contributing to the divergent part of the ground 
state energy are (fT^ : 

a2 = Tr j £x (^-^F^^F'^'' + ^V{x)^ - ^AV(x)^ , (5.145) 

where the trace is to take over spinor indices. 

All the Qn with half-integer n are zero since the background potentials f{r), v{r) 
are supposed to be smooth (the solutions of the NO-equations). 

In terms of yu(r) = ^/(r), v{r) the heat kernel coefficients read: 





Thus the E'^'" reads: 



ai = —Stt j r dr{(x{r) —m^), 



00 

a2 = 87r |rrfr|i^ + i[^'(r)2 + (Mr)^-m,Y]| (5.146) 



^<l^v ^ _ ; _ + In ^ _ 1 r dr^ifir)' - 1)+ (5.147) 



1 , 4/i 



2 



00 



+ In — - — 2 I r dr 
s J 
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5.4.2 The finite and the asymptotic parts of energy 

E^^^-^s) defined by ()5.136|1 is though finite at s — > 0, but it is not possible to 
carry out this limit analytically. Indeed, we have merely the finite sum of two 
terms, each of them is infinite and represented in distinguished form: the first one 
is an infinite sum of infinite integrals, and the second one is generally speaking 
an algebraic polynomial of powers s~^. The possible way to solve this problem 
would be to find some intermediate object possessing both the features of two 
summands in E^'^"', i.e. it should be on the one hand a sum of infinite integrals, 
and on the other hand it must be a polynomial of s containing in particular the 
powers s~^. 

Furthermore we know that the divergence of the expression for E^"^^ is of the 
ultraviolet type. In terms of variables / and k (orbital and radial momenta) it 
means that the divergence appears at large values of them. It suggests that the 
needed object is some kind of asymptotics for k,l oo 

One can construct such an object proceeding from the structure of Eq'^^ . The 
following discussion is restricted on the case of scalar field. It will be applied 
later for other models. 

Let we have an expression of kind ()5.43|1 . Consider the orbital momenta / large 
enough, that it can be treated as a continuous variable. Suppose the function 
\n fi{ik) entering here has an asymptotic expansion in the combined limit 

k oo, / oo, k/l = z = const, (5.148) 

the so called "uniform asymptotic expansion", known e.g. for Bessel functions 
III 

Then the Jost function is expected to have the asymptotic behaviour 

/K^A;)~l + ^ + | + ..., (5.149) 

that follows from the general properties of the Jost function [77j. For In fi{ik) we 
have under assumption of ()5.148p 

1 fCJ.^ C,{z) C,{z) C,{z) 

\nti[tk)r^ — ^ \ \ p h... (5.150) 

Substituted into the expression 

E^'s^CY. / dk{k''-my-'—\nfi{ik), (5.151) 
the ()5.148p provides 

oo 

E^^^^CsY.1'^'-'^ I dz{z'-'^f-^(^ + ^ + ^ + ..\ (5.152) 
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Now, as far as the parameter s is not fixed, we suppose for a moment that s is 
large enough that the integral in ()5.152|) converges. Then we can interchange the 
integration with the sum in brackets. It allows to establish that the first three 
terms in the regularization limit s ^ 0. 

CsJ2l'^'~^^ [dz{z^-^r^(^), j = l,2,3 (5.153) 

can be divergent, and for all further ones up from j = 4: the sum over I converges. 
We denote the divergent terms of the uniform asymptotics (j5.153|) as S^^^"^ 

oo 

oo „ 
gasym = C,J2 ^'^'"'^ / " 

oo 

°° f f) 

C^Y, dk{k'-m'y-'— In frm, (5.154) 

where we come back to the notations of k, and fi'^{ik) denotes now the first 
several terms of the uniform asymptotics of the Jost function in the limit ()5.148p 
which cause the divergence in ()5.153|) . namely the terms up to power. Thus 
we arrive at the statement that the value of the expression 

^/ = ^r-eg _ ^asym (5.155) 

is finite at the limit s — ^ 0. The analytical representation of fi^{ik) has been 
obtained recently in Sec ]5.31 

Since the hifi'^{ik) is defined we make a factitious manipulation with E 
fj5.136p namely, we add the identical zero expression E^^^ — S^^"^ to the E' 



Ten 
ren 



and denote the term in the second bracket as the "asymptotic" part of energy 

Thus we have finally the renormalized energy E^'^"' consistent of two parts - 
the finite and the asymptotic one: 

^^en^^/_^a. (5.158) 

Note, that the both parts are finite. Consider the calculation of each of these 
parts in details. 
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5.4.3 The finite part of the ground state energy 



It is supposed that the finite energy E-^ , defined (for scalar field) by ()5.155jl as 

oo 

= E J - m') — [\nm) - In = 

^ CO 

Cs V / dkik'^ - my~'— In fr\ik) (5.159) 



remains finite in the limit s — > since ah uhraviolet divergences are removed. 
From the explicit construction of \n.fi^^{ik) the asymptotic behaviour of this 
function at / — >• cx) is expected to be at least 

\nfr\ik) ~ 0(i) (5.160) 

In particular, for the spinor field problem the uniform asymptotics of the Jost 
function does not contain in fact any even powers of v because of the symmetry 
properties of this scattering problem (see Sec J5.2j) . Thus the asymptotic be- 
haviour of Infi^^iik) is really of order 1//^, that is also confirmed by numerical 
evaluations. 

The correct subtraction provides indeed a well convergent integral over k for 
each / with following convergence of the sum over z/(or /) 

A couple of numerical examples fFig JlSllBj) in SecEl illustrates how it works. 

1. Ef for the scalar field in rectangular scalar background 

We start with the definition of E^ ()5.155|) given explicit in ()5.159|) 
Using the symmetry of Jost function with respect to sign of I .Sec I5.2tj5.49p . 
we can rewrite the E^^^ 

oo oo °° 

m m 

= eI + E(, (5.161) 

oo 

Ei = Cs^ Jie-m'y-^-^[\nfo{tk)-\nfS%tk)], (5.162) 

m 

oo 

E[ = 2C,-^ E y - m'Y~' — [\nfi{zk) - In frm]. (5.163) 
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These expressions are suitable to calculate immediately. To this end we need 
the uniform asymptotic expansions ln/"'*(iA;) of \nfi{ik), I > and for ln/o(i/c) 
respectively. The asymptotic expansion for ln/"*(zA;) has been obtained from 
the iterations of Lippmann-Schwinger equation , Sec l5.H[ For lnfo{ik) it can be 
done immediately through the substitution of the uniform asymptotics of Bessel 
functions K, I in the exact Jost function ()5.48|) . as given in Math.App ()8.10|) . 

2. for spinor field in magnetic flux n 



The sum over orbital momenta / in is performed for the spinor problem 

by means of the substitution / u, SecE] ()3.31|) . We recall the representation 
SecO dOl: 

oo 

Eo{s) = Cs Yl I dk{e -m^f-'du[\nf^{ik) +\Yif-{ik)l (5.164) 

1^=1/2,3/2,...;; 

where 



2^ 1 r(s — 1) COSTTS 



Cs = -f^^'^T^. ^— — , (5.165) 

and define with the help of ()5.159p the E^ as: 

oo 

Ef = Cs Yl [ dk{k^ -m^Y-'dk[\nf+ ilk) + In f-{zk)-\nf::'+{zk)~ In f::'~{zk)] = 

oo 

Cs [ dk{e -m^Y~'dk\nf^''\ik). (5.166) 

i.=l/2,3/2,...^;; 

It is important to notice, that the Jost function In fi,{ik) in contrary to the 
case In fi{ik) of the scalar problem, possesses no symmetry properties like fl5.49|) . 
so the exact Jost functions ()5.54l I5.55p . Sec 15. 21 for spinor field are different. 
However it has been shown in [20 by explicit calculations, that the uniform 
asymptotics for (j5.541 15.55j) coincides with each other, so that 

ln/;^+(2A;) = ln/;^"(zA;) = ln/;^(2A;), A;, z/ ^ oo, k/u = const (5.167) 

and thus we can define 

In f^^\tk) = \nf+iik)+\nf;itk) -2\n f^%tk). (5.168) 

The numerical evaluations SecO confirms (|5.167|) and provides finite results 
for Ef. 
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5.4.4 The asymptotic part of the ground state energy 

Consider the asymptotic part E""^ of the renormalized energy, defined by ()5.157|1 . 
The divergent part E"^^"" has been defined from the heat kernel expansion and 
will be not further simplified. The S^^"^ allows however advanced transforma- 
tions. We start with the definition of E^"^^ ()5.157|) and perform the analytical 
continuation of S^v^ as a complex function of s in the limit s ^ 0. So we have 

^asym = J ^Hk' " m^) — lu W . (5.169) 



First of all we suppose that initial value of s is large enough to provide con- 
vergence of the integral over k and of the sum over / as well, so that the sum 
and the integral can be interchanged. In the section 15.31 it has been shown that 
the uniform asymptotics of the logarithm of Jost function for all the problems 
considered here can be represented in the form: 



oo 



The functions Xnj = Xnj{r) entering here are calculated in the section f I5.3|) . 

1. Eas for the scalar field in rectangular scalar background 

Let the asymptotic part of the energy be written in the form 

where the contribution of the / = in the sum over / is separated, 

oo 

I (5.172) 

m 

and 

oo °° 

g- = - ^^^^r"(.-tr' ^ / dk{e-m'y-^^\nfrm. (5.173) 

Here the symmetry property ()5.49| established in Sec l5.2|l of the Jost function 
In ff'^{ik) is used. The further calculations are given in details in Math.App l8.3.2l 
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Thus the finite part of E"''^ is represented in terms of integrals not complete 
calculable analytically, but well convergent numerically. 



2. E'^^ for spinor field in square shaped magnetic flux Q 

The procedure is in principle the same as in the case of the scalar field. 
Starting with the explicit formula for the S° 



7 as 

•'oo 



oo 3 



E- = 2C. I dki^-m^Y-^^ I '^Y.J: I Xn,'--E-^ (5.174) 

with the constant Cg = 2^[1 + s( — 1 + 2/n2/i)], what is indeed the constant Cs 
introduced by (j5.165p . expanded at s ^ 0. The exphcit substitution of Xnj 
(|5.114j) in (|5.174j) provides under application of the Abel-Plana formula, but in 
this case it must be summarise over half integer values of u (Math.App. I8.50|) . 
and the sum is represented by integrals as 

E"' = E^' + E^' , (5.175) 

where the resulting expression is separated again into finite E'"* and singular 
(divergent) parts E'f*. The derivation of E^^, i?2* coincides formally with the one 
for the case of spinor field in Higgs background, given in Math.App. (|8.1.5p (up 
to the coefficients X^j, which are given by (|5.114p ) 

The finite part of -E"** ( i^f*) is performed in the similar way as described 
recently for the scalar problem. 

After the partial integrations over u several times, that is necessary to get rid 
of divergent powers, the finite part i?2* results again in the expression, represented 
in terms of integrals not complete calculable analytically, but well convergent 
numerically: for spinor field in square shaped magnetic fiux Q it reads 



oo 

—4 f dr 
Tt J 



= fi2ei(i?i, R2) + fi'e2(/2i, R2) (5.176) 

the gi{x) are given in Math.App. ()8.63|) . The coefficients 61,62, entering here 
are to evaluate numerically, the corresponding results are represented in SecEl 
(FigHHD 
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3. E'^^ for abelian Higgs model with NO-vortex 

The calculation procedure for E'^^ is sufficiently the same as in the case of 
spinor field in magnetic background considered above. Since there is no explicit 
analytical form of potentials /(r), v{r), the integration over r is included in the 
double integral being executed numerically using the same package " NIntegrate" 
of Mathematica. 

The analytical formulae for numerical evaluations is derived and represented 
in Math.App. 18.731 sufficiently in the similar way as it has been done for the 
previous case. 
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5.4.5 Cancellation of divergent contributions in the asymptotic 
energy 

As it was described above, the asymptotical part of the renormalized ground 
state energy is constructed as a difference of two parts i?*^™ and S^, each of them 
is infinite in the limit s ^ 0. On the other hand the complete expression of the 
renormalized energy ()5.136|) was presupposed in the top of Sec J5.4l to be finite. 
Since the part defined above ()5.155|) remains finite due to the cancellation of 
divergences between E"""^^"^ and E^'^^ , (what has been also confirmed numerically), 
the same must take place between divergent contributions in E'^'^^ and E""^"^"^ in 

Fortunately one can see how does it come about analytically, that is a non- 
trivial fact. It could seem to be surprising that the divergent part of the asymp- 
totic expansion of energy S^^ coincides exactly with the divergent energy E'^^^ 
defined from the heat kernel expansion, which are originally different objects. 

The minimal necessary condition, which is to expect is the coincidence of 
terms with negative powers of s in and E'^'^^ in the limit s ^ (speaking 
exactly, there is only the first order pole 1/s). It assures that our renormalized 
result is finite. 

This feature can be also used for explicit calculation of the heat kernel co- 
efficients a„ if the uniform asymptotics of the corresponding Jost function is 
known. Such an example is demonstrated in Sec. 15.61 The E'^^^ ()5.137p contains 
not only pole contributions, which are definitive to be cancelled, but also finite 
terms, which appears also in and provide partial or full subtraction with the 
corresponding ones in E'^^'". 

We consider the mechanism of this subtraction on the example of the spinor 
field in the magnetic background. For this case we have the single term 
contributing to the divergent energy: 

where 02 is the only non-zero heat kernel coefficient a„, < n < 2, contributing 
to the E"^^^ (15.13711 . We recall the aa, obtained in Sec. IOlEa. (l5.14;-i|) : 



00 

a2 = ^n'^ j r drh{rf, (5.178) 



VL and h{r) are the same as introduced in Sec. El Eq. ()3.30[ I3.22|) . Now we 
construct the E^^ in 8^ ()8.(ifij) using ()8.7()|1 as 



°° , 3 3n 

dr 







r 

n=l j=l 
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Here for the problem of pure magnetic string the corresponding coefficients Xnj 
from Sec. 15.31 Eqs. ()5.114j) are to use. By performing of the power expansion at 
the hmit s and collecting of the powers and s° we obtain: 



9^ 
127r 



1 , V 
- + ln-^-2 



dr 



— {a{r)a'{r)r — a{r)a{r)). 



(5.180) 



Then, taking into account, that 



aa 
r 



raa 



aa ^ (a 



i\2 



we rewrite the integral in the (j5.180p as 



dr 



a{r)a"{r)r — a{r)a'{r)) = a{r)h{r] 



r drh{r) 



(5.181) 



From the explicit shape of the potential ()3.211 13. 22^ . SeclHl it follows, that the 
surface term is zero and the £^'^ is equal exactly to E'^^^ ()5.177p . thus the both 
terms in E'^^ cancels each other. 

In particular it confirms that the subtraction was executed properly and the 
terms containing pole at s — >■ cancel each other. 

We could also consider the contribution in ()5.179|) of the coefficients Xnj 
with even n (X22, ^24, -^26)? which are also present in ()5.114|) . The immediate 
calculations provide the term 



81' 



-m- 



fi2 

16 



r dr 



air I 



2a'(r)r 



fi2 a^(r) 



-m- 



16 



(5.182) 



The power of mass suggests the correspondence to the heat kernel coefficient 
03/2. It shows in particular that if the potential a(r) decreases at r — > as power 
^1/2 slower, the non zero contribution is provided. For our potential 1)3.211] 
this is not the case. Generally speaking, we can derive two conditions for the 
potential a(r) at r = to keep the 8^%^^^ be zero: 



a(r)| 



r=0 



d_ 

dr 



a{r) 







(5.183) 



r=0 



The first requirement is clear since a(r) corresponds to the tangential component 
of the vector field A, thus the condition a(0) 7^ would mean for this vec- 
tor potential to have a singularity at r = 0. If we have ■§pa{r)\^^_Q 7^ 0, then 
a(r) being viewed as a surface in cylindric coordinates {r,(f,a{r)} has a conical 
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non-smoothness at r = 0. In both cases the coefficient 03/2 in the heat kernel 
expansion is different from zero. Thus we have an imphcit suggestion, that even 
powers of the orbital momenta does not contribute to the ground state energy 
because of 03/2 = 0. 

It could be suitable to remark, that the exact cancellation takes place only for 
special choice of the uniform asymptotic expansion of the Jost function , namely 
here the choice of the half-integer orbital momenta v is crucial. If we would have 
considered the Jost function depending on the original integer momenta Z, the 

would not coincide with the E'^™ in finite part. 

Further we restrict us on short remarks to the two other models. 

For the spinor field in the electroweak vortex background we proceed 
in the similar way as it has been carried out recently, but in the expression 
of E^^ ()5.179p we have to use the coefficients X„j from ()5.135|) . It brings the 
according to (j8.67t I8.68|) to be exactly the same expression as E'^'^'" defined 
from (j5.177p with the non zero heat kernel coefficients 0.1,02 

(I5.14ei|l . thus the 

cancellation comes about completely. 

In the case of the scalar field in the scalar background, the divergent 
part of E'^ is the E'^^ ()8.56|) . which subtract out exactly from the i?"^™ ()5.140|) . 
The three other constituent parts f^g**, £^2^*5 ^3* pure finite and contribute to 
the finite part of i?"*. 
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5.5 Renormalization and interpretation 

We defined in Sec. ()5.4|l the renormalized ground state energy as 

E'^"" = Eo - E^'\ (5.184) 

where E*^"^"" has been obtained from the heat kernel expansion Sec. ()4.H| I4.5|l . 

The main point of the renormalization in quantum field theory is the redef- 
inition of the theory constants through the appearing counter terms. Certainly, 
here we take in mind renormalizable models. Generally speaking, if one obtains a 
finite effective action, one needs to give an interpretation to appearing (possible 
new) effective constants. 

The physically consequent interpretation of subtraction of the terms contained 
in E'^'^'" is sometimes not quite well defined and not unique. 

For the renormalized ground state energy the asymptotic dependence on pow- 
ers of m at m — >■ oo follows from Eq. ()4.119|) to have the form 

Er^rZooY.^ (5.185) 

n>2 

with some coefficients Cn- 

In accordance with the interpretation of E'^^'^ it must vanish in the limit of 
m — * oo since it is the energy of vacuum fluctuations so that E''^'^ fulfils the 
normalization condition [T7] 

limE'''" = at m^oo. (5.186) 

Through the subtraction of terms containing all non-negative powers of m? 
(which are the terms of heat kernel expansion up to 02) the condition (j5.186|) is 
satisfied automatically. 

Some comments on the subtraction scheme are in order. The scheme had 
been used in a number of Casimir energy calculations. In [TH] it had been shown 
to be equivalent to the so called "no tadpole" normalization condition which is 
common in field theory. It should be noticed that this scheme does not apply to 
massless fluctuating fields, for a discussion of this point see [22] • 

We consider below some further examples of renormalization caused by this 
scheme. 

5.5.1 The counter terms 

Proceeding from the HKE we defined the terms with non negative powers of m 
(all the contributions up to 02) as i?*^™ and subtract it out. Thus the renormalized 
energy is defined by (j5.184p and is now finite. The subtraction of E'^^'" implies the 
addition of corresponding counter terms to the initial lagrangian. Since the whole 
energy of the field system consists of the ground state energy of quantum field 
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and classical background energy, the form of counter terms should be suitable to 
redefine certain parameter of the model, that enters in the classical energy. 



1. The counter terms in the scalar field model with scalar step- 
function background 

The counter term contains the contributions of the heat kernel coefficients ai 
and 02 only since all other are zero for this background fSec l5.4l4.3p . 

The coefficients ai and 02 were calculated in Sec 14. 31 ,Eq. ()4.52f4.54j) in terms 
of the scalar potential V(x). To make the renormalization descriptive we remem- 
ber that V{x) was introduced in SecOlas 

V{x) = a^xy (5.188) 

to be quadratic in the background the interaction constant is thereby di- 

mensionless. 

Now on comparing the classical energy 

= i{ao$<9°<l> + (V$)^ + M2$2 + A$^} (5.189) 
with the resulting counter term 

= -^^(- + In ^ - 1) - ^(- + In ^ - 2) (5.190) 

one can ascertain that the parameter A and M that are the self-interaction cou- 
pling and the mass of $ respectively, become redefined in accordance with: 



r , «^ .1 , 4u^ 
^ ^ ^ = ^ + TT;r^ -+ln-^-2 5.191 



M2 ^ M^ = M^ + ^{- + \n^-l) (5.192) 
so that both parameters has been subjugated the infinite renormalization . 



2. The counter term for Dirac field interacting to the magnetic 
background A^{x) 

The single counter term reads: 



= -^(- - 2 + In ^)o2, (5.193) 
327r"' s 
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where 02 for this problem ()4.71|) 

a2 = —n^ f r dr h{rf = e^— f r dr h{rf (5.194) 
3 y Stt J 

On the other hand we have for the classical energy of the magnetic background 

mi: 

^class ^ iFf^'F^, ^^""^^ " ^^^^^ 

Thus the addition of ()5.193|1 to the ()5.195|1 implies the redefinition of the 
coefficient in front of the Fermi term as follows 



,2 /I , All 

. J- - -2 + ln — 
4 4 3 V s m 



+ e2-(--2 + ln-5^ , (5.196) 



that corresponds to redefinition of the cosmological constant by the factor and 
to the renormalization of the electric charge [H] 
e — s> e', so that 

118/1 , V \ . ^ 

The electric charge becomes in such a way finitely renormalized. 



3. The Dirac field with Yukawa-type coupling to the Nielsen - Ole- 
sen string 

Consider the counter term provided by the divergent energy i?*^™ is represented 
in terms of NO-vortex solutions /(r),t>(r). Sec ]5.4l ()5.148|) . 

On the other hand we have the classical energy of the NO background (recall 
the results of SecO]): 

£^'-^ = 7r j rdri^-^^v\rf + r^'[f\rf + l^{l-v{r)f] + ^{f{rf^ 

(5.198) 

On comparing these two expressions we have to establish that the counter term 
fl5.148p does not allows to carry out the complete renormalization which could be 
consequently interpreted as in two previous cases. Indeed, the term proportional 
to f'(r)^ provides the standard renormalization of the electric charge q in the 
classical action. The divergences proportional to f'{r)^ can be absorbed through 
the infinite renormalization of the scalar coupling constant A. Then the remaining 
freedom is the change of the condensate rj, that cannot be redefined properly 
because the structure of terms proportional to rj in E'^™ distinguishes obviously 



99 



from the one in the E'^^"''^'^. This is not surprising since the model given by ()3.36p 
in SecO] contains a non-polynomial interaction and therefore is algebraically non 
renormalizable. There is however some possibilities to give an interpretation of 
this renormalization procedure by considering the classical background fields to 
satisfy the equations of motions. An example of such on-shell renormalization is 
given in Math.App. 18.41 
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5.6 Higher orders of the uniform asymptotic expansion of 
the Jost function 

and the heat kernel coefficient 05/2 

The magnetic background considered in Sec© has singular surfaces where the 
magnetic field jumps. Here we represent the calculations have been done in |33j . 
The heat kernel expansion for the case of singularities concentrated at surfaces 
has been considered in [j2ll, |Elj; jlEj; |1ZI]- Although the general analysis of 
is valid for our background, an explicit expression for 05/2 has not been calculated 
before. 

We can use our obtained Jost function ()5.54|) to calculate the coefficient 05/2 
in the heat kernel expansion. 

Suppose we have obtained the value of Eq{s) ()5.43|) in the point s = —1/2. It 
follows from ()4.9()j) that 

^ n=0,i,l,... 

and at the hmit of s — | only the term containing 05/2 remains to be nonzero 
in the sum. 

^o(^) = -^(£^«5/2 (5.200) 

From the other hand we have ()5.43|) 

00 

Eois) = Cs Yl I - ^')'~'^/-(^^) = ^Ms). (5.201) 

i/=l/2,3/2,...;; 

Here we substitute the exact Jost function by its uniform asymptotic represented 
in the form 

fua^^f^^^ ^ hn{t^^ ^5202) 
n=l,3,4,5... ^ 

where the coefficients hn are functions of ti = (1 + (^^)^)~^^^, 

^2 = (1 + the power is absent, as noticed above (Sec l5.3p . 

r- /^^^r(.-l) -4sin(7r.) 

- ~40Fr(.-l/2) n ^^-^^^^ 



At the limit s ^ — | it yields 



Eoi-l) = i^^^'s^-i^h{s) (5.204) 
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and therefore we obtain for 05/2 



^5/2 



Res^ 



Ms) 



(5.205) 



To obtain the exphcit form of ReSg^_ih{s) we use the uniform asymptotic 
expansion of the Jost function (j5.54p . The terms /in(^i) ^2) in uniform asymptotic 
expansion can be obtained either by iterations of Lippmann-Schwinger equation 
(see 1201) '^^ t)y using the exphcit form of the Jost function as weU. Ah 
the further terms up from are produced from the exphcit form of the Jost 
function ()5.54j) because of comphcation of the first way for higher orders n (see 
the remark to 15.1141 in Sec ]5.3|) . Namely, we obtain several higher orders l/u of 
uniform asymptotic expansion for special functions I^,, Ky,i Fi which the exact 
Jost function ()5.54l I5.55|) consists of (it can be done starting with the explicit 
form for two first orders and executing the recursive algorithm several times PP), 
then after substitution of each of functions /, K,i Fi by its corresponding uniform 
asymptotic expansion and separation of powers of v we arrive at the form (|8.13|) . 
The coefficients hn(ti,t2),n = 1, ...,4 are given in the Appendix (|8.13|) . 

If the function /in(^i,^2) is a polynomial over ti,t2, so we can consider some 



1,2) of it. Notice, that for /ii(ti, ^2), ^3(^1, ^2) it is not the 

fc-integer, as an infinite 



i+t 



and 



term P ,{t = U, i 

case, but we can treat the terms of kind 
sum of powers t. 

Apart from the construction of t = ti,t2 (!5.202|) they are strictly positive 

00 00 

and less than 1, therefore each of the series in = Yl^'^^ ~ t'^^'^^ converges 

i=0 i=0 

regular and uniform and respecting that we have only finite integer powers k so 
that j P' ] converges as well (the k',i' are integer). The sum over i' can 

\i'=0 J 

be interchanged with the one over u (Princeheim's Theorem) thus the following 
procedure for some integer power P is valid. 
Performing the sum 



h{s) 



[dk{e-m'y-^^ KARuR2)^ (5.206) 



!^=l/2,3/2,...,;^ n=l,3,4,5 

by meaning of ()8.50p we obtain the sum of two parts 



00 00 



h{s) 



jduj dA;(A;2-m2)i-^^^/i„,,(i?i,i?2)^+ (5.207) 



m 

00 



du 1 



00 



p 
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where t = ti^2- The first summand of ()5.208|) gives for each power j oit = {ti, ^2} 
(using the 

m3 r(2-s)r(i±p)r(. + n^) 

2 (i?m)-ir(j72) ^^-^"^^ 

(where R denotes Ri and -R2 respectively) at the hmit s — — | only the terms 
corresponding to n = 2 and n = 4 (which will be calculated below) have a pole. 
For the second part of the Abel-Plana formula ()5.208p we have 

00 

1 [ du r(2-.)r(s + j72-i) 

I J 1 + e2-- ji?4-2«r(j72) ^ ■ ' 



[(zz/)2+^-"[(mi?)2 + (zz/)Y~''/'^' - (-«z/)'+^""[(mi?)2 + (-^zy)2]i-^'/''^] 



X 



It can be seen using the Laurent expansions at s — | of Gamma functions 
entering the (j5.21(Jj) . that only the terms containing -1, 1, and 3 powers of ti can 
contribute to the residuum at s — * — |. But the expression in square brackets in 
()5.210j) can be performed as 

-z/2(zy2 _ (mr)2)[(zz/)J-"e'"(^-^'/2-^) - (_izy)i— e-'"^^-^'/^-^)] (5.210) 

thus for j ~n even the expression in square brackets produces (dropping the non 
sufficient coefficient ±2 or ±2z): 

for j = —1, 1 : sin7r(l — j/2 — s) = costts 

for j = 3 : sin7r(l — j'/2 — s) = — cosTTS 

and for j — ri odd respectively 

for j = ±1 : cos 7r(l — j/2 — s) = ± sin its 
for j = 3 : cos7r(l — j/2 — s) = — siuTTS 

But for s ^ these functions behave as 



2 



siUTTS ~ -1 + ^(i + s) 



2 



COSTTS ~ 7r(| + s) 

and it means that only the contribution of terms with odd powers j—noi iv could 
survive and these corresponds for the possible values of j to z/"^, .... In 

fact the coefficient /i2(s) at i/"^ is zero, and all other possible terms up from 
does not contain any powers of ti lower than 4; one can see it for example in the 
explicit form of the uniform asymptotic expansion of special functions entering 
the ln/^(ifc) ()8.13l Sec l8.1.4|) . Therefore the second summand of ()5.208|) does 
not produce any contribution to the Res^^_ih{s). Thus we have that only the 
contribution from the first summand of ()5.208|) remains, and since the term of 
n = 2 is zero the searched residuum resulting from the term of n = 4 is: 
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where the fj5.208j) and the exphcit form of hi{ti,t2) = ^j^PLj^2)2 [Riiif — ^i) + 
-^2(^2 ~ ^2)] have been used. Finally we have 

and therefore 

5n^/Ti\Ri + R2) 
"^/^ = {Rl - Rlf ^^-'^^^ 

This is the heat kernel coefficient 05/2 for the configuration of the magnetic back- 
ground field as given by Eqs. (j3. 20113. 22p . 

We can calculate the heat kernel coefficient 05/2 in a more general situation 
when the magnetic field jumps on an arbitrary surface S. The coefficients a„ for 
n = 1/2, 2 can be read off rather general expressions of the paper jlHl- Let 
be values of the magnetic field on two sides of S. According the analysis of |46j . 
the coefficient 05/2 must be an integral over S of a local invariant of canonical 
mass dimension 4, which is symmetric under the exchange of B~^ and B~ and 
which vanishes if B~^ = B~ (i.e. when the singularity disappears). There is only 
one such invariant which gives rise to the following expression: 

05/2 = e j{B+ -B-fdii{T.), (5.214) 
s 

where the integration goes over the surface S and -B^ are the values of the 
magnetic field on both sides of S in the given point. The yet undefined constant 
^ can be found using Eq. ()5.213p which constitutes a special case of ()5.214|) . Here 
the surface E consists of two circles in the (X, y)-plane so that 

j {B+ - B-ydfi{J:) = 2n{Ri + R2)B^, (5.215) 
s 

where the jump (5+ — B^) is just the value of -B at r G [-Ri,-R2]- With this, 
Eq. ()5.214|) takes the form 

05/2 = 27re^^B\Ri + R2). (5.216) 
On the other hand, from Eq. ()3.221 13.30p we have 

$ = 5 J Bd^x = TxB{Rl - Rl) (5.217) 

r(i[Ri,R2\ 

and with Ea. (j!On|l from Ea. (|5.213|) it follows that 



2^2 -^1 + -^2 
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comparing ()5.216|) with ()5.218p we get 

« = T^r. (5.219) 
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6 Graphics and numerical results 



6.1 Scalar string 

6.1.1 Finite and asymptotic vacuum energy 
1. Finite energy 

The finite ground state energy for scalar field is calculated numerically by 
using of ()5.161|) . The plot of the In fi^'^{ik) entering here is shown on Fig.® 
for several values of momentum k. It is multiplied by to make the asymptotic 
behaviour at / — oo prescribed in Sec l5.4l visible. The FigO shows some plots of 
the integrand over k in ()5.173|) for several orbital momenta /. It can be seen that 
the shapes of curves provide well convergent integrals for each of /. 




50 100 150 200 250 



Figure 2: ln/f"''(iA;) multiplied by at Ri = 1, i?2 = 2, Vq = 1.1 for A; = 1, 5, 10, 20 




50 100 150 200 250 



Figure 3: klnf[''\ik) at Ri = l,i?2 = 2,Vo = 1.1 for I = 30,50,80,100 
2. Asymptotic energy 
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The asymptotic energy E"'^ contains contributions of orders Vq and Vq of the 
potential respectively. We represent the E""^ to be consistent of these two terms 



= ei(i?i, R2)Vo + e2{Ri, R2)V, 



; 



(6.1) 



and the coefficients ei(i?i, i?2), e2(-Ri, -R2) are shown in FigEl - FiglH for fixed 
value of potential Vq = 33.1. On the first plot the outer radius is fixed and the 
inner one moves from inside to them - the asymptotic energy shows physically 
"well" behaviour, both contributions vanish at infinitely thin potential wall. The 
second plot demonstates the growth of both contributions if the inner radius is 
fixed and the outer increases. 




Figure 4: ei,e2 at Vb = 33.1; Ri = 0.01, 0.010001 < R2 < 0.9999 



e, (R, R, ) 




Figure 5: ei,e2 at Vb = 33.1, R2 = 1., 0.010001 < Ri < 0.9999 



6.1.2 Parametrical dependence of vacuum energy 

We are interesting for the dependence of the ground state energy on the external 
parameter, such as the inner/outer radia -Ri, R2 of the background and the height 
Vo of potential as well. 
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The remarkable properties occur, if the height of scalar potential Vq becomes 
small. Generally, the ground state energy remains positive almost for all values 
of Vq as shown in FiglTOl However, it exists a very thin region of potential 
values between Vq ~ 0.01 and Vq ~ 0.065 where the E^^"" is negative (FigEj) and 
possesses generally a global minimum (FiglTT|l 

While the ground state energy changes the sign from negative to positive one, 
it shows the non-trivial behaviour which is represented by Fig|7| The critical 
curve has two extrema: at i?2 ~ 1.067 (minimum) and i?2 ~ 1-05 (maximum). 
At ~ 0.062 the critical energy approaches the deepest minimum (FiglH}, and 
at V ~ 0.071 the minimum vanishes (FiglH)). 

Since the dependence of classical energy of the background is pure quadratic 
on i?2, this behaviour is kept also for total energy, the positions of extrema will 
be shifted. 

Now consider the dependence of E^^"" on the potential height Vq if Ri,R2 
are fixed. The Fig shows the slight minimum at Ri = 1,R2 = 11.5, that 
corresponds approximatively the global minimum of E^*^" at Vq = 0.003. If the 
potential Vq grows, the dependence becomes obviously quadratic in the leading 
order, (see Table I and Fig |13ll4|) . Since the classical energy is linear on poten- 
tial, this behaviour is kept also for total energy as well. 



Vo 


E^ 






0.001 


3.7765 10-6 


-4.41055 10-^ 


-6.3405 10-^ 


0.101 


-0.00022342 


0.000609497 


0.000386072 


0.201 


-0.00098331 


0.00331244 


0.00232913 


0.301 


-0.0017050 


0.00810441 


0.00639944 


0.401 


-0.0018532 


0.0149854 


0.0131322 


0.501 


-0.0009434 


0.0239555 


0.0230121 


0.601 


0.0014908 


0.0350145 


0.0365053 


0.701 


0.0058864 


0.0481626 


0.0540491 


0.801 


0.0126312 


0.0633998 


0.076031 
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Figure 8: Both extrema of E^^'^ vanish at Vq = 0.071; Ri 




Figure 9: The E'''^" is negative for potential Vq between ~ 0.01 and ~ 0.065 
6.2 Magnetic string 

6.2.1 Finite and asymptotic vacuum energy 
1. Finite energy 

The finite part of the ground state energy, Ef, is used in the form as given 
by Eq. ()5.166|) . Let Ef{u,k) denote the function to be integrated and summed 
over in that expression. In FigEI it is shown as function of k for several values 
of the orbital momentum v. In order the make the behaviour better visible it is 
multiplied there by v'^k^ . These functions are smooth for all values of /c, starting 
from some finite values at A; = 0. 



For large fc, the functions v'^k'^Ef{v,k) shown here tend to a constant thus 
the integral over k is convergent. All integrals have been truncated at = 1500. 
The error caused by this is quite small and does not change the results shown 
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1.2 1.4 1.6 1.8 



Figure 10: Samples of behaviour of E'^'^"' for 0.1 < Vq < 1.1 



20 40 60 80 100 




Figure 11: For small Vq = 0.003 i?^^" demonstrates a global minimum. 

in the Table II. These integrals we denote by Ef{v). They are shown as func- 
tion of V in Fig^l in a logarithmic scale. Again, we multiplied by a power of 
argument, here by z/^, in order to make the behaviour for large u visible. It is 
seen that v'^Ef{v) tends to a constant so that the sum over v is convergent. The 
sum is taken up to = 232.5 and again the remainder is small. The calcula- 
tions have been performed for several values of the parameters. The results are 
displayed in Table. II. The computations are performed with an adapted arith- 
metical precision. In intermediate steps compensations between sometimes very 
large quantities appeared. The precision was adapted accordingly. For example, 
for Ri = 0.99, u = 250.5, k = 1600 as much as 1404 decimal positions have been 
necessary to get at least 16 digits precision of the integrand Ef{v^ k). This was a 
factor causing large computation time. 
2. Asymptotic energy 

The asymptotic part of the ground state energy as given by Eqs. (|5.176|) to- 
gether with (|8.63p can be represented as sum of a part proportional to the second 
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Figure 12: For fixed Ri = 1,R2 = 11.5, _E'^"-dependence on Vq has a sHght 
minimum. 
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Figure 13: , E"', E^^" of Vq for Ri = l,i?2 = 11-5 



power of the flux and one proportional to the fourth power, 



E"' = — [ei(/?i, R2)n^ + e2(i?i, R2)n% 

TX 



(6.2) 



The corresponding coefficients ei{Ri,R2) and e2{Ri, R2) can be calculated nu- 
merically without problems. They are shown as functions of R2 for flxed -Ri in 
FigUi (muhiplied by R^). 

6.2.2 Parametrical dependence of vacuum energy 

Table II The numerical evaluations for several values of Ri and Q 



i2i=0.0001 
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Figure 14: E™" divided by V^. 
-Jk lnff(ik) 



. 002 




Figure 15: for Ri = 0.0001, R2 = 1., Q = 3. the bundle of -ifc In /^^"''(zA;); i < 
< 250 i, multiplied with k'^ 



n 








j^class 


J^tot 


0.5 


-0.0130152 


0.00152886 


-0.0114863 


3.14159 


3.13011 


1. 


-0.052403 


0.00509591 


-0.0473071 


12.5664 


12.5191 


3. 


-0.450093 


-0.052012 


-0.502105 


113.097 


112.595 


6. 


-0.91953 


-1.52936 


-2.44889 


452.389 


449.94 


10. 


4.85954 


-12.9483 


-8.08871 


1256.64 


1248.55 


15. 


46.3333 


-67.3661 


-21.0328 


2827.43 


2806.4 


21. 


215.273 


-261.526 


-46.2534 


5541.77 


5495.52 


30. 


987.62 


-1095.29 


-107.666 


11309.7 


11202.1 



i?i=0.9 
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Figure 16: Contribution of each integral ~ J kin f^'^^{ik)dk (multiplied with v'^ 
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Figure 17: The contribution Eflu, k) multiplied by u'^k^ of the individual radial 
momenta for several orbital momenta to the finite part of the ground state energy 
for i?i = 0.0001, R2 = l,n = 3. 
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j^class 


J^tot 


1. 


-311.182 


-0.425555 


-311.608 


6286.33 


5974.72 



The results are displayed mainly in Table II. For small inner radius of the flux 
the results are close to them of jTHj where the same problem for a flux tube with 
homogeneous magnetic field inside, which corresponds to i?i = here, was con- 
sidered. Especially, it is seen that for large flux Q there is a compensation of the 
(^^-contribution between the finite and the asymptotic parts of the ground state 
energy leaving a behaviour proportional to fl^lnfl as shown in Fig|2ni Here the 
asymptotic part gave an essential contribution. The ground state energy remains 
negative, but numerically small. Only for very large flux it could overturn the 
corresponding classical energy, but these values of the flux are clearly unphysical. 

For values of the inner radius close to the outer one, Ri ^ 1 (where we have 
put i?2 = 1) the picture changes. Here the vacuum energy grows faster than the 
classical one. Generally, both diverge proportional to (1 — Ri)^^, the classical 
energy is equal to ^^ciass = f^^27r(l — Ri)~^. The vacuum energy, multiplied by 
(1 — is shown in Fig J6.3l in a logarithmic scale. It is negative and growing a 
bit faster than the classical one which would be a constant in this plot. Here the 
asymptotic part of the ground state energy becomes increasingly unimportant 
(see Table II). 

6.2.3 Some remarks on the numerical difficulties 

The calculations have been performed for several values of the parameters. The 
results for magnetic string are displayed in Table. 1. The computations are per- 
formed with an adapted arithmetical precision. In intermediate steps compen- 
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Figure 18: The contribution Ef{v) of the individual orbital momenta to the 
finite part of the ground state energy multiplied by u"^ for Ri = 0.0001,i?2 = 1 
and = 3. 




Figure 19: The asymptotic part of energy as function of the outer radius for 
Ri = 0.00001 and Ri ^ 1 resp. 



sations between sometimes very large quantities appeared. The precision was 
adapted accordingly. For example, for Ri = 0.99, z/ = 250.5, k = 1600 as much 
as 1404 decimal positions have been necessary to get at least 16 digits precision 
of the integrand Eflu, k). This was a factor causing large computation time. 
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Figure 20: The ground state energy divided by Vt^hiVL as function of VL. 
6.3 NO vortex 

We investigated numerically the vacuum energy for values of (3 ranging from 
(3 = 0.3 to P = 6 for some choices of the parameters /e and r/. The results are 
represented in Table III. In the Figs. Ifj.Hl and ffi.Hl two examples for the dependence 
of the individual parts of the vacuum energy on the parameter /3 are displayed. 
The constituent parts of the vacuum energy 



Table III 



q 


P 


fe 


V 


-^'class 


gas 




^rcn ^ ~\~ ^ 


0.5 


0.3 


1 


1 


3.902 


-0.006102 


0.02410 


0.01800 


0.5 


0.6 


1 


1 


4.549 


-0.006102 


0.01590 


0.009795 


0.5 


1. 


1 


1 


5.099 


-0.006104 


0.01153 


0.005421 


0.5 


3.5 


1 


1 


6.710 


-0.006146 


0.004224 


-0.001922 


0.5 


6. 


1 


1 


7.507 


-0.006208 


0.001884 


-0.004325 


2. 


0.3 


1 


1 


2.105 


-0.003170 


-0.001883 


-0.005053 


2. 


0.6 


1 


1 


2.399 


-0.002891 


-0.007285 


-0.01018 


2. 


1. 


1 


1 


2.652 


-0.002924 


-0.01182 


-0.01475 


2. 


3.5 


1 


1 


2.652 


-0.002924 


-0.02451 


-0.02743 


2. 


6. 


1 


1 


3.861 


-0.006539 


-0.03022 


-0.03676 


0.5 


0.3 


0.1 


0.1 


0.03902 


-0.00006102 


-9.589 10-6 


-0.00007061 


0.5 


0.6 


0.1 


0.1 


0.04549 


-0.00006102 


-0.00001214 


-0.00007316 


0.5 


1. 


0.1 


0.1 


0.05099 


-0.00006104 


-0.00001421 


-0.00007525 


0.5 


3.5 


0.1 


0.1 


0.06710 


-0.00006146 


-0.00001975 


-0.00008121 


0.5 


6. 


0.1 


0.1 


0.07507 


-0.00006208 


-0.00002216 


-0.00008425 



The vacuum energy i*;'"^" decreases for large values of the vortex parameter 
f3. The contribution of E^'^"' to the total energy of vortex can be positive as well 
as negative, dependent on the magnitude of the electric charge q. Nevertheless it 
remains sufficiently smaller compared to the classical energy Edass- 
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Figure 21: The vacuum energy as a function of (3 for the q = 0.5, fe = 1, rj = 1. 
In order to represent all quantities within one plot the classical energy is divided 
by 200. 
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7 Conclusions and Discussions 



In the preceding sections the ground state energy for a spinor quantum field 
in the background of a rectangular shaped magnetic fiux and of the Nielsen- 
Olesen vortex, and for a scalar field in the rectangular scalar background has 
been calculated. 

The used regularization procedure is a kind of zeta functional regularization . 
The sum over spectrum of hamiltonian is performed in terms of the Jost function 
for the corresponding scattering problem. 

The Jost function for each model has been written down explicitly, also its 
asymptotic part as well. The numerical calculation for several models (spinor 
field) required work with high arithmetic precision. The obtained results pre- 
sume some discussion 



7.1 Scalar field in scalar background. 

This model is the simplest construction which allows to carry out analytic calcu- 
lations as well as numeric evaluations relative easy. 

The regularization procedure with the renormalization of the ground state 
energy and its interpretation is also quite transparent. Nevertheless this model 
looks up a number of interesting results concerning the behaviour of the vacuum 
energy. The vacuum energy can provide positive as well as negative contributions 
to the total energy. The sign of it depends on the magnitude of the background 
potential V. Moreover, there is a thin range of values V, whereat the ground 
state energy shows the critical behaviour for its dependence on the outer radius 
i?2 possessing one maximum and one minimum. It gives reason to suppose an 
existence of stable states of the scalar string, generated by the vacuum corrections. 



7.2 Spinor field in magnetic background. 

The spinor vacuum energy for the constructed magnetic string is investigated 
for several values of potential fl and inner radius Ri. It shows no extrema, so 
the magnetic string constructing in such a way must be unstable. The known 
asymptotic behaviour of the one one-loop vacuum energy for large magnitude of 
background potential in the form In Q is also confirmed fro this model. 

The question whether the vacuum energy for sufficiently small thickness of 
magnetic fiux may become larger than the classical one cannot be answered by 
the numerical results obtained. The problem is that the computations become 
very time consuming because of the increasing precision which is required. Also, 
one has to take higher k and u into account. The weakening of the growth for 
Ri = 0.997 and -Ri = 0.999 (while R2 is fixed to be equal 1) seen in Fig J6.3l mav 
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be caused just by this reason. Here one has to note that the integrand is for 
large k and u always negative (see FiglTTj) so that dropping some part (as we did 
within the numerical procedure) diminishes the result. So, as a result, we cannot 
exclude from the given calculation that the vacuum energy grows for a strong 
background field faster than the classical energy. 



(l-Rj)E 




Figure 23: The vacuum energy multiplied by (1 — in a logarithmic scale for 

n = i. 

Further work is necessary to better understand these questions. An improve- 
ment of the numerical procedure is certainly desirable. It could go along two 
lines. First, in the calculation of the Jost function the compensation of large ex- 
ponentials should be avoided by taking them into account analytically. Second, in 
the compensation between the logarithm of the Jost function and its asymptotic 
expansion in the integrand oi Ef in Eq. (|5.166|) higher orders of the asymptotic 
expansion could be used. However, for this reason one would have to continue 
the procedure invented in |1.9^ for this expansion using the Lippmann-Schwinger 
equation or to find some equivalent procedure. 



7.3 Spinor field in NO vortex background. 

During the analytical regularization of the ground state energy in the model with 
the special interaction term is occurred. Namely the standard zeta functional 
regularization determines the counter terms from the first heat kernel coefficients 
(up to 02). It turned out that to renormalize the vacuum energy one has to 
introduce a counter term which is not present in the initial action and represents 
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a non polynomial interaction. This is caused by the non polynomial interaction 
present in the model itself Sec. ((31). 

In the considered case the background is given purely numerically. This is 
the crucial technical difference to the previous two models where it allowed an 
analytically calculation. However, it has been demonstrated that the similar com- 
putational scheme used here is able to work successfully with such backgrounds. 
This is a step forward to physically really interesting problems. In the considered 
model the stability of the background is given by topological arguments and for a 
realistic choice of the parameters the quantum contribution is small. This small- 
ness has two sources. The one is the smallness of the coupling constants which 
appears in front of the quantum contribution relative to the classical one. The 
second one is a purely numerical smallness of about two orders of magnitude. 
It is present even if the parameters and couplings are all of order one, see for 
instance the Fig l6.3l 

As seen from the table and from the Figs J6. 31 and for the dependence 
of the vacuum energy on the various parameters no general rule can be seen 
so far. Even the sign of the vacuum energy changes. The same applies to the 
relative weight of the individual contributions. Sometimes the asymptotic part is 
dominating, in other cases, however, the 'finite' part is larger. In general, they are 
of the same order. From this one can conclude only that in the given background 
there is no small parameter which could allow for some approximation scheme. 
For instance, if the asymptotic part is dominating one could hope to get a good 
approximation by including higher orders of the uniform asymptotic expansion 
of the Jost function into the asymptotic part of the vacuum energy and neglect 
the 'finite' part of it which is the numerically much harder part of the problem. 
The examples in [13;^ have been a kind of suggesting this way in contrast to the 
example in the present paper. 

For the considered model of a spinor in the background of the NO vortex, the 
vacuum energy due to its smallness has only a very negligible influence on the 
dynamics of the background. Hence in the considered case the vacuum energy is 
of limited physical importance. However, its calculation gave new insights into 
the structure of such calculations and demonstrated the power of the methods 
used. A next step could be to apply them to the Z and electroweak strings where 
the stability is not guaranteed by topological arguments and where the stability 
issue is not finally settled with respect to the fermion contributions [21 EI] • 



7.4 Summary 

Basically, the contribution of the quantum ground state has an influence on the 
macroscopic behaviour of the field system. The degree of this influence depends 
on the nature of the interacting background and fluctuations as well as on the 
kind of interaction and coupling parameters. The vacuum corrections can provide 
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both negative and positive contributions to the total energy. 

It has been shown for scalar and magnetic strings with a simple rectangular 
shape and for the NO vortex. 

In all three cases the general dependence of the vacuum contributions on the 
background properties is not to be formulated certainly. 

For the scalar string, even the alternation of the stabihty properties is possible. 
The large-potential asymptotics known for homogeneous fields is confirmed for 
scalar and spinor strings. 

In the technical aspect, the known methods using the zeta functional regular- 
ization and heat kernel expansion have been applied successfully to more com- 
plicated shapes of backgrounds as before, and in particular to the shape being 
given pure numerically (NO-vortex). 

The obtained results suggest the applicability of the developed methods and 
technics to calculations of the vacuum energy in more complicated background 
configurations, such as monopoles, electroweak strings, skyrmions, sphalerons etc. 

These future investigations would lead to better understanding of the phenom- 
ena of vacuum energy and its relation to the geometry of classical backgrounds. 
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8 Appendix 

8.1 The uniform asymptotic expansions for logarithm of 
Jost function 

8.1.1 The uniform asymptotics of modified Bessel functions 
K and I 

The uniform asymptotic expansion of modified Bessel functions K^^^i (kr), I^^i (kr) 
in the combined limit i^, fc — *• oo, k/u = z = const, can be obtained as series of 
the variable t = 1/ Vl + z'^r'^ P : 

' i = -l ' \ / i = l 



j=— 1 ^ ' 1=1 

i=-l ^ / i=l 

where Sf[a, t] are: 
for Ki,^a, a = ±^ : 

6 1 t, a = m 

"^^ ' ^ t 2 1 +t 

T Int a, 1 + t 

Dn t, a = 1 m 

° ^ ' ^ 2 2 1 -t 

5f = ^(-3 + 12a^ - 12at + 5t^) 

S~ [t, a] = [-8aH + 12a2(-l + 2t^) + a(26t - 30t^) + 3(1 - 6t^ + 5t^)] 
48 

S-[t, a] = {t^ [-375 + 4779t2 - 9945t^ + 5525t^ + 240a^(3t2 - 1)- 
5760 

960a3t(4t2 - 3) - 720at(7 - 22f + 15t^) - 120a2(i3 _ 8At^ + 75t^)] } 
S^[t,a] = 

^— {t^ [-48«^t(5t^ - 3) + 240a^(l - 8t^ + 8t^) - 40a^t(69 - 240^^ + 175t^) + 
1920 

120a2(-7 + 94^2 - 207t^ + 120t^) + a(3561t - 21225t^ - 342255t^ - 16575t^) 
+15(13 - 284^2 + 1062t^ - 1356t^ + 565t^)] } (8.1) 

for I^+a, a = ±^ : 
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S+-. [t al = - + - In 

' ^ t 2 1 + t 

e+r. 1 a 1 + t 

Dn r, a = 111 

° ^ ' ^ 2 2 1 -t 

a] = -^(-3 + 12a2 ^ i2at + 5f) 

S+[t, a] = —t^ [%ah + 12a2(-l + 2t^) + a(-26t + 30t^) + 3(1 - Gt^ + 5t' 
48 

^+[t d = [-375 + 4779^2 - 9945t^ + 5525t^ + 240a^(3t2 - 1) + 

^ ^ 5760 ^ ^ ' 

960a^t(4t2 - 3) + 720at(7 - 22t^ + 15t^) + 120a^(13 - Mt^ + 75t^)] } 
S2[t, a] = 

{t^ U8ah{5t'^ - 3) + 240a^(l - 8t^ + 8t^) + 40a^t(69 - 240t2 + 175t^) + 

1920 

120a2(-7 + 94:f - 207t^ + 120t^) - a(3561t - 21225t^ - 342255t^ - 16575t^) 
+ 15(13 - 284^2 + 1062t^ - 1356t^ + 565t^)] } , (8.2) 

what leads in particular to 



V^v{t) = ] + ln^^ (8.3) 



in dHU EH) 



8.1.2 The saddle point expansion 

The evaluation of integrals of type 



e 





^^("'V(^')f^^', (8-4) 



entering in the expressions of In fu' Sec ]5.3l (j5.102p has been done by using the 
saddle point approach (See. for example |70j). The derivation is described in 
details in [20j. We use the results given in App.[T^ as a starting point: 







oo 

e''^(^')(f){r')dr' = e""^^"^ ^ hk-iv'^ (8.5) 

k=l 



with the first several coefficients of the expansion 

0(r) 0(r)r7"(r) 0'(r) 

ho = — hi - 



ri'{r)' ri'{r)^ ri'{ry 

^ _ r{r) 3<P'{rW{r) 30(r)VXr)^ <P{rW"{r) 
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8.1.3 Scalar field with scalar background 



We proceed from the expansion of the logarithmic Jost function Sec 15. HI ()5.1()2j] 
obtained by iterations of the Lippmann-Schwinger equation 

The terms \n fi^\ik) are obtained as described in Sec l5.3l by substitution of 
CTTTIl in (jnSHI): 



oo 



lnfl '{ik) = J rdrV{r)Ki{kr)Ii{kr) 



oo 

\nfl'^\ik)= J rdrV{r)Ki{kr) r'dr'V{r') [li{kr)Ki{kr') - Ii{kr')Ki{kr)] Ii{kr') 





oo oo 

rdrV{r)Ki{kr)Ii{kr) I r' dr'V{r')Ki{kr')Ii{kr' 



2 





oo 

f(3). 



\nff'{ik)= rdrV{r)Ki{kr) / r' dr'V{r') [li{kr)Ki{kr') - Ii{kr')Ki{kr)\ x 











X I j r"dr"V{r")Ii{kr')Ki{kr")Ii{kr") - j r" dr"V{r")Ki{kr')Ii{kr"f 

vO 

oo oo r' 



i [j rdrV{r)Ki{kr)Ii{kr) j r'dr'V{r')Ki{kr') j r" dr"V{r")[Ii{kr')Ki{kr") 

D 
Il{kr")Ki{kr')]k{kr") + 



vO 

"\ ti '\'\ T n "\ 



rdrV{r)Ki{kr) j V{r')[Ii{kr)Ki{kr') - Ii{kr')Ki{kr)]Ii{kr')x 




oo 

r"dr"V{r")Ki{kr")Ii{kr") I + (8.7) 





oo 

i I rd.rV(r)K,(kr)T,(kr) I r' -^-'^^ ' -'"^ ' ^-'^ ^ ' '-''^ I ^'U^n^r ( n.^n-, r n.^n 



rdrV{r)Ki{kr)Ii{kr) J r'dr'V{r')Ki{kr')Ii{kr') j r" dr"V{r")Ki{kr")Ii{kr") 




The expression In fi^\ik) has already a usable form and In//^ '^\ik) can be 
simplified if we use its symmetry properties in respect of the change r ^ r' . 

The procedure is illustrated by the \n.fl'^\ik). We start with the expression 
(IHSl) 
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In fP{zk) = 

rdrV{r)Ki{kr)Ii{kr) I r' dr'V{r')Ki{kr')Ii{kr') 



rdrV{r)Ki{kr)Ki{kr) / r'dr'V{r')Ii{kr')Ii{kr') 



rdrV{r)Ki{kr)Ii{kr) / r'dr'V{r')Ki{kr')Ii{kr'). 





Then we represent the last term as 



rdrV{r)Ki{kr)Ii(kr) / r'dr'V{r')Ki(kr')Ii{kr') 



rdrV{r)Ki{kr)Ii{kr) 



r oo 



r 



r'dr'V{r')Ki{kr')Ii{kr'), 



and take into account the rule of change the integration order: 

r oo oo 

f{r)dr j g{r')dr' = J g{r)dr J f{r')dr', (8.8) 

Or 

which allows to see, that the last term is equal exactly to the first one but with 
the sign minus in front of them. Thus these terms cancel each other, and we 
have: 

oo 

\nff\ik) = - j rdrV{r)Ki{kr)Ki{kr) j r'dr'V{r')Ii{kr')Ii{kr'). (8.9) 
° 

We do not discuss here the simphfication of In ff\ik)^ because it is not needed 
for minimal subtraction, as pointed out in Sec. 15.31 

Since the contribution of / = in the sum over / in ()5.161l 15.1621 15.163^ 
Sec ]5.41 is separated, the ln/o(ifc) is to perform additionaly. Instead of the uni- 
form asymptotic expansion (because Z = = const) the usual asymptotics at 
/c — > oo is used. 

We use entirely the corresponding asymptotics of modified Bessel functions 



Kniz) ~ e" 



TT 




— 9 
^ 8^ ^ 2(8^)2 ■ ■ ■ 



2z 8z 2(8z) 
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which are to substitute into ()5. 96115. iTM Sec IF!^ with the argument z having the 
values /ci?i, /ci?2, Q'-Ri, (1^2- 

Recall, that q = i/fc^ + Vq, that provides in the considered limit k ^ oo 

Substituted into ()5.48p it provides the leading order for In f^'^^ik): 

\nfS%tk)^-^{R,-R,). (8.12) 

8.1.4 Spinor field with magnetic background 

( To 15.61 Higher orders of the Jost function ) 

The asymptotic of the logarithmic Jost function can be obtained in the form: 
\nU{tk) = E^n(ti,t2),(seeElI2D with: 

n 

hi = ^\^[Ri%i2 + ti)(l + ts)' - SRi^R2\l + ti)%{l + h) + R2\l + + 2t2)] 



h, = -lA2[2i?iV(l+t2)Wi?2'(l + ti)V(l+t2f(3t2'-4)- 
6 

i?2^(l + h)\%3t2^ + Qt2^ - - 8t2 - 2)] + 

^\^[ARi%%A + ti)(l + t2)' - 5i?i^i?2'(l + + t2f + 

lbRi^R2\l + ti)%\l + t2)'(t2' + 2t2 - 1) - 
5i?i'i?2'(l + tl)'t2'(3t2^ + 12t2' + 6t2' - 4^2 - 1) + 
R2\l + tl)'t2'(5t2^ + 20^2' - 4t2 - 1)] 

, _ 5'[R,\\l-t^)+R2H2\l-t2^)] 

4{Ri - Rlf 

^ = {Rl-Rl){l + ti){l + t2) ^^-^^^ 



8.1.5 Spinor field with NO vortex background 

Uniform asymptotics for Jost function in abelian Higgs Background 

We start with the expansion ()5.102|) obtained in |2Dj in the form 
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ln/«(A:) 



oo 

^ j drr $2^'^(r)V^(r)$^(r), 





x$o ^(r')V^(r')$!}(r'), 



(8.14) 



(8.15) 



ln/(3)(A;) 



7r\3 



2i 



In fi'\k) 





x$0,^(r')V^(r')$°,(r')$r(OV^(O<^'°j(r"), 



— I I dr r I dr' r I dr" r" I dr'" r 



2iJ 



X 



(8.16) 





(4$°^^(r)V(r)$°^(r) $°^^(r')V^(r')$|}(r ) 

$°,^(r")V^(r")$j(r")'^'j^(r'")^(^"')'^j(^"') 
+2$°,'^(r)V(r)$°,(r) <l>°,^(r')V^(r' )$?,(/) 

$° ^(r")\>(r")$°j(r")$°j^(r'")F(r'")$°j(r'")l , (8.17) 



As a matter of fact we have to treat three kinds of products, entering in the 
integrals 15.1261 namely 

^Vix)Vi^)^%i^), ^Vi^)Vix)^'j{x), ^Vi^)Vix)^'jix), 

and we denote them as G{hvh}{x), G{jvj}{x) and G{hvj}{x) for the first, second 
and third product respectively. 

Performing the products (|8.18p explicitly we obtain the expressions as follows: 



G{HVH} = fi{r)iC-H.H_ - C+H+H, 
G{JVJ} = fi{r)iC-J^J^ - C+J+J+) - 
G{HVJ} = Kr){C-H^J^ - C+H+J+) 



2v{r) 



VC+CF{H+H_) (8.18) 



2v{r) 



r 
f (r) 



VC+CF{J+JJ) 



(8.19) 



VC+C-{H+J^ + H^J^ 



].20) 



The notations has been used here are 
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=Pq± mo 



(8.21) 

Recall that we have also to construct the Jost function of imaginary momen- 
tum ik. This procedure leads to substitution of Bessel functions by modified 
ones. As a result we have for the products defined by (j8.2U|) : 

2X2 



G{HVH} Q{KVK} = ie ""^ [ — \ X 



ITT 



i^li[zp{K_K^ + K+K+) + mo{K+K+ - K_K_)] - '^^^ + mlK.K+^ 
G{JVJ} G{ivi} = (8.22) 
-z e---" !^fiMl-I- + /+/+) + mo(/+/+ - /-/-)] + + 
_ 2 

y{Hvj} y{Kvi} = — X 

i^fiMK.I. - K+U) - mo(K-/- + K+/+)] - ^^p2 + ml{K+I_ - K_I+)^ 

The notation of indices ± means as previous: 

/± = 4±i(gr) (8.23) 
K^^K^iqr) (8.24) 

The substitution k ^ ik implies also the following redefinition of qo,PQ 

k^ik (8.25) 
Po i p = i \/k'^ — ml (8.26) 



Qo ^ i q = i y p'^ + ml = i y k"^ — ml + ml (8.27) 

respectively. 

It can be seen from ()8.23|) that the Bessel functions which are to be integrated 
in ()8. 14118. l7|) are ordered in products of pairs of kind 

r K^K^, r I^I^, r K^I^, r K^K±, r I^I±, r K^I± (8.28) 
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Now we substitute the uniform asymptotic expansions represented as ()8.23p 
in the ()8. 14118. l7j) and rewrite each of integrands ()8.281 indices are omitted) as 



KI ^ tQki 
K K ^ r e-^'^^OKK 

(8.29) 

The apphcation of the expansion ()8.1.2|) to ()8.17p leads to 



oo 





^ ^ 

1 frg'l, + 2g',, 6{gii + rg'jj)r]" , 3rgij{2r]"Y 2rgnn 



+ 



z/3 I {2ri'f {2ri'Y {2r]'f {2r]') 

oo 

ln/(=^)~-2 (^)Y rdrg 



, 1 r'^gKign , 1 / 'ir'^gKiguri' 

KK \ w ,xo V -7 



y2 ^3 ^ 2r]' 

r'^gKig'ii - r^g'KiGii 

oo 


Here primes mean the derivation with respect to r, and following notations 
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are introduced: 

Gii{r) = 
t 



27rz/ 



ip 



3 5+[t,-l/2] 3 5+[^^i/2], 



+2?;(r)zoi^/(r ) exp(2^ 



(8.31) 



t 

inv 



3 sr[t, -1/2] + -1/2] 3 g-^^^ 1/2] + 1/2] 
ip I exp 2^ --^ exp2^ 



3 sr[t, -1/2] + -1/2] 3 g-^^^ 1/2] + 1/2] 
—m I exp ^ — : — h exp ^ 



-v{r)zQi'/{r'^) 



l + t ^ 



3 57[t,l/2]+5;[t,-l/2] 



l_t 3 5-[^^_i/2]+^+[i^i/2] 
exp 2^ 



QKK{r) = 

-(l)^lf".) 



l + t 



1 - t , 
T^exp(2^ 



3=1 



(8.32) 



57[i,-i/2], 



i=i 



+.(ji±^exp(2i::X£_l^)-Ji^exp(2i: 



57 [t, 1/2], 



1 -t 



2^ ,^57[t,l/2]+S7[t,-l/2], 
-2i;(r)zoz^/(r'')exp(2_^ ■' 



v3 



(8.33) 



Since the functions entering here depends not on r exphcit, but on tir) only, 
we use the derivation with respect to t instead of r, keeping in mind the follow- 
ing obvious relations for derivatives of t with respect to r (index r denotes the 
derivation) : 
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tr = --{l-t^) 

r 

trr = ^{l-t^){2-3t^) 

trrr = ^{l-ty{l5t^ -Q) (8.34) 

Then we can carry out the simple but very laborious procedure, hke the case 
of spinor field in the magnetic background. Namely we use the general form of 
the uniform asymptotics by means of 

g, z/ — i> cx), q/u = zq = const (8.35) 

of the logarithmic Jost function ()5.1()2j) . To this end we proceed from the expres- 
sions of \nf!)^\ik), J = 1,2,3,4 obtained in [IH] and given above in this section 
of Math.App. ()8. 14118. 17j) . They are already expressed in terms of products ()8.23|1 . 
One has now to insert the uniform asymptotics ()8.35p for each of the modified 
Bessel functions entering here ()8.24j) . 

q/u = Zq = const (8.36) 
4±i(gr)^4ii(z/Zor) (8.37) 
K,^i{qr) ^ K^^i{u zor) (8.38) 

which are listed explicit in Math.App 18.1.11 ()8.HI8.2j) . Furthermore, also the 
expansion 

P = ^ff + < = \fM^^o -'oi^-^l (8-39) 

is to substitute. The procedure has been executed with help of Mathematica alge- 
braic package. The final separation of powers u provides the following expressions 
for the first order of u 



In r'{tk) - —i \" ' v{ry - mo/i r tor' - ' V + t ""^^ ^ 

ly J r \ 2r^ 2 2 

+0(1) (8.40) 



Note that this expression must depend on the k- variable, the momentum 
"outside of the vortex core", related to r ^ oo, but it looks not explicit. To 
restore the dependence on k we use the corresponding expressions for 



rriQ — mi 



zo = Zo{z) = q/u = \ z^ + , (8.41) 



Z/2 



to = l/^JlVzy (8.42) 
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we can express it in terms of t = 1/ vT + ^V^ and expand for small 1/z/ , so we 
obtain the contribution to the next non zero order of l/u 

t„~t-tM^^ + 0(^) (8.43) 

Furthermore we omit the imaginary term in I8.4UI since it does not contribute 
to the vacuum energy (because of the symmetry properties of the Jost function 
mentioned in Sec l5.2|) . The terms proportional to odd powers of potentials u(r) 
and /i(r) are not present because of symmetry properties of exact Jost function 
mentioned above Sec J5.2l Finally we have for the order \jv: 



1 

u I r 



dr (t(t'^ — l) , ,a ru(r)n'(r) 
I ^ -v(rY + '^^ ^ ^ 



zH 



+ 0( 



(8.44) 



where /i is redefined as /i(r) = /i(r) — rrie. 

The analogous procedure of collecting all powers of l/z/^ in In fal^ (ik) . . . faf {ik 
with the help of Mathematica algebraic packages one arrives at the contribution 
of l/u^ order in In f°''^{ik): 



1 - ^t\l - t')] v{rf - lr\'{rf - 1(1 - 5t^).(r)^' 



^rh'il - M^)v{rf{pi{rf - ml) - ^h^Z - 5t'){^{rf - ml) - ^vH^f^'ir)'- 



1 



rH'{fi{rr-m-i) 



2\2 



2\2 



{mj - m g) 
4z4 



iA5) 



what leads to the representation in terms of Xij in (|5.135p . 



133 



8.2 Performing of the spectral sum for Bag boundary con- 
ditions 

Consider the expression ()5.4H) for regularized ground state energy Sec 15.11 We 
proceed like the case of scalar field and represent the function Bi~^^\k) in the 
form: 



1 + ) ^ — ^ (- ) (8.46) 



Then by applying of the In on Bl^^\k) and keeping in mind the asymptotics 
of Hankel functions ()5.21|) . we obtain in the limit of i? ^ oo 

ln^(+^)(fc) = ln/+(fc) + ^ 1^ ^ + ^ (^^ - 9 + - ^ - ^) + 

ln[i\/po~+^ + VPo ~ ^\ + + exponential suppressed term] (8.47) 

The last term tends to zero in the limit R ^ oo, under consideration, that k lies 
on the upper branch {k = ^k + ie) . Further, by applying of the derivation with 
respect to k, we get the following terms survived: 



|l.B<-)W ^ + + + (8,48) 



PO+m 



where the phase 9 is defined by ^ = arctg 

The second term is obviously proportional to the empty two-dimensional vol- 
ume. The contribution of this term in E^'^^ is now easy to calculate as 

/.^- r(.-l) ttR' (P + m^)^- 

0Fr(s-i) 8 s-2 ^ ' ' 

that in the large mass limit A; — is proportional to and corresponds to the 
contribution of the coefficient oq in heat kernel expansion ()4.119p . 

The third term is a pure imaginary value and can be omitted since the energy 
is defined to be real. Consider the last term in the large mass limit. Then the real 
part can be integrated immediately and it turns out to be proportional to and 
represents therefore the non-local contribution from the boundary, corresponding 
to the coefficient ai/2- In the considered case it cancels with the corresponding 
term arising from Si ^\k). 
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The crucial fact for this procedure is the vanishing dependence on the back- 
ground through the relation (^^ — | — i7) = 0. 

The performing of Bi~^\k) is analogous, and the final expression results in 

ISM- 
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8.3 Calculation of the asymptotic part of ground state 
energy 

8.3.1 Useful relations 

The Abel-Plana formulae 

The sum over half-integers v has been transformed to integrals according to 
the Abel-Plana formula as follows: 

oo oo 

The summation over integers I can be performed as 

oo oo 

/" „ 1 . ^ /" rf/ F(il) - F(-il) 

/—I ^ 

'-^ 

Integral relations 

The integrations over u and k can be done by using the identities: 

dv I dkik-m^Y = ^ — - — -4 ^— (8.52) 

' ^ ' dku^ 2 (rm)"-ir(|) ^ ' 



oo oo 



m 

and 



/ - m^r-y ^ -,„-^.£(izi)£|£±i^^^^ (8,53) 

8.3.2 Scalar field with scalar background 
Scalar field 

oo 

The sum ^ over / in ()5.173|1 can be replaced by the integration over 
1=1 

dl according to the Abel-Plana formula 1)8.511) . Thus the calculation of is 
reduced to the calculation of the following three parts: 
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£1' 



1 r(s — 1) sinvrs 
2v^ T{s-\) 



2 _2M-sH 



R2 00 2 ^ 





cas 
^2 



1 r(s — 1) sinvrs 

2v^ r(s-i) 



1.55) 



cas 
^3 



1 r(s — 1) sinvrs 

3 7 



-R2 00 

dk{e - m')^-^^ / r dr 



R2 37 

•i n=\ 7=1 



U=0 



2 J ' J 1 - e^^^' 

Ri 



X 



n=l j=l 



{ay 



-iiy 



^.56) 



To calculate the for each of r;,,j we use the identity fl8.52|l . and for 

the expression of type £^2* ^^e ()8.53j] respectively 



327r \ s m 



+ ln^-l)(i??-i?,^) + ^(i + lng-2)(i??-i?^, 



^.57) 



cas 
^2 













[64m 


V^2 





m 



(i?2 — Ri] 



327rm 



(i?2 - Ri) (8.58) 



The results for E^^ contain expressions of kind 



sin(7rs + — J 



1-e' 



2ls+i/2 



(8.59) 



which are not convergent in the limit s ^ for several n,j and must be fur- 
ther performed. The partial integration executed several times over / allows to 
abrogate divergent powers of [Z^ — (mr)^]. The similar procedure is described in 
details also in this section below, for the problem of spinor field inn NO-vortex 
background. 

All the remaining surface terms turn out to be zero. Finally the results 

in 



£^3^(1^0, R2) 



2tx 



fl - g/2 - ^/3 - 



\/2 

247rm2^^' 
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where 



and 



fi 

h 



mRi 

oo 



mRi 

oo 



oo 



mRi 
oo 



Fi{l)Gi{l,Ri)dl- 
F2{l)Gi{l,Ri)dl- 
F^{l)Gi{l,Ri)dl- 
FS)Gi{l,Ri)dl- 
F2{l)G2{l,Ri)j 



mRi 



2-kI ) 



m = ^(Kkt3^)')', 



m,i?2 
oo 



m/?2 

oo 



oo 



Fi{l)Gi{l,R2)dl 
F2{l)G,{l,R2)dl 
F;{l)Giil,R2)dl 
Fi{l)Gi{l,R2)dl 

F2{l)G2{l,R2)f, 



mi?2 



^2(0 = /( 



1 



(8.60) 



/(I -e2-') 



)' (8.61: 



11 /3 
^4(/) = /(t(t( 



Gi(/,x) =ln 



j_ 

mr 



Thus the finite part of -E"* is represented in terms of integrals not complete 
calculable analytically, but well convergent numerically. 



8.3.3 Spinor field with magnetic background 

The procedure of calculation of i?"* for spinor field is described completely 
in [201, the similar problem is considered for the spinor field in NO- vortex 
background, so we omit all details for this problem and give the final result only: 



oo 

-4 f dr, 



[n a{r) gi{rm) — Q r a'{r) g2{rm) + Q a{r) g3{rm)] 



n^eiiRi,R2) + n^e2{Ri,R2) 



(8.62) 
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here gi are 



9i{x) = - ^/i,3(a;) + ^/3,3 - ^/3,5W + y/3,7(a;) - ^/3,9(a;) 

11 7 5 

- 2^5x(--/3,3(a;) + -/3,5(a;) - -/3,7(a;)) 

+ l^dli^fsA^) - ^/3,5(x)), 

93{X) = -^(/3,3(x)-6/3,5(x) + 5/3,7(x)), (8.63) 

with fij are 



/l,3(x) 
/3,3(X) 



1 



1 + e^'^^ 



1 + e 



2lTX ■ 



1 + e 



2ttx ' 



3x1 + e^^^ 

The derivation of this resuh can be viewed as a particular case of more general 
problem, which is given in details below. 

8.3.4 Spinor field with NO vortex background 

During the calculation of the E"'^ one has to evaluate the following types of 
integrals over u: 

°° 3 3n 

1 3 Jm J 1 -^1 1 3 



"=1 i=l u= 



2 ' 2 



(8.65) 

Treated by means of the Abel-Plana formula ()8.5()|1 it can be written as 



= £f + £f (8.66) 
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where 



CCLS r)/~1 

0^2 — g 



dv 



iil + e 



(8.67) 



It is known from the experience of previous calculations that the i^f* con- 
tributes to the finite and the Si'' to the divergent part of S^. Proceeding from 
the explicit form of J^(z^) ()8.65p one left with the evaluation of the integrals of 
kind 



(8.68) 



in and 



-IV] 



.69) 



in 82'^ respectively. It comes about by using of the indentities ()8.52|8.5Hp . The 
results for ()8.68I8.69|) corresponding to ^ for each of tt, = 1, 3 and j = 1, 3, 5, 7, 9 
are: 



Til = Ii3 = -Cgm'm-''V{2 - s)V{s - 1) 

J33 = X35 = X37 = X39 = -^--^^C,m-6T(2 - s)T{s - 1), (8.70) 



for 8^^ and 



dv 



^n-4C.cos(..)^r(2-.)r(.-l)^ 



dv 



(z/2 - [mr]2)^-i/2^i'3('^) 
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„2s-2 



Jsa = 8CsCos{tts] 



-r(2-s)r(s + 



^2.-2^(^2 -s)r(s + 3/2) 



2' J (i/2 - [mr]2)--i/2 



1^37 = 



/3,5(^) 



32Cs cos(7rs) 
1/39 = 

64Cs cos(7rs) 



^2s-2r(2-s)r(s + 5/2) 



150F(3 + 2s)(2s + 1)(1 - 2s) J (z/2 - [mr]2)^-i/2 



dv 



.2s-2r(2 - s)r(s + 7/2) 



105v^(5 + 2s)(3 + 2s)(2s + 1)(1 - 2s) J (z/^ - [mr]2)«-V2 



du 



(8.71) 



for ^^2^ respectively. 
Here the relation 

[(^^)2 + 3,2]-(s+,/2)| 



x=—ii/ 



-2(mr)2^+^'[w2-(mr)2]-(^+-''/2) sin[7rs+j7r/2] (8.72) 

has been applied, and the partial integration several times has been executed to 
abrogate the singular denominator [(mr)^ + x'^]~^'^~^^^'^^ at s — > 0. 

The cancels against the E'^^^ (Sec J5.4p . and we arrive finally at the finite 
part of the asymptotic energy: 



E"' = £. 



as 
■2 



00 00 

— du X 
TT y r-^ J 

X 



!^v{rf - fi3{u) - i/33(z/) + ^/35(^^) " ^fsiil^) + ^/39(^^) 



f (r)' 



111 

^/33('^) - ^f35{l^) + Y^fsril^ 



+ r'^{^{rY — ml) 



1 1 
-^fssi^) — 2/35! 



fu{l^) + ^/35('^) - ^/37('^) 



+ 

(8.73) 



where x means mr and functions /ij(z^) are as usual 
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/37(^) = \l 



1/(1+62^") 

1 r 1 1 ^ 



1 + e 



2'KV 



1 + e 



2-KV 



1 


1 


'1 


1 


1- 


r 


1- 


V 


1/ 




_1 + e^'''' _ 
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8.4 The renormalization of the NO-vortex on-shell 

We start with the classical energy density of the NO-vortex per unit length for 
the particular case n = 1, Sec 01 ()5 .19811 



£ = 71 / r dr 



v'{r) +7] 



(8.74) 

Here the expression of S is given in terms of r- variable, (but not p = qif). 
There are various possibilities to rewrite the classical energy "on shell" using 
the equations of motion 



fir) + 



fir) Jir) 



n 



[l_^(^)]2 + Ar/2(l-/(r)2)/(r) = 0, 



v"(^r) - ^ + gV/(^)'[l - vir)] = 0. 



^.75) 



Let us transform the mixed term in the ()8.74|) : 

i(l-^(r))V(r)^; 



(8.76) 



from the second equation of ()8.75|) we have: 

1 



^(1 - .(rrf 



v"(r) — 



Mr) - 1). 



Now we rewrite the vortex equations in the form 



/(r) fir) + 



fir) 



+ Xr,\l-f{r)^)f{rf 



v"{r - 



v'{r) \ V — 1 f{ry^ 



(8.77) 



(l-t;(r))2, (8.78) 



so that right-hand sides of the both equations are equal and therefore the left- 
hand sides as well. Further, if we rewrite the ()8.74|1 using ()8.77j) 



8 = TT I r dr <^ -^-^ [t>'(r)^ + (t>(r) — l)(t>"(r) — v'{r)/r)~\ + rj^ f'{rY+ 



(fir)' - 1) - \v\f{r)' - 1) 



(8.79) 



and use a dummy trick (represent the coefficient at rf'f'{r)'^ as 1 + ,^ — ,^), we 
obtain: 
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8= n I r dr\-^ [v'irf + {v{r) - l){v"{r) - v\r)/r)\ + rfif{rf+ 



Xr] 



4 



(/(,)4 _ 1) _ Ar^4(/(r)^ - 1) + r^\l - 0/'(r)^ • (8.80) 



2 

Note, that the integral of full derivative 



dr^Hfiryrr (8.81) 

is zero, and \.{f {r) f {r))' = f'{rY + f{r) + Then it follows from 

(JH2H1), that 

firf = -fir) (^r(r) + = (8.82) 

W - ^) ^ + A,^(l - /(r)^)/(r)^ 
So the classical energy can be rewritten in the form: 

£ = n jrdr |er/V'(r)^ + (2^ - l)^(/(r)^ - 1) + iWil - f{rf) + 

1 



[v'iry + e(t;(r) - l)iv"ir) - i;'(r)/r)] ^ . (8.83) 



We remark, that the integral of the full derivative 



OO 



drjUm^.mX (8.84) 





provides a finite contribution —2a under assumption that v{r) — >■ ar^ + br + 
at r — > 0, and we find that 

Tc r dr { -7r^^iv{r) - l)iv"{r) - v'{r)/r) \ = -2aTC— - n rdr { — r— - 
J { q'^r'^ J Q J [ 

(8.85) 

Finally we arrive at the expression: 

8 = -2a7r-^ + (8.86) 



TT I rdr{ ir,'f{Tf + (2^ - l)^(/(r)^ - 1) + - f^rf) + ^^'(0' 
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Now we have a number of possibilities to define ^. We require for example: 



Then 



:vVir? - 1) + 4^;T^^^^'(r)^ 



2(A- 1) ' ' ' ' gV2 2(A 

Finally, we can give an interpretation of counter terms entering in S'^K To this 
end it is suitable to represent the divergent energy i?*^™ fl5.148j) . Sec 15.41 as 



E' 



div 



1 



f 4 4 

r dr^(/(r)2-l) + 



(8.89) 



1 , 4^2 

- + In - 2 



r dr 



4 3 



Now the subtraction S — E provides the ordering of coefficients as follows: 



cl 



E 



div 



1 

47r 



r dr { {f{rf - 1) 



lnl/£-2 



+ 



2(A 
47r2 A 



2 2A-1 



+ 



+ 



1 , 4/i2 

- + ln^ 



+ In 



4/i2 



2 



+ 



+ 



2(A-1) 
47r2ry2^ 



+ 



2(A 



Thus, the possible treatment of the renormalization scheme can be established 
as follows: If we state, that the first term, proportional to (/(r)2 — 1) corresponds 
to the finite redefinition of the coupling constant A, and the second term with 
t''(r)2 implies the finite redefinition of the U(l)-coupling (electric charge), then 
only one free parameter remains, namely the Higgs condensat ?], which must be 
fixed by two remaining terms. This fixing condition is then easy to obtain: 



2^2 



exp 



7r2A /4 



- 1 



^.90) 
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Here we recalled, that 



m = me = ^. (8.91) 
v2 

(as introduced in Secl^l)- 

This example shows, that one of the free parameter of this model must be fixed 
anyway, since we have four relations for three parameters at the given structure 
of the counter terms. 
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